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EFFECTS OF ROCK COATING ON REFLECTANCE SPECTRA OF ROCK SAMPLES. H Akiyama 1 and S. 
Sasaki 2 , ^ishimatsu Construction Cn , Ltd. 2-4-10-202 Fukami Nishi Yamato-city Kanagawa-pref. 242-0018 Japan 
akiyama@dream.big.jp, 2 Dept.Earth & Planetary Sci., Univ. Tokyo, Tokyo 113-0033 JAPAN sho@eps.s.u- 
tokyo.ac.jp. 


Introduction: Optical theory of “rock coating 

(dust film type)” is important, because rock surface on 
planetary bodies is frequently covered with another 
materials such as dust and regolith Rock coating is 
defined as the state that the surface ■ »f the rock is cov- 
ered by material of differing optical properties [1]. 

Hapke [2] applied a two-layer model to the rock 
coating. However, validity of the two-layer model has 
not been discussed enough. Hapke compared meas- 
ured values [3] with his theoretical values only at a 
specific wavelength (750nm). Then: is no guarantee 
that measured value agrees well with the theoretical 
value in other wavelengths. We tried a different 
method to match the measured value for wide range of 
wavelength. 

The “adding method” is used in the studies of vege- 
tation canopy and of solar emission In this method, 
multiple scattering between coating layer and rock sur- 
face is calculated and reflectance of coated rock sur- 
face is estimated. W T e here applied the adding method 
to the rock coating. 

Experimental Method: As for dust film-type rock 
coating, it is important to measure the optical proper- 
ties of the film and the reflectance of the coated sur- 
face. In this paper, olivine powders were used as a 
material of dust film. We scattered powders thinly and 
homogeneously on the surface of the rock, and simu- 
lated rock coating (Fig.-l). Thickness of the olivine 
layer at each survey point was measured by a laser dis- 
placement gage, and the reflectance of each survey 
point was measured by the bidirectional spectrometer. 

Moreover, the reflectance and transmittance of 
the olivine layers were measured using an integrating 
sphere, since the reflectance and transmittance of the 
covered layer should be known in the adding method. 



Albite 

Olivine 

Anorthite 

£nstatite 

Oiopsitfe 



Fig.-l Mineral palette (simulate rock surface) and Rock 
Coated surface 


Results and Discussion: We compared the meas- 
ured values and theoretical values by Hapke' s. Com- 
parisons between the measured values and both theo- 
retical values for the optical depth are shown in Fig.-2. 
When optical depth is large, the measured values agree 
well with both theoretical values. However, when opti- 
cal depth is small, the Hapke’s values gives large dis- 
agreement approaching at most 70%. On the other 
hand, disagreements of the new theoretical values re- 
main less than 10%. 

In Hapke two-layer model, the optical properties of 
each coating layer are estimated implicitly assuming 
the layer thickness is sufficiently large for multiple 
scattering. If the layer is optically thin, the theory can- 
not be applied. Multiple scattering is not expected 
when the layer's optical depth is less than 1. As for 
Fig. -2, when optical depth is smaller than 1, Hapke 
theoretical reflectance values disagrees from the meas- 
urement value greatly. 
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Fig.-2 Disagreement of measured reflectance values 
from theoretical values. 


Doubling Method: Using a so-called doubling 

method, we can calculate reflectance and transmittance 
of any thickness of layer if we have reflectance and 
transmittance data of one optically-thin layer. How- 
ever, it is not clear that doubling method can be applied 
to the rock coating. We performed measurements to 
confirm the applicability of the doubling method. Us- 
ing optical properties of spherical glass beads and ir- 
regularly-shaped olivine powder, theoretical values of 
the film for thickness of approximately 0.4mm was 
calculated from the measured values of 0.2mm thick- 
ness. Fig.-3 shows the result. The disagreement of the 
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measured values and the theoretical values are less than 
20% when the particle size is under 90pm as for about 
0.5mm thickness layer. The disagreement is larger than 
40% in the particle size over 125pm. 


reflectance and emittance spectroscopy 1 ', University Press, 
Cambridge. [3] Wells, E. et al. (1984) ICARUS, 58, 
331*338. [4] Douglas L. Smith (1997) Engineering & 
Science, 3, 8-17. 


l aw of doublia0 MttK>d<r*fl«ct«noa| 



Fig. -3 Disagreement of reflectance between theoretical 
0.4mm thickness values by the doubling method and meas- 
ured 0.4mm thickness values. 

Application: There are various applications of the 
new rock coating theory. MARS Pathfinder, which 
landed and surveyed MARS surface, observed two 
types of rock coating. First, IMP (Imager for MARS 
Pathfinder) observed colored rock surface. Some rocks, 
like “Yogi”, they have two colors; half side of them is 
bluish gray and half side of them is red. Bluish gray 
side is original rock surface color, and red is the color 
of the oxidized iron contained in the circumference 
regolith [4]. Using these surfaces spectral data and 
Fe 2 03 spectral data, we estimated the thickness of this 
coating layer. It can be estimated that that thickness is 
about 0.08mm from this experiment if the material, 
which covers Barnacle Bill, is supposed to be the 
Fe 2 0 3 particles of 0.1pm size. 

Second, MARS Pathfinder had two radiative cali- 
bration targets (RCTs). Their reflectance were chang- 
ing because airborne dust coated their surfaces. IMP 
observed change of target reflectance. We used these 
data, and estimated the thickness of the air-borne dust- 
coating layer. As for this estimate, we used Hematite 
reflectance to airborne dust reflectance. Martian Sol 10 
data agrees well with 0.13mm thickness theoretical 
data, and Sol 20 and Sol 35 data agree well with 
0.25mm thickness theoretical data. And Sol 52 data 
agrees well with 0.50 thickness theoretical data. 

However actually, the covered particle size and 
their chemical compositions are different. This discus- 
sion is only a preliminary examination. 


References: [l]R.Dorn (1998) “Rock Coating”, 
Elsevier, Amsterdam. [2] B. Hapke (1993) “Theory of 
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On the Distribution and Implications of Mantled and Exhumed Terrains on Mars. R. E. Arvidson 1 , P. R. Chris- 
tensen 2 , J. Garvin 3 , M. Mellon 4 , department of Earth and Planetary Sciences, McDonnell Center for the Space Sci- 
ences, Washington University, St. Louis, MO, 63130 (arvidson@wunder.wustl.edu), 2 Arizona State University, 
Phoenix, AZ, 3 NASA Goddard Spac e Flight Center, Greenbelt, MD, 4 University of Colorado, Boulder, CO. 


Introduction: The intent of tins paper is to show 
the relationships for Mars among albedo, thermal iner- 
tia, roughness inferred from MOLA pulse width spread 
data, and geology inferred from photogeological analy- 
ses. Mapping of surface units using these parameters 
and approaches, in combination with analysis of hyper- 
spectral image data from ISM, TES, OMEGA, and 
CRISM observations, will maximize our understanding 
of the distribution and nature of surface units on the red 
planet. Results will directly impact the selection of 
landing sites that exhibit geological records needed to 
understand planetary habitability. 

Data Sets and Approach: Our approach is to 
compile and examine a global-scale image cube con- 
sisting of MOC WA Geodesy Campaign mosaics [1], 
MOLA-based elevation [2], MOLA pulse width data 
that have been corrected for regional scale tilts (CPW) 
[3], bolometric albedo (over solar wavelengths) from 
TES observations [4], and thermal inertias modeled 
from TES emission data [5]. These observations were 
examined together (e.g., Fig. 1) to infer surface proper- 
ties and to map geologic units, emphasizing the nature 
and extent of mantling of preexisting terrains. The 
work is preparatory to establishing units for analysis of 
TES spectral emission data, ISM spectral reflectance 
data, and upcoming OMEGA and CRISM observa- 
tions. In addition, crater counts were used to help un- 
derstand the relationships between surface properties, 
morphology, and the extent of mantling. 

Results: Previous work has shown that there is a 
negative correlation on a global basis between albedo 
and thermal inertia [e.g., 5], In particular the Tharsis 
Plateau, Elysium region, and Arabia Terra have high 
albedoes and low thermal inertias, consistent with dust 
mandes. The classical dark areas have high inertias 
and there is an intermediate surface that has been inter- 
preted as due to exposures of weakly indurated soil 
deposits of the type documented at the Viking landing 
sites. 

We add a number of new' points to this type of 
analysis: 

• The dust and duricrust surfaces have elevated 
CPW values relative to classical dark areas. Fur- 
ther, examination of MOC WA and NA images for 
the dust covered areas in Arabia Terra shows that 
the deposits are tens of meters in thickness in most 
cases, thinning to approximately ten meters at the 


edges. Crater size frequency data for the mantled 
areas are consistent with this interpretation and 
with a steady accumulation of dust at a rate of a 
fraction of a micrometer per year. Arabia has been 
accumulating dust for billions of years. Similar 
trends are emerging for dust-covered areas on the 
Tharsis Plateau and within Elysium. 

• A number of units on the Tharsis Plateau and 
within the Elysium region exhibit very high CPW 
values along with intermediate albedoes and low 
inertias. These materials clearly correspond to 
wind eroded ash deposits. 

• Dark, high thermal inertia terrains have low CPW 
values as compared to the dust mantled areas. This 
is an interesting pattern that presumably is associ- 
ated with reworking of surfaces by aeolian activity. 
We offer two possibilities: local reworking of dust 
into bedforms and surfaces or beveling and grad- 
ing of the dark areas by wind. 

• High inertia and CPW surfaces can also be found 
locally within craters, ejecta deposits, and within 
fracture and valley systems, particularly in Valles 
Marineris. 

• Exhumed terrains can be found at the margins of 
the dust mandes, e.g., in Terra Meridian!. In this 
region a widespread set of exhumed layered de- 
posits can be found, and each major horizon 
mapped from MOC WA and NA images exhibits 
unique albedo, color, inertia, CPW, and spectral 
emissivity properties. This includes the hematite 
deposit mapped from TES data [6] and the under- 
lying etched unit that has morphologic patterns 
similar to the young platy ridged pahoehoe flow's 
mapped in the Athabasca region. On the other 
hand, dust mandes show uniform spectral signa- 
tures [4]. 

Implications: We restrict implications to those as- 
sociated with targeting OMEGA and CRISM data for 
high spectral and spatial resolution mapping. This tar- 
geting is for spectral analyses of the sites to understand 
planetary habitability and to help select landing sites 
for surface exploration of terrains that would provide 
strong evidence of past climatic or hydrothermal condi- 
dons and thus habitability. The bottom line is that the 
dust mantles should be avoided since these regions 
have accumulated up to tens of meters of aeolian cover. 
Careful mapping of exhumed or non-mantled terrains is 
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needed to find places that preserve and expose evi- 
dence of ancient aqueous conditions. These targets 
will form the core of the list for detailed hyperspectral 
mapping. 

References: [1] Caplinger and Malin, (2001) JGR , 
1 06, 23595-23606. [2] Smith et al., (2001), JGR , 106, 


23689-23722. [3] Garvin et ah, (1999), GRL, 26, 381- 
384. [4] Christensen et ah (2001), JGR, 106, 23823- 
23871. [5] Mellon, M. et ah, (2001), Icarus , 148, 437- 
455. [6] Christensen et ah (2001), JGR, 106, 23873- 
23886. 




Figure 1: Perspective view (top) and nadir view (bottom) covering all longitudes and from 50° N to 60° S in latitude. 
Red color is modulated by TES-based bolometric albedo [4], green by MOLA topography-corrected pulse widths 
[2], and blue by TES-based thermal inertia [5]. Colors have been overlain on a shaded relied map generated from 
MOLA elevation data (2). The prime meridian is located in the middle of the frame. Relief exaggeration for the per- 
spective view is 200:1. Orange areas correspond to regional-scale dust deposits that are bright, have low inertias, 
and are rough. These include Arabia Terra, the triangularly-shaped area in the middle, top part of the map. Yellow 
areas correspond to wind-eroded volcanic ash deposits that have intermediate albedoes, low inertias, and are very 
rough. Blue areas are the classical dark areas with high inertias and low roughness values. Purple regions have in- 
termediate brightness and thermal inertia, roughness values similar to regional-scale dust deposits, and have been 
interpreted as consisting of duricrust. 
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TESTING THE HAPKE MODEL FOR THE OPTICAL EFFECTS OF SPACE WEATHERING. D. T. Ble- 
wett 1 and P. G. Lucey 2 , 'NovaSol, 1 iOO Alakea St., 23rd Floor, Honolulu, HI 96813 (dave@nova-sol.com), 2 Hawaii 
Institute of Geophyiscs & Planetology, Univ. of Hawaii, 2525 Correa Rd., Honolulu, HI 96822 (lu- 
cey @higp.hawaii.edu). 


Introduction: In the Apollo er.i, Hapke and col- 
leagues used simulation experiments and theoretical 
considerations to argue for the importance of vapor- 
phase coatings in causing lunar soils to darken and 
redden with age [1]. During the 1990s, direct detection 
of vapor-deposited iron on actual lunar soil grains 
came from the work of Keller and associates [2]. 
These workers discovered that most soil grains have 
coatings, some of which differ from the interior com- 
position. The morphology of the coatings is such that 
they must have condensed from a vapor. The coatings 
and grain rims typically contain abundant submicro- 
scopic iron particles (SMFe). While SMFe is largely 
responsible for the reddening of lunar soil, Hapke’s 
work suggests that it alone does not account for the 
changes accompanying soil maturity. For example, 
Hapke [3] found that the inclusion of dark glass was 
required in addition to submicroscopic iron to account 
for weathering trends. Using a different formalism, 
Shkuratov et al. [41 emphasized the importance of 
known particle size evolution in maluration modeling 
to match observed space weathering I rends. Therefore 
Hapke's ground-breaking work must be followed up to 
show how the model glass content and composition and 
grain size evolution required to fit observed data com- 
pare to measured soil properties. 

Spectral Modeling: We have developed computer 
code to produce model spectra and model fits to labo- 
ratory or telescopic spectra using the equations in [3]. 
The program uses 1 1 parameters: the abundances of 
olivine, orthopyroxene, clinopyroxcne, anorthite, il- 
menite, glass, large-grained metal, and SMFe, the 
grainsize, the Mg#, and the FeO content of the glass. 
Figure 1 shows a telescopic near- infrared reflectance 
spectrum of Aratus, a 10-km diameler highland crater 
located in the Hadley- Apennine region of the Moon 
[5]. The test fit shown in the plot for Aratus was done 
with the Mg#, grainsize, and glass FeO content held 
constant. The fit parameters were of 11%; opx, 17%; 
cpx, 4%; an, 52%; ilm, 2%; glass, 14%, large-grain Fe, 
0%, and SMFe, 0.08%. This produces an excellent fit 
to the original spectrum. We then t >ok this composi- 
tion and held it constant except for the SMFe content, 
which we varied in equal steps from 0.005 to 0.23 
wt.%. This produces a simulated space- weathering 
trend. The results are shown in Figure 1. The model 
spectrum with very little SMFe (top spectrum in Figure 
1) is similar to spectra of powdered lunar rocks, with 


strong absorption bands and a flat continuum. As opti- 
cal maturation proceeds, the overall albedo drops, the 
continuum reddens, and the absorption bands are at- 
tenuated. 

We have undertaken a study to use the Hapke 
model, including grain size and glass content evolution 
(representing accumulating agglutinate glass), to fit the 
Apollo 16 soil maturity series of Fischer and Pieters [6] 
and compare the model results to the known correlation 
of glass content, particle size and SMFe content. The 
Apollo 16 samples have nearly constant composition 
and mineralogy, but differ in maturity, i.e., they have 
variations in grainsize and the abundance agglutinitic 
glass and SMFe. Bidirectional reflectance spectra for 
Apollo 16 samples 60051, 61141, 61161, 61221, 
61241, and 67461 were retrieved from the RELAB 
web site at Brown University [7]. The laboratory spec- 
tra and model Fits are shown in Figure 2. As expected, 
the amount of SMFe in the model fits increases with 
the maturity of the samples. 

We plan to compare the model predictions of glass 
content, particle size and SMFe content to the actual meas- 
ured values. Ultimately our goal is to couple these 
parameters in a systematic and physically realistic way, 
and apply these systematics to other highland and mare 
soils in the RELAB, telescopic and Clementine data- 
base. We also seek to understand the compositional 
controls on the OMAT optical maturity parameter [8] 
that has been used to remotely map relative variations 
in soil maturity. 

Acknowiegement: DTB thanks Bruce Hapke for 
serving as a mentor and teacher early in his career. 

References: [1] Cassidy, W. and B. Hapke, Icarus 
25, 371-383, 1975; Hapke, B., Moon 7, 342-355, 
1973; Hapke, B., Phys. Earth Planet. Inter. 15, 264- 
274, 1977; Hapke, B. et al.. Moon 13, 339-354, 1975. 
[2] Keller, L. P. and D. S. McKay, Science 261, 1305- 
1307, 1993; Keller, L. P. and D. S. McKay, Geochim. 
Cosmochim . Acta 61, 2331-2340, 1997; Keller, L. P. 
et al., in Workshop on New Views of the Moon , p. 44- 
45, 1998. Wentworth, S. J. et al., Meteorit. Planet. Sci. 
34, 593-603, 1999. [3] Hapke, B., J . Geophys. Res . 
106, 10,039-10,074, 2001. [4] Shkuratov, Y. G. et al. 
Icarus 137, 222-234, 1999. [5] Blewett, D. T. and B. 
R. Hawke, Meteorit. Planet. Sci. 36, no. 5, 701-730, 
2001. [6] Fischer, E. M. and C. M. Pieters, Icarus 111, 
475-488, 1994; Fischer, E. M. and C. M. Pieters, J. 
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Geophys. Res. 101 , no. El, 2225-2234, 1996. 
[7] http://lf314rlds.geo.brown.edu/. [8] Lucey, P. G. et 
al. f / Geophys . Res. 105 , no. E8, 20,377-20,386, 2000; 
Grier, J. A. et al., 7. Geophys . Res. 106 , no. El 2, 
32,847-32,862,2001. 



Figure 1 . Telescopic spectrum (points with error 
bars) and model spectra (solid lines) for the crater Ara- 
tus. From top to bottom, the model spectra represent a 
space-weathering trend produced by increasing the 
amount of SMFe, with other parameters held constant. 
Telescopic spectrum is scaled to the Clementine 0.75 
pm albedo of the area covered by the spectrometer 
aperture. Model plots were done in the range 0.75 - 
2.0 pm because of limitations in lab data for pure min- 
erals in the UV, and uncertainties in the optical con- 
stants of Fe at long wavelengths. 



Wavelength (microns) 

Figure 2. Reflectance spectra for six Apollo 16 soil 
samples with similar compositions but differing expo- 
sure degree (maturity) are shown with the symbol (+). 
Model fits for each spectrum computed from Hapke's 
equations are shown as solid lines. 
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PHOTOMETRIC MODELING OF A COMETARY NUCLEUS: TAKING HAPKE MODELING 
TO THE LIMIT. B. J. Buratti ( 1PL), M. D. Hicks (JPL), L. Soderblom (USGS), J. Hillier (Grays 
Harbor College), D. Britt (U. Term.) 


In the past two decades, photometric models 
developed by Bruce Hapke have been fit to a 
wide range of bodies in the Solar System: The 
Moon, Mercury, several asteroids, and many icy 
and rocky satellites. These models have enabled 
comparative descriptions of the physical 
attributes of planetary surfaces, including 
macroscopic roughness, particle size and size- 
distribution, the single scattering albedo, and the 
compaction state of the optically active portion 
of the regolith. One challenging type of body to 
observe and model, a cometary nucleus, awaited 
the first space based mission to obtain images 
unobscured by coma. 

The NASA-JPL Deep Space 1 Mission (DSI) 
encountered the short-period Jupiter-family 
comet 19/P Borrelly on September 22, 2001, 
about 8 days after perihelion. Prior to its closest 
approach of 2171 km, the remote -sensing 
package on the spacecraft obtained 25 CCD 
images of the comet, representing the first close- 
up, unobscured view of a comet's nucleus (1). 

At closest approach, corresponding to a 
resolution of 47 meters per pixel, the intensity 
of the coma was less than 1% of that of the 
nucleus. An unprecented range of high solar 
phase angles (52-89 degrees), viewing 
geometries that are in general attainable only 
when a comet is active, enabled the first 
quantitative and disk-resolved modeling of 
surface photometric physical parameters. 

The geometric albedo of Borrelly is 0.029 ± 
0.006, comparable to the dark hemisphere of 


Iapetus, the lowest albedo C-type asteroids, and 
the Uranian rings. The Bond albedo, 0.009 ± 
0.002, is lower than that of any Solar System 
object measured. Such a low value may 
enhance the heating of the nucleus and 
sublimation of volatiles, which in turn causes 
the albedo to decrease even further. A map of 
normal reflectance of Borrelly shows variations 
far greater than those seen on asteroids. The 
two main terrain types, smooth and mottled (2), 
exhibit mean normal reflectances of 0.03 and 
0 . 022 . 

Bonelly's physical photometric parameters 
are typical of other small dark bodies. The 
nucleus exhibits significant variations in 
macroscopic roughness, with the oldest, darkest 
terrain being slightly smoother. This result 
suggests the infilling of low-lying areas with 
dust and particles that have not been able to 
leave the comet. The surface of the comet is 
backscattering, but there are substantial 
variations in the single particle phase function. 
One region scatters isotropically (like cometary 
dust), suggesting that its regolith is controlled 
by native dust rather than by meteoritic 
bombardment. 

Work performed at NASA’s Jet Propulsion 
Laboratory, California Institute of Technology. 

(1) Soderblom, L. et al. (2002). Science 296, 
1087. 

(2) Britt, D. et al. (2002). LPSC Abstract 1686. 
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THE DIVERSE SURFACE COMPOSITIONS OF THE GALILEAN SATELLITES. R. W. Carlson, Jet Pro- 
pulsion Laboratory, California Institute of Technology, Pasadena, CA 91 109 (rcarlson@Hvely.jpl.nasa.gov). 


The galilean satellites represent a diverse collection, ranging 
from the volcanic moon Io, with a surluce that is changing 
yearly, to Callisto, with a dark, ancient surface overlying ice. 
The composition of these surfaces are also quite different due 
to a variety of processes and influences, including tidal heat- 
ing, radiolysis, gardening, a magnetic field (Ganymede K and 
meteontic infall. 

lo’s surface contains large quantities of sulfur dioxide (S0 2 ) 
and colorful sulfur allotropes, both originating in plumes and 
flows from the tidally driven volcanoes. A broad, 1-pm band 
is found at high latitudes and may be due to absorption by 
long-chain sulfur polymers produced by SOi radiolysis, al- 
though iron and iron sulfide compounds are candidates. An 
unidentified 3.15 pm absorber is equatonally distributed 
while a 4.62 pm band, perhaps due to a sulfate compound, 
exhibits a non-uniform distribution. Hot spots are generally 
dark, and some exhibit negative reflectance slopes (i.e., 
blue). The composition of these lavas has not been estab- 
lished spectroscopically, but the high temperatures of some 
volcanoes suggest ultramafic silicates or perhaps more re- 
fractory material such as oxides. 

Europa’s icy surface contains the radiolytic product hydro- 
gen peroxide (H 2 O 2 ). This molecule is rapidly destroyed by 
solar ultraviolet radiation, so replenishment of H 2 O 2 by high- 
energy charged particle impact must be occurring continu- 
ously. The same process may produce annospheric and surfi- 
cial 0 2 and the hydrated compound that was found on Eu- 
ropa, suggested to be hydrated sulfuric acid. The hydrate is 
observed to be associated with visually dark material, whose 
spectra resemble sulfur polymers, and with surficial S0 2 . A 
radiolytic sulfur cycle continuously creates and destroys 
these compounds, and their observed relative abundance 
agrees with predictions based on laboratory radiolysis ex- 
periments. The equilibration time for this cycle is ~ 100 
years. Sources of sulfurous material include Iogenic plasma 
implantation, preferentially deposited on the trailing hemi- 
sphere where the hydrated acid, visually dark material, and 


S0 2 are indeed found to concentrated. Thermal processes 
such as diapiric heating, accompanied by H?0 sublimation, 
can further concentrate the sulfur compounds and produce 
the variegated surface patterns, while gardening and asyn- 
chronous rotation modify the global patterns. An alternative 
suggestion for Europa’s hydrate is hydrated sulfate salts from 
an internal ocean, but the radiolytic lifetime of these com- 
pounds at the surface is < 4000 years so they would be rap- 
idly assimilated into the H 2 S0 4 , S0 2 , and S x cycle. NIMS 
spectra are best matched with sulfuric acid hydrate. Surficial 
carbon dioxide is also observed on the leading hemisphere of 
Europa. This C0 2 may arise from the radiolysis of carbona- 
ceous meteoritic material that is deposited preferentially on 
the leading hemisphere or it may be an outgassing product. 

Ganymede and Callisto are both covered by a dark de- 
posit overlying ice, the latter exposed on ridges, slopes, and 
within craters. Both satellites show similar spectroscopic 
features, generally more pronounced in Callisto’ s spectra. 
C0 2 is observed to be concentrated on Callisto’s trailing 
hemisphere, suggestive of magnetospheric influence, in par- 
ticular radiolysis. A corresponding C0 2 atmosphere is found 
on Callisto. The atmosphere is rapidly depleted by ionization 
and sweeping, so there must be a continuous source of C0 2 
such as radiolysis or outgassing. A band at 3.88 pm has been 
suggested to be an SH stretch band, but radiolytically- 
produced H 2 C0 3 provides an equally good fit. A 4.57-pm 
band has been suggested to be from a CN stretch transition 
but the shape and position is better fit by carbon suboxide 
(C 3 0 2 ), a radiolysis product of C0 2 . Sulfur dioxide has been 
suggested as a surface constituent on both Ganymede and 
Callisto based on a 4-pm infrared feature. However, the cor- 
responding feature is not observed on Europa even though an 
S0 2 ultraviolet feature is present in Europa’s trailing side 
spectra. Ganymede and Callisto’ s 4-pm feature may repre- 
sent something other than SO 2 
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EVOLVING PERSPECTIVES ON SPACE WEATHERING OF ASTEROIDS. Clark R Chapman, 
Southwest Research Institute, Suite 400, 1050 Walnut St., Boulder CO 80302. 


Efforts to measure the colors of asteroids 
date back at least as far as the 1930s, when Bo- 
brovnikoffs spectroscopy actually yielded 
Vesta’s rotation period based on cyclical 
changes in its reflectance spectrum. There were 
several efforts in the 1940s and 1 950s to com- 
pare asteroid colors with laboratory measure- 
ments of meteorite colors. They proved incon- 
clusive for several reasons, including (a) the 
blandness of mineral spectra in the short- 
wavelength bands commonly used in astronomi- 
cal photometry and (b) poor preparation of often 
terrestrially weathered meteorites. 

In the early 1970s, both laboratory measure- 
ments of properly prepared meteorite samples 
and astronomical observations of visible and 
near-IR asteroid spectra had advanced and were 
sufficient to address the fundamental questions 
concerning asteroid mineralogy and relationship 
to meteorites. Two basic approaches were em- 
ployed: (a) attempts to assess asteroid mineral- 
ogy from first principles, based on absorption 
band centers for instance, and (b) spectral 
matching, in which it was simply assumed that if 
an asteroid spectrum resembled the spectrum of 
a particular class of meteorite it could be in- 
ferred that the two had similar mineralogy, even 
if the specific spectral features were not under- 
stood or were not unique. 

Application of both approaches yielded a few 
successes, particularly in associating Vesta with 
basaltic achondritic meteorites, but they gener- 
ally revealed an unexpected disparity: despite 
sharing some spectral similarities, the most 
common asteroid type (not corrected for obser- 
vational selection) could not be reconciled with 
the most common meteorite type (again, not cor- 
rected for atmospheric filtering processes). Thus 
emerged the conundrum in which S-type aster- 
oids apparently yielded few meteorites while 
ordinary chondrites (OCs) had essentially no 
observable parent bodies in the asteroid belt. 

Several prominent cosmochemists assumed 
that the S-types must be associated with OCs 
and that the disparities must be due to unknown 
technical issues associated wiih the remote- 


sensing techniques or interpretations. But aster- 
oid specialists (including myself) reached a 
nearly unanimous consensus from the late 70s 
through the early 90s that it was not “scientific” 
to adopt unknown space weathering processes to 
resolve the discrepancy. They turned instead to 
adopting equally unknown, very selective dy- 
namical processes (which had little basis in the 
published literature of asteroid physics concern- 
ing dynamical or collisional processes) to derive 
OCs from small, rare, generally unseen asteroids 
while excluding meteorite delivery from the 
abundant S- type asteroids, which even domi- 
nated the Near-Earth Asteroid (NEA) popula- 
tion. This attitude of asteroid specialists was 
very frustrating for cosmochemists since the de- 
gree of space weathering required to resolve the 
conundrum was in fact modest compared with 
the extraordinary (if not specifically understood) 
space weathering effects that clearly were re- 
sponsible for the dramatic differences between 
telescopic reflectance spectra of the lunar sur- 
face compared with laboratory measurements of 
lunar samples. 

Once again, despite the in-your-face evidence 
for lunar space weathering, most specialists re- 
fused to accept any extrapolation to the asteroids 
until the actual physical processes responsible 
for lunar space weathering were identified. In 
the mid-90s, several types of asteroid observa- 
tions (e.g., Galileo studies of Gaspra and Ida and 
spectral studies of small NEAs) demonstrated 
that space weathering processes (of some un- 
known sort) must be occurring on asteroids. Si- 
multaneously, laboratory studies finally identi- 
fied the predominant causes for lunar space 
weathering as well as showed that micrometeor- 
ite “zapping” of asteroidal minerals can modify 
their spectra in the sense that S-types differ from 
OCs. Finally, in the last few years, advances in 
asteroid dynamics (primarily due to enhanced 
capabilities of numerical simulations) have 
demonstrated that the combination of minor 
resonances and the Yarkovsky Effect necessarily 
leads to general, representative sampling (rather 
than highly selective sampling) from most of the 
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inner half of the main asteroid belt. So there can 
no longer be any doubt about the association of 
OCs with many S-type asteroids and the reality 
of space weathering processes (probably associ- 
ated with the solar wind and micrometeorite 
bombardment) as the dominant cause of the dis- 
crepancies. 

In retrospect, one wonders why it took more 
than a quarter of a century for the early advo- 
cacy of space weathering (by Bruce Hapke, 
among others) to gain wide acceptance. Some- 
how, the demand for explicit physical “proof’ of 
specific space weathering processes (even as 
NASA refused to fund the required research on 
space weathering) delayed for three decades the 
acceptance of what we now realize is correct. 
Such attitudes can now be seen to have been 
intellectually unsupportable: the asteroid spec- 
troscopists generally refused to take seriously 
the evidence from fields outside their discipline, 
specifically physicists, whose evolving under- 
standing of collisional and dynamical processes 
never supported the extreme selectivities on 
which the asteroid remote-sensing community 
relied. 

Asteroid science was similarly led astray for 
a few years (in the late 60s and early 70s), by 
analogous attitudes on the part of dynamicists, 
such as George Wetherill, who refused to be- 
lieve that any meteorites could be derived from 
the asteroid belt because — in that pre -chaos 
epoch — no specific dynamical processes had 
yet been identified that could bring asteroidal 
samples into Earth-crossing orbits. They criti- 


cized cosmochemist Ed Anders for not taking 
dynamics seriously when Anders argued that the 
meteorites must somehow come from the aster- 
oid belt even if the dynamicists had not yet fig- 
ured out how. 

In both cases, the respective specialists in- 
sisted that it wasn’t “scientific” to rely on un- 
known processes until the specific processes 
were proven to operate. Even as recently as last 
year I heard loud objections to the term “space 
weathering” because of the speaker’s belief that 
we do not yet know for sure, from first princi- 
ples, exactly how space weathering works. Yet, 
in both cases, history has demonstrated that the 
non-specialists who basically argued that any- 
thing was possible in the spheres of other fields 
so long as it wasn’t specifically disproven turned 
out to be right. I would not extrapolate these 
cases to every scientific dispute and suggest that 
hand-waving generalists or non-specialists who 
ignore conservative specialists outside their 
fields of expertise will always prove to be right. 
But I think there are lessons to be learned from 
these case studies in asteroid science. Sometimes 
intuitively plausible hypotheses may prove to be 
true — provided they don’t violate fundamental 
physical laws — even if we haven’t yet figured 
out specifically how they could be true. We still 
operate in considerable ignorance and have 
much to learn, as I think Beth Clark 1 s following 
talk will exemplify. I expect that studies of as- 
teroid space weathering will still hold some sur- 
prises for us. 
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THE NATURE OF SPACE WEATHERING EFFECTS ON ASTEROID 433 EROS. Beth Ellen Clark, Physics Department ; 
Ithaca College , Ithaca NY 14850, USA, tbdark@ithaca.edu ), Paul Lucey, University of Hawaii, HIGP, Honolulu HI 96822, USA, 
Paul Hel fenstei n, Astronomy Department, Cornell University, Ithaca NY 14853, USA. 




Figure 1: Near-Earth Asteroid Rendezvous (NEAR) observa- 
tions of 433 Eros have revealed evidence for space weathering 
on an asteroid surface. This oblique view across the top of 
Psyche crater on Eros shows the steep northern slopes. Steep 
slopes allow gravitational processes to compete successfully 
with surface darkening processes. 


Combined Multi -Spectral Imager (MSI) (0.95/im) and 
Near-Infrared Spectrometer (NIS) (0.8 to 2.4/im) observations 
of Psyche crater on 433 Eros are analysed to investigate the 
nature of space weathering effects on Eros. At 5.3 km in diam- 
eter, Psyche is one of the largest craters on Eros which exhibit 
distinctive brightness patterns consistent with downslope mo- 
tion of dark regolith material overlying a substrate of brighter 
material (Figure 1). At spatial scales of 620 m/spectrum. Psy- 
che crater wall materials exhibit albedo contrasts of 32-40% 
at 0.946/jm. Associated spectral variations occur at a much 
lower level of 4-8% (i2 — 4%). 

In seeking to constrain the effects of space weathering pro- 
cesses on Eros, we have investigated several possible causes 
for the observed albedo and spectral trends: grain size differ- 
ences, olivine, pyroxene, and troilite variations, and lunar-like 
optical surface maturation. 

Spectroscopic data used for this analysis have been pho- 
tometrically corrected [1]. Photometric correction procedures 
correct each spectrum for variations due to illumination and 
viewing geometry. The spectra we used were constrained to 
phase angles less than 100° , incidence and emission angles 
less than 75° , and spatial resolution 400 - 800 m/spectrum. 

We use mineral mixture models to simulate the observed 
spectra. Next, we modify our mineral mixtures to simulate 
the spectral consequences of variations in grain size, compo- 
sition, and optical maturity. We constrain our model mixture 
variations to be consistent with the observed albedo contrasts 
of 32-40%. Finally, we compare model predicted spectral 
variation trends with the measured variations. 

We use an intimate mixing model based on the work of 


Figure 2: (left) The nominal mixture model spectrum of bright 
materials (top solid line spectrum) in Psyche is compared with 
mixture models made by increasing the proportions of both the 
neutral dark and the SMFe endmembers. (right) Each compar- 
ison spectrum is ratioed to the nominal spectrum and normal- 
ized to 1.0 at 1 3pm. Also plotted are the observed dark/bright 
material spectral ratios (solid dots) to aid in comparison. This 
figure shows that the albedo contrasts and overall amplitude of 
spectral variation trends observed in Psyche crater can be sim- 
ulated with a mixture model incorporating enhanced neutral 
dark and SMFe endmember components. 


Hapke to quantitatively explore the effects of input param- 
eters. Hapke [2] presents the methods necessary to com- 
pute the optical effects of submicroscopic iron which serves 
an important role in lunar space weathering optical effects. 
Our model allows the computation of the reflectance spec- 
trum of a meteorite-like assemblage consisting of intimately 
mixed olivine, ortho- and clino-pyroxene, plagioclase feldspar, 
troilite, Fe-Ni metal, and submicroscopic Fe metal. We can 
compute the spectrum of an arbitrary mixture of these compo- 
nents at any grain size. 

The model begins with a definition of the relative abun- 
dances of the components and the mineral chemistries of the 
components For the transparent minerals (olivine, pyroxene, 
plagioclase and glass) we compute the optical constants from 
defined chemistries. Chemistry -dependent optical constant co- 
efficients for olivine and pyroxene are from Lucey [3]. The 
chemistries of olivine and orthopyroxene were linked to con- 
form to the relationship between olivine and pyroxene Mg- 
number in ordinary chondrites from Keii and Fredrikkson [4]. 
Optical constants for glass as a function of iron content were 
computed by Lucey [3]. We used the Hapke-based methods 
presented in Lucey [3] to compute the iron dependent optical 
constants of plagioclase from reflectance spectra of plagioclase 
in the USGS Denver spectrum library. 

Next, using the methods outlined in Hapke [5, 6, 7, 8] mod- 
ified to include the effects of submicroscopic iron presented 
by Hapke [2] we compute the single scattering albedos of the 



14 Solar System Remote Sensing 


SPACE WEATHERING ON EROS: Beth Ellen Clark et al. 


transparent components from the optical constants, grain sizes 
and the defined abundance of submicroscopic iron, assuming a 
coating of metal for the minerals, and a uniform distribution of 
metal throughout the grain for glass. We note that lunar and as- 
teroidal submicroscopic iron deposition may be quantitatively 
different because the asteroidal submicroscopic iron may also 
contain nickel. Our models assume that the submicroscopic 
coatings are pure iron, which appears to be a good assump- 
tion on the Moon Hapke [2]. The single scattering albedos 
are combined using equation (17) of Hapke [5]. The single 
panicle phase functions are combined using the methods of 
Helfenstein et al. [9,10]. The mixture single scattering albedo 
is converted to reflectance using equation (37) of Hapke [2J. 

We find that the albedo contrasts in Psyche crater are not 
consistent with a cause due solely to variations in grain size, 
olivine, pyroxene or lunar-like optical maturation. A grain 
size change sufficient to explain the observed albedo contrasts 
would result in strong color variations that are not observed. 
Olivine and pyroxene variations would produce strong band- 
correlated variations that are not observed. A simple lunar-like 
optical maturation effect would produce strong reddening that 
is not observed. The contrasts and associated spectral variation 
trends are most consistent with a combination of enhanced 
troilite (a dark spectrally neutral component simulating optical 
effects of shock) and lunar-like optical maturation (Figure 2). 

These results suggest that space weathering processes may 
affect the spectral properties of Eros materials, causing surface 
exposures to differ optically from subsurface bedrock. How- 
ever, there are significant spectral differences between Eros’ 
proposed analog meteorites (ordinary chondrites and/or prim- 
itive achondrites), and Eros’ freshest exposures of subsurface 
bright materials. 

Perhaps most surprisingly, after accounting for all differ- 
ences in the measurement units of our reflectance compar- 
isons, we have found that the bright materials on Eros have 
reflectance values at 0.946/im consistent with meteorites, but 
spectral conti nua that are much redder than meteorites from 
1.5 to 2.4/jm. Most importantly, we calculate that average 
Eros surface materials are 30-40% darker than meteorites. 

In summary, we find 

(a) The background surface of Eros in and around Psyche 
crater is 32-40% darker than the bright materials exposed on 
steep slopes in the crater walls at spatial scales of approxi- 
mately 620 m/spectrum. 

(b) The dark material is 2-4% redder than the bright ma- 
terial at 1.5-2.4/im, and has a 1 -micron band which is 2-4% 
broader at 0.8- 1 .0/im. The total spectral variation for the high- 
est albedo contrast area found on Eros is 4-8% peak-to-peak 
from the normalized mean at 1 .3pm. 

(c) The albedo and contrasts in Psyche crater are not con- 
sistent with a cause due solely to either grain size or composi- 
tional variations. 

(d) The contrasts and associated spectral variation trends 
observed in Psyche are consistent with a lunar-like optical mat- 
uration which combines accumulation of submicroscopic iron 
with an enhancement of the dark spectrally neutral component 


of Eros’ surface materials (presumably shocked silicates). 

(e) Bright materials in Psyche appear optically less mature 
than dark materials, but do not spectrally match meteorites. 
The bright materials are approximately the same brightness as 
meteorites, however their spectra are 33% redder from 1.5 to 
2.4pm. 

(f) We calculate that average Eros surface materials are 30- 
40% darker than meteorites. We also find that when differences 
in reflectance measurement units and photometric geometry 
are accounted for, most meteorites are 2-3 times brighter than 
asteroids. 

(g) If Eros is compositionally similar to ordinary chon- 
drites, and if the darker materials on Eros are due to space 
weathering (optical maturation of a fresher brighter substrate), 
then the weathering process must be a strange two-step mech- 
anism whereby materials are reddened before darkened (bright 
materials are redder than meteorites), and darkened more than 
reddened (albedo contrasts are stronger than color contrasts). 

(h) At the scale of our observations Psyche crater did not 
excavate materials which are spectrally similar to meteorites. 
Nor does Psyche have any bright crater ejecta deposits, indi- 
cating that since the formation of the crater the surface of Eros 
has been weathered. 

(i) These findings suggest that Eros’s surface may be com- 
pletely weathered due to its exposure to the space environ- 
ment. Our results suggest that space weathering processes 
affect Eros’ spectral properties, but do not fully explain the 
differences between Eros and its proposed analog meteorites, 
the ordinary chondrites and/or the primitive achondrites. 

(j) Differences in optical maturation trends between Eros 
and other S-type asteroids studied by spacecraft may indicate 
different rates of the competing processes of surface maturation 
and impact cratering between objects in Near-Earth orbits and 
objects in the main asteroid belt. 
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Introduction: ANTS (Autonomous Nano- 
Technology Swarm), a mission architecture 
consisting of a large (1000 member) swarm of 
picoclass (1 kg) totally autonomous space- 
craft with both adaptable and evolvable heu- 
ristic systems, is being developed as a NASA 
advanced mission concept [1, 2], and is here 
examined as a paradigm for lunar surface ex- 
ploration. 

As the capacity and complexity of hard- 
ware and software, demands for bandwidth, 
and the sophistication of goals for lunar and 
planetary exploration have increased, greater 
cost constraints have led to fewer resources 
and thus, the need to operate spacecraft with 
less frequent human contact At present, 
autonomous operation of spacecraft systems 
allows great capability of spacecraft to ‘safe’ 
themselves and survive when conditions 
threaten spacecraft safety. To further develop 
spacecraft capability, NASA is at the fore- 
front of development of new mission archi- 
tectures which involve the use of Intelligent 
Software Agents (ISAs), performing experi- 
ments in space and on the ground to advance 
deliberative and collaborative, autonomous 
control techniques. Selected missions in cur- 
rent planning stages require small groups of 
spacecraft weighing tens, instead of hundreds, 
of kilograms to cooperate at a tactical level to 
select and schedule measurements to be made 
by appropriate instruments onboard. Such 
missions will be characterizing rapidly un- 
folding real-time events on a routine basis. 
The next level of development, which we are 
considering here, is in the use of autonomous 
systems at the strategic level, to explore the 
remote terranes, potentially involving large 
surveys or detailed reconnaissance. 


One Previous Application: The use of 

this approach has already been discussed [1,2] 
for mainbelt asteroids, a population consisting 
of more than a million small bodies. Although 
a large fraction of solar system objects are 
asteroids, relatively little data is available for 
them because the vast majority of them are 
too small to be observed for more than single 
point measurements except in close prox- 
imity. Determination of the systematic 
distribution of physical and compositional 
properties within the asteroid population is 
crucial in the understanding of the solar 
system formation. In addition, there has been 
interest in asteroids as sources of exploitable 
resources. Far more reconnaissance is 
required before either goal is achieved. A 
traditional mission approach (to explore 
larger asteroids sequentially) is not adequate 
for studying the asteroid population 
systematically. Our approach involves the use 
not of a smart spacecraft with ‘drones’, but of 
a totally autonomous truly distributed 
intelligent network of sensors. Each sensor 
will have specialized instrument capability 
(e.g., advanced computing, imaging, 
spectrometry, etc.) and heuristic systems that 
are both adaptable and evolvable, to 
simultaneously evaluate a broad cross-section 
of the entire population. Spacecraft use solar 
sails, which require no expendables, as the 
propulsion systems needed to fly from the 
Earth directly to asteroids 1 kilometer or 
greater in diameter. Clusters of spacecraft 
operate simultaneously within a broadly 
defined framework of goals to select targets 
from among available candidates and to de- 
velop scenarios for studying targets, and se- 
lected members of the swarm return to Earth 
with data and are replaced as needed. 
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Application to Lunar Exploration: The 

ANTS concept has broad application to 
planetary targets in the inner solar system 
where detailed reconnaissance is a next logi- 
cal step in exploration, but the cost and capa- 
bilities of lander or rover technology are still 
limiting factors. The Moon would be an ideal 
candidate for such a totally autonomous in- 
vestigative approach using stationkeeping or 
low- periapsis solar sailors, each with a spe- 
cialized instrument or capability and placed 
optimally for its viewing requirements. The 
many lunar remote sensing investigations, 
limited in coverage or instrumentation or 
both, which have already have been done 
provide an essential context for the step pro- 
posed here. A comprehensive simultaneous 
multi-sensor survey for evidence of the de- 
posits of ice or volcanic material inferred to 
be buried or just beneath the surface, could 
provide essential clues on the formation and 
history of, for example, polar regions, antipo- 
dal magnetic swirls, KREEP, the farside ba- 
sins, to name a few candidate sets of features. 
Or, clusters of spacecraft could be used to do 
a systematic highest possible resolution sur- 
vey of the entire lunar surface. Sensors to be 
considered should include visible and near 
infrared red imager/spectrometers, X-ray, 


Gamma-ray and neutron spectrometers, al- 
ti meter/ ranger/ sounders, magnetometers . 

Spacecraft could work in same sensor vertical 
or planar arrays for 2D or 3D imaging or field 
determination, for example, using magne- 
tometers to investigate antipodal magnetic 
swirls. Alternatively, multisensor clusters 
would allow optimal data fusion [3] when 
prospecting for a particular material, such as 
neutron, gamma-ray, and X-ray spectrome- 
ters, along with radio sounders, all with rele- 
vant compositional information to allow in- 
terpretation of measurements in terms of the 
presence and distribution of ice. The lunar 
environment is more familiar (to humans), 
and Earth/Moon communication is near real 
time. As a result, the Moon presents a lower 
risk than other targets and would be an ideal 
place for ANTS to train and for the totally 
autonomous distributed intelligence concept, 
and human/machine interfaces generally, to 
undergo development and testing without the 
expense of sending humans into space. 
References: References: [1] Curtis et al, 
(2000) IAU Proceedings, IAF-00-Q.5.08. [2] 
Rilee, Clark, and Curtis, (2002) IAU Proceed- 
ings, in preparation. [3] Clark and McFadden 
(2000) JGR, 105, 4291. 
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Introduction: Reflectance spectroscopy is now a tool 
widely used for remote sensing applications, both in 
planetary and terrestrial observation Among the large 
number of possibilities offered by this technique, the 
investigation of the surface mineralogy of a target is 
often explored. However, with the considerable 
increase of the spatial resolution of optical experiments 
on and near the planetary surfaces, it becomes more 
and more important to understand the physics that 
control the light interaction process with soils, in order 
to better interpret the physical (nature, texture, surface 
roughness, maturity, degree of cristallinity) and 
mineralogical properties of the surface at subpixel 
scale. 

Over the years, a number of semi -empirical models 
have been developed for analyzing the bidirectional 
reflectance of particulate surfaces on the basis of the 
scattering and absorption properties of minerals and 
rocks. In particular, Hapke [1J developed an 
approximate analytical solution, which provides a 
description of the photometric behavior in terms of 
physically meaningful parameters. 

This work uses a method, founded on a genetic 
algorithm, for a determination of the global set of 
parameters involved in Hapke’ s model for planetary 
surface analogs when dealing with a set of angular 
conditions representative of the usual range of 
observation in planetary exploration for spacebome 
optical instruments. 

Experimental measurements: Based on experimental 
measurements produced by a new laboratory wide- 
field multispectral imaging facility [2], the method is 
applied to three soils: two of them are light materials 
(single scattering albedo about 0.7) and one is dark 
material (single scattering albedo about 0.3). Those 
soils are sorted to four grain sizes, and are 
representative of the spectral variability of 
observations realized on planetary surfaces. 

The grain size sorting is the following: less than 75 pm 
(labeled “Gl”), from 75 to 250 pm (“G2”), from 250 
to 500 pm (“G3”) and from 500 pm to 2 mm (“G4”). 

The domain of geometries of obseivation investigated 
is controlled by the observational situations the most 
frequently met in planetary exploration by spacebome 


optical instruments (e.g. MOC/MGS, HRSC/Mars- 
Express, AMIE/Smart-1 , LIC/Lunar-A, LISM/Selene) 
The incidence angle is selected between 0° and 55° 
and the emergence angle between -70° to 0° (azimuth 
angle 0°) and 0° to 70° (azimuth angle 180°) with a 
minimum phase angle of 20°. 

In this range, ten measurements are produced in order 
to span the photometric variability. 

Results: The Hapke’s formalism using the parameter 
values determined by the genetic algorithm well 
models the variations in bidirectional reflectance 
related to particle grain size distribution for the three 
materials. It establishes that all the steps in our 
methodology are sound, including experimental 
measurements, the derivation of photometric quantities 
by means of Hapke’s model and the use of genetic 
algorithm for an optimized determination of the Hapke 
parameters [3]. 
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Figure 1: comparison between modelled and measured spectra 
for one of the twelve samples. 

Furthermore a by-product of this application is the 
derivation of Hapke parameters for samples 
representative of materials found on planetary soils, 
which may be used in future reflectance spectroscopy 
data interpretations: 

Phase function parameters (double Henyey-Greenstem 
function, [4]): 

• b keeps similar values for all samples, between 
0.4 and 0.5. 

• c generally increases with grain size. 
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Those parameters are well constrained, showing that 
the angular domain available in this study (phase angle 
between 20° and 120°) is large enough to determine 
this function. 

Backscattering function parameters ; 

• B 0 values are very small, for our dark material 
(the backscattering effect is negligible) and for 
our light materials, B 0 values increase with grain 
sizes from 0 for G1 to 0.6 for G4. 

This shows that materials with larger grains have a 
structure containing more tunnels, creating a stronger 
backscattering effect. This result is due to our target 
experimental design where the surfaces of samples 
with small grain sizes have a smoother aspect than 
those with larger grain sizes, implying that, when 
grains become smaller, the textural roughness of the 
observed surface decreases. 

• h value ranges about 0.15 for the Palagonite 
sample and about 0.55 for the Tephra sample, 
when the backscattering effect is not negligible 
(Bo greater than 0.2). 

Given the limited mesoscale topography range (see 
below), those values could not be interpreted as it is 
usually done, because there is an interaction with 
topographic shadowing function (see in this volume 
the companion abstract by Stankevich et al.). 
Backscattering function is less constrained than phase 
function. In fact to better describe backscattering 
effect, smaller phase angle measurements are needed. 
Roughness function : 

• 0 values are highly correlated with grain sizes for 
our light material samples: 14° for Gl, 20° for 
G2, 23° for G3 and 26° for G4. Darker material is 
linked with smaller 0 values. 

In the case of real planetary observations, we may 
consider three levels of scale (cf. table I): first, the 
largest one is the topographic scale representing the 
surface topography variations greater than the pixel 
size, which ranges typically between 5 and 100 m 
scale for current and existing measurements (e.g. 
MOC/MGS, HRS C/Mars -Express). Then, the smallest 
scale describes the particulate medium and is 
controlled by the grain size distribution, which goes 
from a few tens of microns up to a few hundreds 
microns. Intermediate to these two scales, exists the 
mesoscale structure that is linked to the topographic 
textural distribution inside each pixel, and is 
represented in Hapke's theory by the 0 parameter, 
referred to as the surface roughness. 

In this study, the pixel size is about 3 cm and grain size 
range goes from a few tens of microns (for Gl) to 2 
mm (for G4). As a result, the mesoscale domain is 
reduced in scale (from mm to a few cm) in our 
investigation. 



Particulate 

Medium 

Grain Size 

Textural 

Mesoscale 

Topography 

Observation 
Scale 
Pixel Size 

Planetary 

Observations 

1 0 - 200 fim 

1 mm - 1 m 

5- 100 m 

Laboratory 

Measurements 

10 - 200 |im 

1 mm - 1 cm 

3 cm 


Table I: Three considered levels of scale 
for planetary observation and laboratory 


The values found by the genetic algorithm are in the 
same range as those determined in McEwen study 
(1991) for planetary soils: it establishes that the 
algorithm gives plausible values despite the change in 
scale of observation implying a more reduced range 
for the mesoscale features (cf. table I). Indeed, our 
experiment shows that in case of large grain sizes in 
the millimeter range, the textural topography variation 
inside the pixels is dominated by the “rocky" aspect of 
the surface, this is not yet the case in planetary 
observation. However with the increase of the spatial 
resolution in orbital imaging or for in-situ 
measurements, this situation could arise. 

Conclusion: Hapke’s model implicitly considers that 
the scale of surface topography is much larger than the 
size of regolith particles handled as a particulate 
medium. However, in laboratory and in-situ planetary 
measurements one deals with surface topographies that 
have the characteristic scale comparable with the size 
of regolith particles, and often such mesoscale 
topographies are formed by rock fragments. This study 
shows that despite the change of scale, the Hapke’s 
formalism may describe the photometric behavior of 
the surface, but the meaning of 0 and h parameters 
becomes more difficult to handle with the possibility 
of ambiguous physical interpretation. Theses results 
point at the need for an improved theoretical modeling 
for describing the real case of natural regolithic 
surfaces, for which the mesoscale organization of the 
textural relief influences the surface photometric 
behavior (see in this volume the companion abstract by 
Stankevich et al ). 

Acknowledgement. 
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To interpret remote spectral observations, scattering and 
absorption in a particulate surface are simulated via radiative 
transfer models. The standard model for this purpose among 
the planetary science community is the Hapke model (e.g., 
Hapke 1981, 1993, 2002). This model (like many others) 
uses two parameters to characterize the optical behavior of 
individual grains in a particulate surface, the single-scattering 
albedo w and phase function p(g). These terms describe, re- 
spectively, the quantity and the angular distribution of light 
scattered by an individual grain. Unfortunately, these parame- 
ters are strictly optical. They can be rather difficult to interpret 
in terms of more interesting particle properties such as grain 
sizes, shapes, and compositions, that a remote sensing experi- 
ment might seek to discover. 

An equivalent slab approximation is i ypically used to relate 
w to the grain size and optical constants of the material. This 
approach can mimic the wavelength-dependent absorption be- 
havior of irregular grains, as long as the imaginary index k is 
much less than 1, the shape is equant, and the grain size D is 
much larger than the wavelength A. Unfortunately, the equiva- 
lent slab approach provides no information about p{g), which 
also has a strong dependence on optical constants and particle 
form (e.g., Grundy et al. 2000). Instead, p(<?) is commonly 
approximated with simple analytical functions such as Legen- 
dre polynomials or Henyey-Greenstein functions (e.g., Hapke 
1993), effectively discarding the compositional and textural 
information carried by the wavelength dependence of p(g). 
Often p(g) is assumed to be independent of wavelength (e.g., 
Cruikshanlc et al. 1998,Quincoefa/. 1999, Dout6e/ al 1999), 
thereby forcing the single- scattering albedo to account for all 
observed spectral features. Such an assumption can lead to 
significant errors, especially in simulating layered media (e.g., 
Johnson & Grundy 2001). 

There are many ways of computing or empirically mea- 
suring p(g) for an irregular grain having specified optical con- 
stants and form (e.g., microwave analogs, Gustafson 2000; 
discrete dipole approximations, Draine 2000; Monte Carlo 
ray-tracing, Macke 2000; etc.). Unfortunately, none of these 
techniques is fast enough for use in iterative radiative transfer 
models. Using equivalent spheres simulated via Mie theory 
offers the advantage of providing w and p(^) simultaneously, 
but such a model also has difficulties. It is still too compu- 
tationally expensive, and its application to grains in contact 
with one another (as in a planetary surface) is somewhat inap- 
propriate, since diffraction is included without considering the 
effect of neighboring grains. 

A new, much less computationally expensive, approxi- 
mate model is urgently needed for calculating w and p(<?) 
from optica] constants and particle form within iterative re- 
mote sensing radiative transfer models. A new single particle 
model would be especially valuable in newer models which 
are more accurate in their accounting for the angular distribu- 


tion of multiply-scattered light (e.g., Dout£ & Schmitt 1998, 
Cheng & Domingue 2000, Hapke 2002). 

We will present a status report on our ongoing effort to 
remedy this situation by developing a new model of w and 
p(g) for large, irregular grains. Previous studies have estab- 
lished that w and p(<j) exhibit consistent trends as functions 
of particle characteristics (e.g., McGuire & Hapke 1995, Hart- 
man & Domingue 1998). We are capitalizing on this regular 
behavior by using it to constrain simple models. We began 
by performing a Monte Carlo ray-tracing study of scattering 
by large, irregular particles (Grundy et al. 2000). Figure 1 
shows results from this study, in which the Henyey-Greenstein 
parameters of p{g) are seen to vary with absorption qD , re- 
fractive index n, and particle surface roughness. Our studies 
with the Monte Carlo model continue, with the results being 
used to constrain and test candidate models of w and p(p) 
as functions of grain size, form, and optical constants over a 
broad range of parameter space for large irregular particles. 



For each Monte Carlo run, we record distinct populations 
of rays separately; i.e., those reflected from the grain’s exte- 
rior, those transmitted once through the grain, those transmitted 
with one internal reflection, those transmitted with two inter- 
nal reflections, etc. In addition to the angular distribution for 
scattering from each sub-population, we record the mean in- 
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temal path lengths, polarizations, etc. Each of these scattered 
populations, as shown in Fig. 2, has a characteristic amplitude 
and angular distribution, which is directly dependent on par- 
ticle properties such as optical constants, shape, and surface 
texture. 



Our objective is to represent scattering by each of these 
separate populations using simple functional forms such as 
Henyey-Greenstein functions or Legendre polynomials. We 
are attempting to develop semi-analytical/semi-empirical ex- 
pressions for the amplitudes and angular distributions of these 
distinct contributions to p(g) , which can then be summed to ob- 
tain computationally inexpensive synthetic p(^) for arbitrary 
irregular grain parameters. From the estimated amplitude of 
each contribution, we can also compute w from the mean in- 
terior path lengths for the different ray classes, combined with 
the absorption coefficients and the Beer-Lambert law. 

It remains to be determined how many separate ray classes 
must be considered to achieve acceptable levels of accuracy. 
For weakly-absorbing particles (especially those with high re- 
fractive indices), the contribution from rays having been inter- 


nally reflected numerous times can be significant. Fortunately, 
this highly -scattered component of p(g) tends to have a rela- 
tively uniform angular distribution. We are investigating the 
possibility of merging all the highly- scattered components into 
a single isotropic or back-scattering term, having an amplitude 
estimated from the interior scattering efficiency Si and parti- 
cle opacity aD in a manner analogous to the equivalent slab 
approximation. 

Acknowledgements: Support for this project is being 
provided by NASA Planetary Geology and Geophysics grant 
NAG5- 10159 and by the Programme National de Plan^tologie 
delTNSU. 
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Introduction; To interpret the opposition phe- 
nomena observed for regolith surfaces optical proper- 
ties of individual particles should be studied. Even 
under neglecting diffraction, the problem of light scat- 
ter by particles of irregular shape is rather complex. 
Using the computer ray-tracing method, see, e.g, [1,2], 
we investigate here brightness and polarization back- 
scatter effects of particles with different shapes in the 
geometric optics approximation. 

Model description: We use a computer Monte 
Carlo ray tracing approach [2] that allows us calcula- 
tions of the particle scattering indicatrix with account- 
ing for polarization. A particle is generated in the com- 
puter memory. It is presented by a succession of trian- 
gular flat facets, for which the laws of reflection and 
refraction are applicable. Thus one can put a normal to 
each facet of any particle and calculate all necessary 
local angles. Then, the particle is "illuminated” with 
great enough number of rays. These rays are muluply 
reflected and refracted on their way Besides, they can 
be partly absorbed within the particle, if it is semi- 
transparent. The interaction with the particle surface is 
described by means of the Fresnel and Snell laws. 

In this work we concentrate on the comparison of 
optics of particles of perfect shapes widely used in 
various models of scattering to data for their non- 
perfect analogs and randomly shaped particles. Several 
classes of irregular particles are examined: (1) faceted 
spheres that have a different number of flat facets, N, 
which varies in our studies from 100 to 30000; (2) 
cubes of “spoiled” forms, when the length of cube sides 
slightly varies; (3) randomly shaped particles generated 
by an auxiliary random gaussian field. The latter is 
different from so called gaussian particles described in 
[1], where they are assumed to be stochastically de- 
formed spheres, the radius of which follows to the log- 
normal statistic law. In order to generate the particles 
of this class we use an auxiliary 3-D random gaussian 
field (RGF). After generating such a field, its dissec- 
tion at a certain level is made that gives a 3-D medium 
consisting of panicles with random shapes and sizes 
that could be called RGF particles [3]. 

We have thoroughly tested the model. In particu- 
lar, the algorithm was verified on the reciprocity prin- 
ciple expressed in the vector form. The quality of aver- 
aging-out over particle orientations has been checked 
by the criterion of smallness of the elements F l3 , F j4 , 
F 23 h F 24 and by the condition of symmetry F n = F 2 \ 
and anti-symmetry F 34 = -F 43 . Although here we pre- 


sent only intensity and polarization curves, our model 
allows calculation of all non-zero elements of scatter- 
ing matrix [3]. 

Results: The element F n and ratio -F 12 /F H of scat- 
tering matrix vs. phase angle a for faceted spheres hav- 
ing different numbers of facets and RGF particles are 
shown in Fig. 1 at n = 1.5 (refractive index) and k = 0 
(absorption coefficient). Results were averaged over 
orientations and shapes (in case of RGF particles). A 
sphencal on average particle formed with facets, a 
number of which is more than approximately 10000, 
can be considered as a rather perfect sphere. When N is 
lesser, great changes are observed. Thus, for F u at 



0 30 60 90 

Phase angle, deg 


Fig. 1 . Phase functions of brightness and linear po- 
larization degree for spherical particles approximated 
by different number of facets N and for RGF particles 
at n = 1.5 and k = 0. 
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N= 1000 the glory surge disappears and, moreover, 
instead of the surge an opposition “anti-spike” is de- 
veloped. At N = 100, neither the glory nor the first 
rainbow are observed. Weakening of the features with 
the decrease of N is clearly seen for the ratio -F n!F n - 
This ratio corresponds to the definition of linear polari- 
zation degree of scattered radiation, when particles are 
illuminated with unpolarized light. So, as one can see 
in Fig. 1, the small negative polarization branch, which 
is observed for the perfect spheres near the backscatter 
direction disappears for roughly faceted particles. The 
surge of positive polarization corresponding to the first 
rainbow goes down quickly with decreasing N. 

The phase function for RGF particles is almost con- 
stant in the concerned range of phase angles. It does 
not show any presence of the rainbow or opposition 
spike. Light scattered by this kind of particles is posi- 
tively polarized in backward direction and the polariza- 
tion degree is much less than that for spheres. 

In Fig. 2 the results for perfect and randomly de- 
formed cubes are represented. To obtain particles with 
such shapes, we used perfect cubes deforming them in 
the following way. Each vertex of a cube is randomly 
moved in a random direction. The shift has the zero 
average and the certain amplitude 5 that is given in the 
length of the cube edge. We varied & from 0 (perfect 
cube) to 0.2 (very irregular particle). As in the case of 
other particles the irregular cubes are approximated 
with triangle facets. The ideal cube reveals the bright- 
ness spike as well as strong negative polarization 
branch at small a. These results are in a good agree- 
ment with [4]. The backscatter spike is caused by the 
component of scattered radiation corresponding to 
double internal reflection and is quickly degraded with 
deviation of particles from cubical shape. This compo- 
nent also contributes very strongly to the negative po- 
larization that is rapidly relaxed with growing 6. 

Conclusions: Scattering properties of faceted 
spheres with number of the facets less than 10000 dis- 
tinct significantly from those of the perfect spherical 
particles. The phase angle curves of the element F n 
and ratio -F^/Fn (and other scattering matrix ele- 
ments too) of roughly faceted spherical particles, e.g., 
with N - 100, diverge very much from those of perfect 
sphere. The irregular particles show smoother curves. 
Even small deviations from the perfect shape kill the 
opposition spikes and negative polarization branch. 
The ideal cube exhibits the brightness spike as well as 
strong negative polarization at small a. When the 
cubes are deformed the spike are reduced and become 
wider; the negative polarization branch vanishes. 



Fig. 2. Phase functions of brightness and linear po- 
larization degree for particles with shapes randomly 
deviating from cube at n = 1.5 and k = 0. 

Our modeling shows that there is qualitative simi- 
larity of the phase angle curves of intensity F u and 
ratios -Fi 2 /F n (as well as of F 34 IF n , F 22 /F u , F 33 /F n , 
FxJFyi ) for the roughly faceted spheres (ellipsoids), 
RGF particles, and irregular cubes with 6 = 0.2. Thus 
strongly irregular particles of all studied classes reveal 
much more resemblance in scattering properties, than 
their perfect representatives. 

References: [1] Peltoniemi J. et al. (1989) Appl 
Opt 28, 4088-4095. [21 Stankevich D. et al. JQSRT \ 
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3051-3060. 
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Introduction: The Goldstone Solar System Radar 
(GSSR) has successfully collected radar echo data from 
Mars over the past 30 years. As such, the GSSR has 
played a role as a specific mission element within Mars 
exploration. The older data provided local elevation 
information for Mars, along with radar scattering informa- 
tion with global resolution (e.g. [1,2]). Since the upgrade 
to the 70-m Deep Space Network (DSN) antenna at 
Goldstone completed in 1986, Mars data has been 
collected during all but the 1 997 Ma *s opposition. Radar 
data, and non-imaging delay-Doppler data in particular, 
requires significant data processing to extract elevation, 
reflectivity and roughness of the reilecting surface [3]. 
The spatial resolution of these experiments is typically 
some 20 km in longitude by some 150 km in latitude The 
interpretation of these parameters while limited by the 
complexities of electromagnetic scattering, do provide 
information directly relevant to geophysical and geo- 
morphic analyses of Mars. 

Landing Site Assessment with Radar Data. The use- 
fulness of radar data for Mars exploration has been 
demonstrated in the past. Radar data were critical in 
assessing the Viking Lander 1 site [4, 5] as well as, more 
recently, the Pathfinder landing site [6, 7]. In general, 
radar data have not been available to the Mars explora- 
tion community at large. A project funded initially by the 
Mars Exploration Directorate Science Office at the Jet 
Propulsion Laboratory (JPL), and later funded by NA- 
SA’s Mars Data Analysis Program has reprocessed to a 
common format a decade’s worth of raw GSSR Mars 
delay-Doppler data in aid of landing site characterization 
forthe Mars Program. These data will soon be submitted 
to the Planetary Data System (PDS). The radar data used 
were obtained between 1988 and 1995 by the GSSR, and 
comprise some 63 delay-Doppler radar tracks. Of these, 
15 have yet to be recovered from old 9-track tapes, and 
some of the data may be permanently lost. The available 
data used for Mars Exploration Rover landing site 
analysis for example, are listed in Table 1. 

Enhanced Radar Observations: Two modem radar 
techniques offer the opportunity to improve the spatial 
resolution and the global coverage of available radar 


data :random-long-code delay- Doppler observations, and 
interferometric delay-Doppler observations. 

Random-long-code delay-Doppler radar. Instandard 
delay-Dopplerobservations, Marsisoverspread[8] This 
can be overcome by using a pseudo-random, non- 
repeating code that is much longer than the observing 
run. This technique was successfully applied to Marshy 
Harmon et al. [9] to image the Martian surface at 39 km 
resolution with the Arecibo Observatory at 12.6 cm 
wavelengths. The GSSR used the technique at 3.5 cm 
during the 1999 opposition [10]. Dual-polarization coded 
long pulse data were collected. The 1999 opposition 
provided the opportunity to observe Mars at high 
northern sub-radar latitudes, providing a good view of 
the north polar region. The campaign resulted in 14 days 
of data over the course of a month. The GSSR 25 kHz 
baud rate used in 1999 had the potential to produce 12 
km resolution cells in the polar regions. For signal-to- 
noise reasons, the actual cell size was 48 km The 
random-long-code technique is particularly useful for 
obtaining radar information away from the sub-Earth 
radartrack at moderate to high latitude. The north-south 
ambiguity inherent in delay-Doppler is still manifest, and 
so the technique really lends itself to observations of 
bright features on the surface. 

Interferometric delay-Doppler radar. Observingthe 
radarecho with more than one receiver provides a means 
to remove the north-south delay-Doppler ambiguity. Ten 
of the 63 radar tracks in the recently re-assembled radar 
dataset are interferometric radar tracks. The 
interferometric information has never been processed, 
because the signal to noise is insufficient to constrain 
both the phases and the radar scattering parameters. The 
new topographic data from the Mars Orbiter Laser 
Altimeter (MOLA) on the Mars Global Surveyor (MGS) 
spacecraft offer the best means to analyze these unused 
data to make radar maps that extend the radar properties 
coverage some 3 to 4 degrees beyond the sub-earth radar 
track. This will be a significant expansion ofthe dataset, 
and is all the more warranted as the radar spatial resolu- 
tion improves away from the sub-Earth track (smaller 
range ring-Doppler strip intersections away from the sub- 
Earth track). At the outer edges the radar resolution cell 
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GSSR Mars Data: A. F. Haldemann et al. 


is of the same order ofsize as the landing siteellipses for 
future mission (approximately 20 km diameter). Progress 
on this new technique is reported by Larsen etal. at this 
meeting, in particular as regards data acquired during the 
2001 Mars opposition. 
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23 Jan 1993* 

4.9 

186 

238 

28 Jan 1993* 

4.4 

137 

179 

31 Jan 1993* 

4.0 

99 

172 

3 Feb 1993* 

3.9 

102 

168 

7 Feb 1993* 

3.8 

109 

187 

14 Feb 1993* 

3.7 

8 

113 


* Interferometric data recorded 


Table 1 . Existing GSSR Mars Radar Data for MER 


Date 

Lat. 

(deg.) 

West Longitude (deg.) 
Rise Set 

14 Sept 1990 

-6.0 

242 

289 

22 Sept 1990 

-4.8 

160 

214 

29 Sept 1990 

-4.0 

108 

131 

2 Oct 1990 

-3.8 

52 

109 

12 Oct 1990 

-3.3 

41 

89 

25 Oct 1990 

-3.6 

227 

242 

27 Oct 1990 

-3.7 

234 

286 

20 Nov 1990 

-7.6 

65 

82 

15 Dec 1990 

-12.0 

53 

92 

17 Dec 1990 

-12.2 

23 

77 

24 Dec 1990 

-12.8 

319 

20 

28 Dec 1990 

-13.0 

290 

6 

30 Dec 1990 

-13.1 

263 

313 

15 Dec 1992* 

11.0 

214 

238 

22 Dec 1992* 

10.0 

192 

218 

31 Dec 1992* 

8.6 

158 

207 

2 Jan 1993* 

8.2 

150 

188 

3 Jan 1993* 

8.0 

90 

175 

5 Jan 1993* 

7.7 

19 

68 

10 Jan 1993* 

6.8 

15 

109 

14 Jan 1993 

6.0 

308 

353 

21 Jan 1993 

5.1 

239 

273 
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Introduction: The thermal region ot the spectrum is 
one of special interest in planetary science as it is the 
only region where planetary enussion is significant. 
Studying how planetary surfaces emit in the thermal 
infrared can tell us about their physical makeup and 
chemical composition, as well as their temperature 
profile with depth. This abstract will discuss a model of 
thermal energy transfer in planetary regoliths on airless 
bodies which includes both conductive and radiative 
processes while including the time dependence of the 
solar input function. The applications of this model to 
the Moon, Mercury, and Io, as well as the rings of Sat- 
urn, will be discussed. 

The classic model [1] of thermal energy transport in 
planetary regoliths assumed all heat transfer within a 
regolith occurred by conduction only, and all radiation 
occurred at the surface only. However, for bodies 
without atmospheres, radiation can occur between 
grains of a regolith and can escape into space from the 
immediate subsurface. These processes are important 
as they can affect the brightness temperature of a body 
as well as its temperature profile with depth. Models 


which account for radiative processes, but do not in- 
clude time dependence were developed first [2,3], and 
showed that gradients could develop in the topmost 
layers of a regolith which could affect the position of 
the Christiansen feature, and which also had implica- 
tions for the actual temperature of the planetary sur- 
face; namely it could possibly be colder than previ- 
ously thought. This model extended that work by ac- 
counting for the time dependence of the solar source 
function. 

The assumptions of our model are as follows: we as- 
sume a particulate regoltih with grains close packed 
and much larger than the wavelength. Particles are as- 
sumed to scatter isotropically, and have no wavelength 
dependent properties (i.e. we assume the regolith is 
spectrally "gray”, with no strong absorption features). 
We also assume that particle properties do not vary 
within a defined layer (i.e. properties do not vary con- 
tinuously with depth). We also treat the visible and 
thermal parts of the spectrum separately. This model 
has several free parameters, most of which (single scat- 
tering albedo, for example) are fixed by the object. 
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A NEW METHOD FOR ESTIMATING THE SINGLE SCATTERING PHASE FUNCTIONS OF REGOLITH GRAINS. 
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Introduction: Hapke’s [ 1,2,3] photometric model has been 
widely used in solar system remote sensing applications for 
nearly two decades. Recently, Hapke [4] extended his model to 
describe the coherent-backscatter opposition effect and 
multiple-scattering by particles with anisotropic single particle 
phase functions (SPPF’s). A practical difficulty for retrieving 
Hapke’s model parameters from typical planet, satellite, and 
asteroid photometry data sets is that the model employs a large 
number of adjustable parameters (at least eight) that can be 
reliably constrained only for a small nunilier of planetary data 
sets in which both disk-resolved and whole-disk observations 
are available from opposition to very large phase angles [5]. 
The present work aims to reduce the number of adjustable 
parameters and preserve (or even enhance) the model’s 
accuracy and usefulness by expressing Hapke’s parameters in 
terms of more fundamental physical properties on which they 
mutually depend. The most difficult part - if this task, described 
here, is to develop a simple method for computing the effective 
SPPF for structurally complex regolith grains from optical 
constants, grain-size distribution, and avei age regolith porosity. 

The development of light-scattering models for 
irregularly shaped particles is a large, complex subject and 
many sophisticated methods, such as Discrete Dipole Approxi- 
mation (DDA) and Monte-Carlo simulations, have been 
explored elsewhere [6]. Many of these methods are 
computationally intensive and probably impractical for routine 
substitution in Hapke’s model. Here, 1 report progress in 
developing a practical, semi-empirical method for estimating 
the directional scattering behavior (i.e SPPF) of irregular 
regolith grains. The method employs Optical Transfer Function 
(OTF) techniques to model how the structural complexity of 
regolith particles broaden and attenuate the angular distribution 
of scattered light relative to that expected from ideal spherical 
particles of equivalent size and composition 

Methodology: In solar system remote sensing, OTF 
filtering is often used as an image processing technique for 
restoring spatial details in photographs obtained by imaging 
cameras. It relies on measurements made in the engineering 
laboratory and m-flight of how the actual light transmission 
and scattering behavior of a camera’s optical system deviates 
from the predicted behavior for a corresponding hypothetical 
ideal camera (one that lacks optical impel fections, internal light 
scattering, diffraction scattering, etc.) For spatial restoration of 
images, the OTF is a convolution filter that degrades the 
theoretical brightness distribution of a scene that would be 
imaged by an ideal camera to that which is observed by a real 
camera. In image-processing practice, a picture obtained from 
an imaging camera is first Fourier transformed, the measured 
OTF for the camera optical system is then deconvolved, and the 


result is inverse-Fourier transformed to produce a picture in 
which spatial detail and photometric fidelity has been restored. 

The aim of extending OTF methods to the problem 
of scattering of light by small, irregularly shaped particles is to 
provide a simple technique to adjust single particle phase 
functions (SPPF) computed from Mie theory to realistically 
approximate the SPPF of an average, naturally occurring 
regolith grain. A theoretical, perfect spherical grain can be 
considered as an analog to an ideal camera lens. The angular 
distribution of scattered light (SPPF) from the ideal spherical 
particle can be exactly computed from Mie theory given grain 
size, optica] constants, and the wavelength of light [7]. An 
average regolith grain having the same size and optical 
constants as the Mie panicle is then considered as an analog to 
a real camera lens. The regolith gram’s shape irregularities, 
internal structure, and packing with other grains are akin to 
optical defects and diffraction limit of a real camera optical 
system. A convolution function analogous to the OTF can be 
found that transforms the Mie scattering phase function to 
more realistically approximate the corresponding SPPF of the 
irregular regolith grain. Because SPPF’s are most often 
expanded in terms of Legendre polynomials (for example, to 
treat multiple-scattering of light among particles [8,4]), it is 
advantageous to substitute Legendre transforms in place of the 
Fourier transforms used in optical engineering. 

To derive OTF’s for specific examples of rough 
panicles observed in the laboratory, their measured SPPF’s and 
corresponding SPPF’ s predicted from Mie theory are expanded 
in Legendre series and a transform filter is found that maps the 
distribution of Legendre coefficients from the Mie SPPF to 
those of the rough panicle SPPF. Following the analogous 
approach used in laboratory derivation of OTF’s for optical 
systems, we treat separately each major component of 
scattering degradation (surface scattering, internal scattering 
and absorption, and diffraction) by the irregular particle. To 
first-order, the effects of surface scattering can be isolated 
from internal scattering and absorption by analyzing linear 
polarization phase curves for a given irregular particle type and 
resolving the SPPF contributions from rays scattered parallel 
and perpendicular, respectively, to the scattering plane. 
Separate OTF’s are then derived for each component of 
polarization. In cases where it was not possible to accurately 
measure the diffraction lobe from the rough particle, Allen’s 
[9] approximation for diffraction from rough particles is used 
to model it. Mie theory describes diffraction for small spherical 
particles that are isolated and independent. However, the 
amplitude of the diffraction component in an average regolith 
grain’s SPPF strongly depends on the packing density of 
particles. Hapke [10] provides the following expression for 
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REGOLITH PARTICLE PHASE FUNCTIONS: P. Heifenstein 


computing the diffraction scattering efficiency, Q D , for 
particles in proximity to one another 

i 11 

Q d - 1 - 1.209(G (p) 3 where G = 4 a 2 /3 n~ v 

v and o are, respectively the volume and geometric cross- 
section of each grain, and 0is the packing factor. The ratio of 
Q d values computed for the irregular grain and Mie particle, 
respectively, is used to adjust the amplitude of the diffraction 
contribution from the rough grain prior to extracting the 
diffraction OTF. 

Preliminary Results and Conclusions: Preliminary tests 
of the OTF method with published intensity and polarization 
phase curves for rough, irregular particles [11,12,13] are 
encouraging. The results suggest that, in many cases, 
satisfactory OTF’ s for use with Hapke’s model can be obtained 
by simply re-scaling the power spectrum of Mie SPPF’s to 
lower Legendre orders and reducing the amplitude of the re- 
scaled low-order Mie Legendre coefficients to accommodate 
energy losses due to absorption. This is consistent with the 
intuitive notion that grain surface irregularities and internal 
particle structure distribute scattered light over broad ranges of 
scattering angle. 

For the OTF method to be realized in Hapke’s model, it 
will be necessary to conduct a systematic investigation of the 
Legendre-spectrum scaling rules and their possible dependence 
on grain-size, optical constants, and particle packing. 
Ultimately, it is hoped that simple empirical relationships will 
be found for the scaling rules so that Hapke’s current totally- 
empirical SPPF (i.e. a double Henyey-Greenstein function 
[ 1 1,4]) can be replaced by Mie phase functions that have been 
OTF-filtered to realistically model SPPF’s of irregular regolith 
grains, A special advantage of this approach over the use of the 
empirical Henyey-Greenstein function is that it provides 
physically-based a priori information about how a realistic 
regolith grain scatters light at phase angles that may not be 
available in planetary photometry data sets from which model 
parameters are being retrieved. Physically-based modeling of 
SPPF’s in Hapke’s equation would greatly reduce the risk of 
obtaining meaningless solutions for parameters that are 
retrieved from data sets with limited phase coverage. A wider 
suite of simultaneous laboratory measurements of polarization 
and intensity phase curves for isolated, irregular grains of 
different sizes, compositions, and structural complexity are 
badly needed to complete this study. Such measurements 
would be especially valuable if made from the fine particle 
fraction of returned lunar soil samples [14]. 
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Introduction: We present resubs from an ongoing 
program to perform UV measurements (215.0 and 
237.0 nm) of the Moon at varying solar phase angles to 
understand the lunar phase curve al ultraviolet wave- 
lengths. We use new observations Irom the Ultraviolet 
Spectrometer (UVS) aboard the Student Nitric Oxide 
Explorer (SNOE) combined with existing observations 
from the Galileo UVS. The lunar UV phase curve can 
be used to further understand the scattering properties 
of the lunar surface. The Moon’s scattering properties 
at visible wavelengths are well understood [1]; study- 
ing scattering properties at shortei wavelengths may 
provide insight into the roles of volume scattering vs. 
surface scattering and how weathering processes may 
affect scattering properties [2]. The UV lunar phase 
curve can also be helpful for UV observers, as the 
Moon is often used as a UV calibration source, but the 
UV brightness variation with phase angle has not been 
well understood. 

SNOE UVS Moon Observations: The SNOE sat- 
ellite has been in Earth orbit since March 1998. It was 
built and is currently operated at the Laboratory for 
Atmospheric and Space Physics at the University of 
Colorado. The primary objective o?' SNOE is to meas- 
ure nitric oxide in the Earth’s atmosphere, as well as to 
measure solar soft x-rays [3]. The SNOE spacecraft is 
in a low altitude, polar Earth orbit. The UVS on the 
spinning spacecraft is in a configuration such that it 
serendipitously observes the Earth's moon once per 
month. The angle between the orbital plane and the 
sun varies with time; this beta angle was initially -25°. 
Thus the lunar measurements geneially occur 1-2 days 
before full moon. The offset angle between the orbital 
plane and the UVS viewing direction can be adjusted, 
thereby varying the solar phase mgle at which the 
moon is observed. The UVS slit is 0.75 o x0.071°, so as 
the slit scans over the moon, 7-9 samples are obtained. 
The UVS has two spectral channels, 215.0 nm and 
236.5 nm, each with a bandpass of 3.6 nm. The UVS 
has made measurements of stellar sources, and cali- 
bration of both channels is well understood. 

The SNOE lunar data set conststs of both seren- 
dipitous and dedicated observations. The initial moon 
observations performed by SNOE were purely seren- 
dipitous. During the first month of the mission, the 
moon was observed, ~2 days before full moon, as 
SNOE passed over the north pole. 

The first dedicated Moon observation was per- 
formed in February of 2000. The observational timing 
was changed during the nighttime portion of the orbit 


so that, rather than scanning the limb, the UVS 
scanned the portion of the sky (in the zenith direction) 
where the Moon was projected to be, based on the 
knowledge of the Moon’s right ascension and declina- 
tion. Other dedicated lunar observations have been 
performed, to observe the Moon at a smaller phase 
angle than the orbital beta angle. 

The solar phase angle during each Moon observa- 
tion is determined using data from the Solar Angle 
Sensor on the Solar X-Ray Photometer, which meas- 
ures the angle between the sun and the viewing direc- 
tion of the instruments. 

Galileo UVS Moon Observations: To increase the 
phase angle range of the Moon observations, we have 
combined the SNOE data with data from the Galileo 
UVS encounters with the Moon. The phase angle 
range from the Galileo flybys is ~19°-70°, providing 
some overlap with the SNOE data set in addition to 
valuable larger phase angle coverage. The Galileo 
UVS F-channel covers the 200-320 nm wavelength 
range; in this analysis we use only the data from the 
ranges covered by the SNOE UVS channels. 

Analysis: The phase angle range of the current data 
set is 3.5-70°. We model the data using the Hapke 
model (see the figure below) to derive a determination 
of the ultraviolet photometric parameters of this rela- 
tively well-understood surface. These Hapke parame- 
ters can then be used for comparison with visible re- 
sults to understand changes in the scattering properties 
of the surface with decreasing wavelength. 
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Introduction: In this study, we use ultraviolet ob- 
servations from the International Utlraviolet Explorer 
(IUE) and the Galileo Ultraviolet Spectrometer (UVS) 
to compose the ultraviolet solar phase curves of the icy 
galilean satellites. Broadband rotation phase curves 
from 0,26 to 0.32 microns are constructed in order to 
examine the rotational behavior of the icy Galilean 
satellites in the ultraviolet. After normalizing the rota- 
tional variations, modeling of the solar phase varia- 
tions are compared to comparable studies in the visi- 
ble. 

Background: Studies of the rotational phase curve 
variations at visible wavelengths for the icy Galilean 
satellites indicate distinctive albedo changes from 
leading to trailing hemisphere. Processes such as the 
preferential bombardment of the leading hemisphere 
by micrometeorites and/or the preferential bombard- 
ment of the trailing hemisphere by co-rotating ions 
within Jupiter’s magnetosphere have been proposed as 
mechanisms for creating these albedo variations. Ul- 
traviolet observations by IUE [1] and the Hubble 
Space Telescope (HST) have measured absorption 
features attributed to surface chemistry resulting from 
ion interactions with the surface. 

Discovery observations with IUE in the late 1970’s 
showed an absorption feature at 0. 18 microns on Eu- 
ropa whose strength correlated witli the albedo varia- 
tions observed on this satellite, and whose source was 
proposed to be the implantation of magnetospheric 
sulfur ions into the water-ice surface [2]. Later meas- 
urements with HST [3] were compared to laboratory 
spectra of irradiated water-ice and SO 2 /H 2 O layered 
frosts [4], to show that the observed spectral feature is 
better matched by the layered frosts, implying an en- 
dogenic origin. 

Ganymede’s rotational albedo variations are similar 
to Europa’s, with the leading hemisphere being gener- 
ally brighter than the trailing hemisphere. Ganymede’s 
trailing hemisphere, however, has a different absoption 
feature than Europa’s. Ganymede’s absorption feature 
has been found to be ozone, and is most concentrated 
in the trailing hemisphere polar regions and at large 
solar zenith angles., commensurate with production by 
ion bombardment and destruction by solar photon in- 
teractions [5, 6]. HST visible observations of Gany- 
mede have also detected oxygen on the trailing hemi- 
sphere at low latitudes, consistent with production via 
ion bombardment processes [7]. 


Callisto’s rotational albedo varations are opposite 
to those seen on Europa in the visible, where the lead- 
ing hemisphere is brighter than the trailing hemisphere. 
This is also the case in the ultraviolet. IUE and HST 
observations in the ultraviolet have also detected a 0.28 
micron absorption featuer, linked to S0 2 [8, 9]. This 
absorption feature shows no correlation to the rota- 
tional albedo variations, however the IUE observations 
show some spatial correlation to possible neutral wind 
interactions suggestive of neutral sulfur implantation 
into Callisto’s water-ice surface. The strongest 0.28 
micron absorption seen in the IUE data was in the re- 
gion containing Valhalla basin. Galileo UVS also de- 
tected this absorption feature on the leading hemi- 
sphere, but found no definitive correlation with surface 
feature [10]. 

Data Set and Analysis: Rotational phase curves 
were derived for each satellite for each wavelength bin 
(0.26, 0.27, 0.28, 0.29, 0.30, 0.31 and 0.32 microns), 
using IUE data in the 10° phase angle range. These 
rotational phase curves were then used to correct the 
entire data set (including additional solar phase angles) 
for each satellite to particular longitudes for photomet- 
ric analysis. The Europa data set was the most exten- 
sive and allowed us to analyze all four longitude re- 
gions: the leading and trailing hemispheres, as well as 
the Jovian and anti-Jovian hemispheres. The Gany- 
mede and Callisto data sets were more limited; for 
these satellites, we used the rotational phase curves to 
correct to 90° W longitude, and derived a global set of 
photometric parameters. This global set of photometric 
parameters was also derived for Europa, for compari- 
son with Ganymede and Callisto. 

Results: In this presentation, we discuss the rota- 
tional and solar phase curves we have derived for the 
three icy galilean satellites. We have fit Hapke pho- 
tometric models to the solar phase curves for Europa, 
Ganymede and Callisto, and to four longitude regions 
on Europa. We intercompare these to look for varia- 
tions related to differences in exogenic processes as 
discussed above. Although the data sets are not ideal - 
the data are sometimes noisy and phase angle coverage 
tends to be sparse, this work represents the first in- 
depth ultraviolet analysis of the photometric parame- 
ters of the icy galilean satellites, and the results are 
useful for understanding variations in scattering phe- 
nomena with wavelength. 
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Introduction: The spin-vector of Venus H was 
investigated by various Earth- and spacecraft-based 
techniques [1-3] (see also refs, in [3]) but until now 
no experimental data have been obtained on varia- 
tions of both in magnitude and direction because of 
insufficient accuracy and too long measuring interval 
(much longer than the period of variations). In this 
work a new 1 on principle ground-based radar interfer- 
ometric technique named Radar Speckle Displace- 
ment Interferometry (RSDI) is projiosed to measure 
precisely instantaneous vector compments of H and 
their variations with time. 

The essentials and theory of RSDI were devel- 
oped in [4-6] (see also refs, in [6]). The technique is 
based on a so called far coherence (speckle displace- 
ment) effect for speckled radar fields scattered from 
randomly rough surfaces of moving objects and aims 
at precise measurement of their instantaneous rota- 
tional-progTessive motion parameters. 

Far Coherence or Speckle Displacement effect: 
Let us illuminate Venus from Goldslone (USA) by a 
monochromatic radiowave at wavelength X = 13 cm 
(Goldstones 3.5 cm is absorbed by the Venus’ at- 
mosphere). The randomly rough surface of Venus can 
be considered as a reflecting antenna with a random 
speckled backscattering pattern. The average size of a 
speckle near Earth will be 1 = XR/a. ~ 3.5 km, where 
R - 40 million km - the Earth-Venus distance at infe- 
rior conjunctions, a t - 1500 km - the effective size of 
subradar reflecting area in the direction orthogonal to 
the line-of-sight. The Fraunhofer pattern does not 
depend on distance while the Fresnel pattern changes 
with the distance to Venus. When Venus rotates with 
respect to Earth all its Fresnel speckle pattern range 
layers rotate as well remaining “free zed” for a rather 
long time. One of these “freezed” speckle layers cor- 
responding to the current Earth Venus distance 
sweeps across the Earth surface. It was shown in [4] 
that the speckle pattern displacement in a ‘Treezed” 
state can be characterized by a length L = 2R(A/a T )° 5 
- 130000 km, where 2 ^ ~ 50 km - the effective size of 
subradar reflecting area along the line-of-sight. It 
follows from L » 12800 km (the diameter of Earth) 
that the speckle pattern remains “freezed” while its 
displacement all over the surface of Earth. In other 
words high correlation can be observed between radar 
echoes received at opposite sides of Earth separated 
by thousands kms and this effect can he referred to as 
far coherence of scattered speckled iadar fields [7-9 J. 


RSDI: The main principle of RSDI as applied to 
a planet’s rotation measurement is as follows. If two 
antennas are placed along with the speckle displace- 
ment direction they can see the same speckles passing 
through them and cross correlation of the received 
echoes will give the velocity v and direction v/v 
(the receiving baseline orientation changes with time 
slowly due to rotation of Earth) of speckle displace- 
ment related to the transversal (orthogonal to the line- 
of-sight) vector component of the planet’s instan- 
taneous spin- vector £1 [4-151. It was shown in [5] that 
the above RSDI-principle coincides with the proce- 
dure of optimum processing of scattered random 
speckled radar fields under monochromatic illumina- 
tion and the minimum of measuring time to estimate 
instantaneous rotational-progressive motion parame- 
ters and their fluctuations with time. This means that 
RSDI may reach the limit accuracy of planet’s instan- 
taneous spin -components estimation when the meas- 
uring time is no more than a few minutes and there 
exists no another Earth-based radar interferometer 
which could have accuracy better than RSDI in these 
conditions. 

RSDI accuracy: Let us calculate the expected 
one-look accuracy for the radar interferometers Gold- 
stone-Green Bank (USA) and Evpatoria (Crimea, 
Ukraine)-Effelsberg (Bonn, Germany)-Bear Lakes 
(Moscow, Russia) in Europe (X = 6 cm, values are 
given below in brackets) in accordance with the 
methodology [5-10, 12, 13]. The potential limit of Sit 
estimation was derived in [5] where it was shown 
using the Maximum Likehood criterion that the limit 
one-look relative accuracy o of ft t magnitude (index 
1) and direction (index 2) estimation by a two- 
element radar interferometer can be expressed ana- 
lytically as 

lv.o v 

a 1>2 = — 10' 6 (~ 1.5 TO' 6 ) (1) 

qb v n 

where the baselength b ~ 3000 km (~ 2000 km), l v is 
used for Oi and 1 0 for o 2 , l v ~ 1 0 ~ 1 due to nearly 
spherical shape of Venus, v n is a part of v caused 
by rotation only (under V = 0), q = Q° 5 , Q - power 
signal-to-additive noise ratio (snr) at the correlator 
output, defined as 

Q = N/(l/Q, + 1/Q 2 + I/Q 1 Q 2 ) - 7 • 10 5 (~ 1.7 • 10 5 ) 
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N = FT ~ 450 (~ 450), T ~ 30 s (- 14 s) - correlator 
integration time, Q u - power snr at the correlator 
input within Doppler bandwidth F ~ 15 Hz (~ 32 Hz) 
for Green Bank or Effelsberg (index 1) and Goldstone 
or Evpatoria/Bear Lakes (index 2), Qi = P, / P n ~ 
4600 (~ 1060), Q 2 = P 2 / P n - 2300 (~ 535), noise 
power P„ = k B T s F ~ 5.4 • 10' 2 ' W (~ 2.3- 10‘ 20 W), 
T s - 25 K (~ 50 K) - system temperature, k B - 
Boltzmann’s constant, from the radar equation [3] 

P,- G 12 - G 2 - A. 2 * S 

Pij= 

(4 nf • R 4 

Pi - 2.5 • 10' 1 ' W - echo power at Green Bank, P 2 ~ 
1.2iO i7 W - echo power at Goldstone, P t - 450 kW 
(- 100 KW) - transmitted power, Gj ~ 1.8- 10 6 (~ 
8.4 ■ 10 6 )Green Bank (Effelsberg) antenna gain, G 2 ~ 
9.1- 10 5 (~ 4.3- 10 6 ) - Goldstone (Evpatoria/Bear 
Lakes) antenna gain, S - 10 7 km 2 - Venus’ radar 
cross section. Eq. (1) is consistent with the well 
known results for time delay estimation of random 
signals [5]. 

The optimum processing algorithm leading to the 
limit accuracy Eq. (1) is too complicated and a qua- 
sioptimum one can be used in practice. Computer 
simulation showed that loss in accuracy while this 
transition is 25... 30% [6] and instead of the numerical 
values in Eq. (1) we will have 

o~ 1.3 ■ 10" 6 (-2 1CT 6 ) (2) 

Values (2) mean that one-look (for - 1 min obser- 
vation) accuracy in determination of Venus’ polar 
axis orientation as projected onto a transversal plane 
orthogonal to the instantaneous line-of-sight with one 
of the above two-element interferometers will be 
300. ..400 mas (milliarcseconds). If one interferometer 
is used n times (=days) within a single conjunction 
the accuracy can be - n° 5 times better. With m com- 
parable interferometers accuracy will be another m 05 
times higher. For n = 16 and m = 4 the resulting accu- 
racy can be as fine as 40... 50 mas (!). It is not easy to 
achieve these limits even through space missions to 
the planet. In addition variations of Venus’ polar axis 
(if present) can be studied. The resulting relative ac- 
curacy in rotational velocity of Venus can approach ~ 
2 * 10" 7 again with the possibility to study variations. 
RSDI-observations during several successive inferior 
conjunctions may lead to a few tens of mas in polar 
axis orientation and ~10 7 in rotational velocity. More 
powerful transmitters should give further improve- 
ments. 


Experiments: The effect of far coherence 

(speckle displacement) discovered theoretically in 
[4,5] is now proven in fact both in Europe [15] and in 
the USA [16,17]. The RSDI data obtained on Mer- 
cury in accordance with [4,5,10,11] in the USA in 
May- June 2002 [17] confirm the theory developed in 
[4-15] in main part and are now processed to deter- 
mine Mercury’s obliquity and librations [18]. RSDI 
accuracy is expected to be several orders of magni- 
tude better than that obtainable through other Earth- 
based radar techniques [5-12,14,15] and can be com- 
parable [19] with the accuracy of the two upcoming 
space missions to Mercury: MESSENGER (USA) 
and BepiColombo (Europe-Japan). 

RSDI experiments on Venus based on the above 
theory are possible both in the USA and in Europe as 
early as this 2002 fall. Further opportunities for 
maximum precision measurements will take place in 
spring/summer 2004 while the next inferior conjunc- 
tion of Venus with Earth. Future RSDI experiments 
can be aimed at substantial improvement in knowing 
the Venus’ polar axis orientation in the sky and pre- 
cise estimation of the spin-vector variations with time 
both in magnitude and orientation. 
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Introduction: Planetary rotations and related pa- 
rameters have been measured throughout years by 
radar and optical imaging techniques remotely from 
Earth or spacecraft. It is well understandable from our 
everyday life that spin components of rotating bodies 
can be estimated by obtaining their images at differ- 
ent moments of time and then comparing them to 
each other. Our eyes employ this principle for thou- 
sands years with Sun as a source of illumination. Note 
that Sun is an incoherent source while laser or radar 
are coherent. At what extent the imaging approach is 
good in these “coherent” conditions and if there exists 
a better approach? An answer was given in [1,2]. 

Theory: A new approach does not necessarily 

lead to a substantial improvement in accuracy but to 
improve accuracy by several orders of magnitude as a 
rule one must use a new physical principle. As ap- 
plied to planetary rotation measurements by Earth- 
based radar a new principle can lx: based on a so 
called “freezed” speckle pattern displacement effect 
or simply speckle displacement effect (SD-effect) [1]. 
Let us illuminate Mercury by a monochromatic ra- 
diowave at wavelength X = 3.5 cm Due to coherent 
illumination and rough surface of Mercury the back- 
scattering pattern is speckled like those well known 
and easy observable in laser optics experiments [3]. 
The mean effective transversal (orthogonal to the 
line-of-sight) size of a speckle near Earth is l t ~ kR/a t 
- 3 km, where the Earth-Mercury approach R - 90 
million km, a t - 1000 km - the effective transversal 
diameter of subradar backscattering area on Mercury. 
The mean radial (along the line-of-sight) size of a 
speckle can be estimated as l r - 1,71 ~ 260,000 km 
characterizing the thickness of transversal speckled 
layers inside which the Fresnel backscattering pattern 
does not change with range. 

What happens with the speckle pattern when the 
object rotates was considered in [ 4 | for a particular 
case of cylinder rotation with the spin-axis orthogonal 
to the line-of-sight. For arbitrary rotational- 
progressive (spin-orbital) motion and arbitrary shape 
of an object the problem was solved in [1] where it 
was shown that in the most general case excluding 
some theoretical geometries and motions nearly im- 
possible in practice of Earth-based radar astronomy 
the reflected speckle pattern first sweeps across the 
receiving aperture plane in a nearly “freezed” state 
for a rather long distance L ~ 2R(Xyar) 05 ~ 150,000 
km, where for Mercury a r ~ 50 km - effective radial 


size of the subradar backscattering area, and then 
decorrelates (decorrelation starts immediately but 
becomes visible after the speckle pattern overcomes 
the distance L; also there are particular geometries 
free of decorrelation, e.g. when the transmitter is 
placed at the center of the receiving baseline). L is 
much longer than the diameter of Earth so the speckle 
pattern while its displacement in any case must re- 
main “freezed” all over the Earth’s surface. If this 
“freezed” SD-effect really exists in the global Earth 
scale than the correlation of radar echoes obtained at 
the two antennas, which form a two-element interfer- 
ometer with the baseline that due to rotation of Earth 
can be instantaneously colinear with the SD-direction, 
can provide two observables, the time delay and the 
epoch assotiated with maximum correlation, both of 
which may constrain Mercury’s spin state to accuracy 
as high as one part in 100,000 (a few arcseconds) [5- 
13]. 

Experiments: To verify the SD-effect in practice 
of radar astronomy a number of radar interferometric 
experiments were proposed in Europe [5,6,12], Euro- 
Asia [2,5,6] and in the USA [5-13]. Alternative radar 
techniques, e.g. repeat orbit radar interferometry 
[14,15] were proposed and discussed as well (for 
comparison of the techniques see [16]). During a year 
long discussion of the above SD-approach with S. J. 
Peale he organized a group of radar astronomers who 
successfully carried out the Goldstone - Green Bank 
SD-experiment on Mercury in May- June 2002 [17]. 
The theory and potential limits of SD-experiments 
including Goldstone - Green Bank can be found in 
[8,9,12,13]. A large cross-correlation between Green 
Bank and delayed Goldstone radar echoes was clearly 
detected in accordance with the above theory at least 
on two days: on May 22 and June 02. One of these 
cross-correlation peaks caused by the SD-effect is 
shown in [17]. It follows from the existence of the 
peak that the effect of “freezed” speckle pattern dis- 
placement discovered theoretically in [1,2] is now 
proven in fact over the base length ~ 3000 km and 
further Goldstone - Green Bank data processing (now 
in progress) may lead to the first in history (!) deter- 
mination of Mercury’s obliquity and librations as 
discussed in [7-9,12,13]. 

Discussion: The intermediate experimental result 
[17] allow estimation of the expected accuracy which 
can be compared with the potential limits [2,7- 
9,12,13]. In fig. 1 a correlation graph is shown sche- 
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maticly with necessary parameters estimated from the 
experimental graph [171 (a periodic disturbance of the 
amplitude ~ 0.07 with the period - 15... 20 times lar- 
ger than the width of the correlation peak seen in the 
graph [17] was removed). 



Fig. 1. Goldstone - Green Bank cross- 
correlation caused by the SD-effect 


As it follows from the theory [ 1-2,5-13] the 
maximum of the Goldstone - Green Bank cross- 
correlation peak should take place at the time delay t 
~ 10... 13 s needed by the speckles to come from 
Green Bank to Goldstone and its width can be ex- 
pected within Ax ~ 10... 15 ms (both parameters were 
not given in [17]). The maximum value of cross cor- 
relation is - 0.73 and the noise amplitude is - 0.04 
(fig. 1) so the signal-to-noise ratio is q - 18 and this 
may lead to the relative accuracy [2,5-13] in instanta- 
neous rotational velocity Oi ~ Ax/qx ~ 6 10‘ 5 and in 
obliquity 0 a 2 ~ 12 arcsec. The obliquity itself can be 
expected to be as small as a few arcminutes [16] so 
the relative accuracy in 0 will be ~ 10 % (for 0-2 
arcmin). The physical libration amplitude (p is ex- 
pected within 20... 60 arcsec [16] so the one look rela- 
tive accuracy in (p in accordance with the methodol- 
ogy [9,13] will be - 50 % (for cp - 30 arcsec). So al- 
ready the intermediate result [17] is outstanding al- 
lowing the first in history measurement of 0 with 
good accuracy. At the same time potential limits for 
the May-June 2002 Goldstone - Green Bank condi- 
tions are ~ 10 5 and a 2 ~ 2 arcsec [9,12,13] and we 
all do hope that further quasioptimum processing of 
the echoes will reach these limits allowing determina- 
tion of both 0 and cp (in case of damped free libration 
[16]) with fractional uncertainties - 3 % and - 15 % 
respectively. 
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Introduction: Emission lines of atmospheric con- 
stituents of the Moon and Mercury (especially sodium 
and potassium) must be absolutely calibrated before 
the abundances can be obtained. It is convenient to use 
the intensity of the surface at an adjacent wavelength as 
the standard, and to use Hapke theory to obtain this 
intensity. This procedure and some of its pitfalls are 
discussed here. 

Methods Used: Observations of these emission 

lines in the presence of a strong background from the 
surface must be made at high spectral dispersion. In 
spectra of Mercury both lines and continuum are pre- 
sent and the intensity of the light from the surface can 
be obtained simply by use of a wavelength nearly ad- 
jacent to the emission line. The lunar atmosphere is 
usually observed above the bright limb. Although 
there is also a strong component of light from the sur- 
face, scattered into the instrument by the Earth’s at- 
mosphere, it cannot be used for an intensity calibration. 
Instead a second exposure is taken, with the telescope 
pointed directly at the surface. The intensity from this 
region may then be computed if the Hapke parameters 
have been obtained and tabulated. I or Mercury we use 
average parameters for the whole surface, taken from 
the review chapter by Veverka et a). [1]. We compute 
an array of rough reflectances rp and then an array of 
intensities for the entire illuminated area of the planet. 
These can then be used to compute a synthetic image 
of the planet, smeared by various assumed amounts of 
atmospheric “seeing”. Cuts are then taken across these 
images at the location of the spectrograph slit, and the 
one that fits best selected. The procedure is described 
in detail by Sprague et al. [2], who find that the 2- 
dimensional Gaussian used to represent the seeing has, 
in one example, a half-with parameter of 1.5 arcsecond. 
This large value is the result of observing the planet 
close to the horizon. Another group [3] uses similar 
procedures, though they have not been described in as 
much detail. 


Discussion: A major advantage of the method de- 
scribed is that the calibration is obtained (for Mercury) 
under the same geometrical and sky conditions as the 
measurement; for the Moon they are very nearly the 
same, differing only by the time taken to move the tele- 
scope. The accuracy does depend on the photometry 
from which the Hapke parameters were obtained, and 
on the accuracy with which these parameters represent 
the data. 

Intensities computed with Hapke theory tend to be 
brighter near the backscatter direction than those ob- 
served. If the instrument is calibrated under this condi- 
tion, it will appear to be less sensitive than it really is. 
The sodium data set of Potter and Morgan [4], which 
was taken over a wide range of lunar phase angle, 
showed a marked minimum of intensity near full phase. 
The suggestion of these authors was that the location of 
the Moon in the geomagnetic tail shielded its surface 
from sputtering by solar- wind ions. But a later meas- 
urement [51 by a different method, calibrated against a 
source inside the instrument, did not find any diminu- 
tion of intensity compared with what was found at 
smaller phase angles. It seems possible that the result 
of [4] stemmed from a computed surface intensity that 
was too bright near the full phase. A useful remedy 
would be to check the derived calibration factors of the 
instrument, which should normally be stable over peri- 
ods of days to months unless there are big changes in 
sky conditions.. 
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Summary: Multispectral imaging of the Mars Path- 
finder (MPF) radiometric calibration targets (RCTs) by 
the Imager for Mars Pathfinder (IMP) camera provided 
sufficient data to model the photometric changes asso- 
ciated with variable deposition of atmospheric dust on 
the RCTs. Sky models were used to correct for varia- 
tions in diffuse/direct illumination conditions during 
the mission, and two-layer Hapke radiative transfer 
models were developed to constrain variations in spec- 
tral properties associated with dust deposition on the 
RCTs. The results are consistent with dust deposition 
during the MPF mission amounting to dust optical 
depth (tau) values < 0.5 (i.e., <10 am effective thick- 
ness of 10-20 pm particles for porosity 0.80), although 
variability in the dust thickness with time may result 
from erosive/depositional phases of aeolian activity. 
Background: The IMP camera was a multispectral 
stereo imaging system covering the 440-1000 nm re- 
gion in 12 wavelengths. Two RCTs were placed on the 
MPF lander to provide a direct reference for radiomet- 
ric and color calibration of IMP images (Figure 1) [1]. 
All surfaces at the MPF landing site were subject to 
dust deposition from the atmosphere throughout the 83- 
Sol mission. Images of the RCTs were acquired regu- 
larly to assist in calibration of images to relative reflec- 
tance. Contamination of the RCTs with dust could 
therefore influence image calibration. Radiative trans- 
fer models of sky illumination that incorporate opacity 
measurements have been developed to help compen- 
sate for variations in illumination and the relative con- 
tributions of diffuse and direct skylight on both scene 
and RCT images [2-5]. Two-layer radiative transfer 
models of dust-coated surfaces have been developed to 
distinguish dust and substrate spectral characteristics 
and dust thickness [e.g., 6]. 



Figure 1. (a) IMP image of lower radiometric calibration 
target (RCT) showing shadow post, white, gray, and 
black rings, and color chips; (b) IMP image of the upper 

RCT, (spacecraft cabling blocking portion of view). 

Methodology. We initially approached this problem 
using images of the RCTs acquired early in the mission 


(Sol 3) to calibrate to relative reflectance all RCT im- 
ages taken at nearly the same local solar time (LST) 
throughout the mission. Using this approach we dem- 
onstrated that dust indeed coated the RCTs and red- 
dened its spectrum with time [7]. However, we treated 
these as preliminary results owing to the effects of scat- 
tered light from the lander and ongoing improvements 
made to the sky and two-layer models. Our current 
approach used calibrated radiance data from the RCTs 
in combination with improved radiative transfer mod- 
els. We again used the Sol 3 observations as a baseline 
to observe later variations. Knowing the photometric 
properties of the RCTs from laboratory measurements 
[1], we used the sky model to determine the expected 
radiance of the targets under Mars illumination condi- 
tions. These values were compared to data extracted 
from calibrated radiance images of the targets. Out- 
liers between modeled and observed data found to re- 
sult from observational problems (e.g., dropped image 
data, IMP cabling in the field of view; Figure 1) were 
removed from the dataset. The resulting data set was 
provided to the two-layer model. 

Two-layer model. We used a two-layer version of the 
Hapke model [8, section 9.D.3 ] to compute bidirec- 
tional reflectance for both the clean substrate (from lab 
measurements) and RCTs throughout the mission. We 
first studied cases where dust accumulation was as- 
sumed to be linear as a function of time over the dura- 
tion of the mission. Accumulation rates derived from 
these models typically ranged from 0.001 to 0.002 op- 
tical depths per sol. For dust grain sizes between 5 and 
20 pm and 80% porosity, these tau rates correspond to 
mean accumulation of 0.02 to 0.1 pm per sol. We ran 
models with isotropic single scattering phase functions 
as well as with phase functions based on atmospheric 
dust scattering properties [e.g., 2-4]. For some runs, 
we used dust absorption coefficients determined for 
palagonites from previous studies [6,9], while in others 
we left the dust single scattering albedos as free pa- 
rameters. We also ran test cases where dust thickness 
for each sol was a free parameter, rather than assuming 
a simple linear accumulation with time. 

Results. Figure 2 shows representative observed 
spectra of the white RCT ring on Sols 3, 42, and 83 
compared to model results using an isotropic phase 
function for constant dust accumulation rate. A de- 
crease in 443 nm reflectance values with time resulting 
from dust deposition is apparent. Sol 83 model values 
fit more poorly, potentially due to greater variability in 
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atmospheric opacity compared to earlier in the mission 
[2-4]. Figure 3 shows results from the same model as 
a function of time, where the RCT white ring's 443 ran 
reflectance decreases with time and the black ring’s 
671 nm reflectance increases as dust accumulates. 



Slight variations between the observations and 
linear dust accumulation rate models suggested that 
models in which dust thickness ( tau ) was allowed to 
vary on a daily basis were warranted. Figure 4 dem- 
onstrates the results from one such model, in which 20 
pm particles and an isotropic phase function were used 
to model the dust. The total dust accumulation of ~0. 1 
tau (corresponding to ~6 pm in this model) over 83 
Sols is consistent with previous estimates of 20-25 pm 
per year dust accumulation based on atmospheric fall- 
out rates after dust storms, laboratory dust-coating ex- 
periments, and observations of dust obscuration at the 


MPF site [10-12]. Daily variations in tau may be real 
or the effect of noise in the observational data, the sky 
model, and/or the two-layer model. We are currently 
studying how these different error sources may vary 
among the RCT surface types (white, gray, black) and 
how those uncertainties influence model results. 



Figure 4. Computed dust optical depth tau variations 
(and corresponding coating thickness for assumed poros- 
ity 0.80 and grain size D= 20 pm) as a function of Sol us- 
ing a model with an isotropic phase function and with a 
variable dust accumulation rate. 

Conclusions. Radiometric ally corrected IMP images 
of the Mars Pathfinder RCTs can be used to track the 
amount of deposited (or eroded) atmospheric dust at 
the MPF site. Two-layer radiative transfer models 
provide constraints on the dust opacity as it relates to 
grain size, porosity, and dust thickness, and may help 
constrain the spectral properties of the deposited dust. 

Future work will incorporate new laboratory 
photometric measurements of palagonite-coated sur- 
faces (both rocks and spare RCT materials to be used 
with MER Pancam instrument). More temporally- and 
radiometric ally-control led observations of the MER 
RCT target are planned for that mission. 
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Introduction: Space weathering is a term used to in- 
clude all of the processes which act on material ex- 
posed at the surface of a planetary or small body. In 
the case of the moon, it includes a variety of processes 
which have formed the lunar regolith, caused the 
maturation of lunar soils, and formed patina on rock 
surfaces. The processes include micrometeorite impact 
and reworking, implantation of solar wind and flare 
particles, radiation damage and chemical effects from 
solar particles and cosmic rays, interactions with the 
lunar atmosphere, and sputtering erosion and deposi- 
tion. Space weathering effects collectively result in a 
reddened continuum slope, lowered albedo, and at- 
tenuated absorption features in reflectance spectra of 
lunar soils as compared to finely comminuted rocks 
from the same Apollo sites [1,2]. However, the re- 
golith processes that cause these effects are not well 
known, nor is the petrographic setting of the products 
of these processes fully understood. 

Although it was initially proposed that impact vit- 
rification was the main darkening agent acting in the 
lunar regolith, Hapke and co-workers showed that this 
process had only a minimal effect on regolith optical 
properties [3,4]. In a very important set of experi- 
ments, Hapke et ai showed that thin coatings produced 
by vapor deposition of thermally evaporated analogs, 
and by the deposition of vapors generated by solar 
wind sputtering have optical properties that resemble 
those of the lunar fines [5]. In both of these deposi- 
tional environments, the postulated darkening agent 
was submicroscopic Fe metal grains which are strong 
absorbers of visible wavelengths [5]. Although there 
existed extensive theoretical data to support an abun- 
dance of vapor-deposited material in the lunar regolith, 
the concept was resisted because of a lack of direct 
evidence of this material in lunar soils. 

An interesting confluence of ev ents occurred in the 
early 1990’s. First, came the discovery of vapor 
deposited coatings on lunar regolith grains by Keller 
and McKay [6], who showed that amorphous coatings 
from 50-100 nm thick containing fine-grained Fe metal 
(1-10 nm in dia.) were common in the fine size fraction 
of several mature lunar soils. The other discovery was 
the recognition that the optical properties of lunar soils 
were dominated by fine grain sized material (<45 urn 
fraction) by Pieters and coworkers [7]. These discov- 
eries led to coordinated studies that looked at the min- 
eralogy, chemistry, and optical properties of lunar soils 
as function of composition, maturity, and grain size 


fraction [8,9]. One of the major revelations from these 
studies was the recognition that much of the nanophase 
Fe metal is surface-correlated especially in the finest 
size fractions, and that it was this nanophase Fe that 
dominated the optical properties of the soil. 

What we know: Nanophase Fe° has been strongly 
implicated as one of the main agents in the modifica- 
tion of optical properties of lunar soils, associated with 
space weathering Recent work has shown that it the 
size, abundance, and petrographic setting of the nano- 
phase Fe° that determines the magnitude of the optical 
effects [10,11]. Nanophase Fe° occurs in two main 
petrographic settings in fine grained lunar soils: 1) as 
inclusions within agglutinitic glass produced by mi- 
crometeorite impact, and 2) as inclusions in thin amor- 
phous rims surrounding soil grains that are formed by 
vapor deposition and irradiation effects. The average 
size of the nanophase Fe° in agglutinitic glass is a fac- 
tor of ~2 larger than that observed in the amorphous 
rims. The average size of the nanophase Fe°, regard- 
less of whether it is surface- or volume-correlated, ap- 
pears to have a major effect on the optical properties of 
the host grain. The occurrence of Fe° inclusions that 
are >10 nm in diameter results in strong darkening in 
the corresponding reflectance spectra, whereas it is the 
smaller Fe° grains <5 nm diameter that give rise to 
strongly reddened slopes in the spectra [12]. It has 
been shown that in the fine size fraction of mature 
soils, nearly every grain has nanophase Fe metal within 
100 nm of the grain surface [13]. The same situation 
exists for patinas on lunar rocks, which also contain 
abundant nanophase Fe metal on exposed surfaces 
[14]. 

All of the soil particles that we have analyzed to 
date have some accretionary material (e.g. splash glass, 
vapor deposits, sputter deposits, etc.) on their surfaces, 
although the amount of material can be highly variable. 
A common characteristic of the accretionary material is 
the occurrence of nanophase metal as randomly ori- 
ented inclusions or in layers. In the coated particles, the 
metal grains are concentrated in thin (50-150 nm thick) 
rims surrounding mineral grains (mainly plagioclase 
and augite). The large compositional difference be- 
tween rims and their hosts led [6] to propose that the 
inclusion-rich rims on mineral grains likely result from 
depositional processes in the lunar regolith, either con- 
densation of impact-generated vapors or sputter 
deposition. This accretionary material can have 
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sition. This accretionary material can have profound 
effects on the optical properties of the soil grains. 

Radiation effects can also produce surface- 
correlated nanophase Fe metal in Fe-bearing minerals 
such as ilmenite and olivine. Transmission electron 
microscope (TEM) analysis of the ilmenites shows that 
the grains are surrounded by altered rims up to 200 nm 
thick where Fe has been preferentially removed from 
the ilmenite surface, nanophase Fe metal grains have 
been produced, and the Ti has been partly reduced to 
Ti 3+ . While these altered rims contain submicroscopic 
metal grains, they do not have a large effect on the 
reflectance data. A more dramatic effect is observed in 
olivine grains. We have observed altered rims on oli- 
vine grains where much of the Fe within ~50 nm of the 
surface has been reduced to nanophase Fe metal sur- 
rounded by Mg-silicate glass. Spectra from the “al- 
tered” olivine grains are much darker than pristine oli- 
vine of similar composition. 

Recent work has demonstrated that these space 
weathering effects (amorphous rims with nanophase Fe 
metal) can survive lithification, compaction and brec- 
ciation, at least for some lunar breccias [15]. 

What we do not know: An extensive observa- 
tional database now exists on the space weathering 
effects in the lunar regolith in terms of changes in opti- 
cal properties, surface microstructure, and chemistry. 
What is sorely lacking is an experimental database that 
will allow for quantification of the observational ef- 
fects in terms of the processes involved. For example, 
there is still no consensus on the degree to which sput- 
ter deposition and vapor deposition (of material vapor- 
ized by micrometeorite impacts) have contributed to 
the deposited material observed on lunar soil grains. 
Similarly, it is clear from recent debates in the litera- 
ture that certain minerals in the lunar regolith (ilmenite 
for example) behave differently than others (most sili- 
cates). Even within silicates, there is the open question 
of how Fe-bearing silicates undergo space weathering 
as compared to Fe-poor or Fe-free silicates. Finally, it 
is generally recognized now that the nanophase and 
submicroscopic Fe metal beads in melt glass, and in 
grain rims are the primary cause of the optical effects 
of space weathering. However, little is known about 
the specific processes involved - Is a hydrogen satu- 
rated surface required? Does sputtering produce sig- 
nificant nanophase Fe? What are the kinetics of these 
processes, i.e. how long does it take to form a vapor- 
deposited rim 50 nm thick? Can the experience we 
have gained from studying space weathering effects in 
the lunar regolith be extrapolated to understanding 
space weathering in asteroidal regoliths [e.g. 16]? 
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Figure 1 . High-resolution transmission electron micro- 
scope image of a nanophase (-7 nm dia.) Fe metal 
grain embedded in agglutinitic glass from 79221. 
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Introduction: Planetary targets have been ob- 

served with radar since the late 1050’s when it was 
first used for ranging experiments with the Moon [1]. 
As telescope size and power increased, it became pos- 
sible to observe more distance targets (Venus, Mars, 
and the outer satellites). Inherent to radar observations 
is the uncertainty as to the source of the reflection, 
there being two points where range and Doppler rings 
intersect on a sphere. The use of interferometric 
methods, first used on the moon with two stations [2] 
and later on Venus and Mars [3], solved this problem. 
We extend the method through the addition of a fourth 
receiving telescope (thus doubling the number of pro- 
jected baselines) and integration of the newly available 
Mars Orbiter Laser Altimeter (MOLA) topographic 
datasets. 

Observations, 2001: During the two and a half 
months surrounding the 2001 Mar> opposition, four- 
teen nights of observations were conducted using four 
telescopes of the Goldstone Deep Space Communica- 
tion Complex (GSDCC). The observations sampled a 
significant swath of the Martian equatorial region, cov- 
erage dependant on planetary geometry, sampling a 
diverse collection of geologic settings and conditions. 
For the purpose of developing and testing the four 
station interferometric technique, we have focused our 
efforts on a single night’s data, that of July 14, 2001. 
This night’s data has the advantages of a high signal to 
noise ratio on all four telescopes and a strong return 
from the moderately rough surface of Chryse Planitia. 
The extent of coverage on the Martian surface is de- 
pendant on the geometry of the Earth-Mars system, 
which determines the longitude range, and the distance 
to the first null phase point, covering approximately 
three degrees of latitude north and south. 

Method: There are two phases to the processing 
of the data. The first phase prepares the data for further 
processing by calibrating the raw' data stream and con- 
verting it from the time domain to range/delay-doppler 
space. The details of this phase are well known and 
covered in more detail elsewhere [4]. Figure 1 is an 
example of the data at this stage of the processing and 
shows the relative signal to noise ratios between the 
four stations. We are primarily interested in the sec- 
ond phase of the processing, the interferometric 
method and modeling of the reflection coefficient. 
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Figure 1: A typical set of range/delay-doppler images 
for July 14, 2001, approximately 10° East longitude 
and 7° North latitude. The frames have been similarly 
stretched to highlight the differences in the signal to 
noise between the four stations. 


Past planetary interferometric radar experi- 
ments have successfully used data from two or three 
stations to map the surface reflectivity (or radar back- 
scatter coefficient) [3]. The addition of a forth receiv- 
ing telescope doubles the number of projected inter- 
ferometric baselines. The equations for the observed 
signals (voltage) at each of the four telescopes at each 
point on the surface of Mars are as follows: 

= J/v + fs + n \ (!) 

V 2 = V N e iuh ' + V s e iuh ‘ + n 2 ( 2 ) 

V 3 =V N e ivh ’ +V s e ivh ’ +« 3 (3) 
V 4 =V N e^ +V s e‘ wh ' +n, (4) 

Where V x is the complex voltage measured at 
each telescope, n x is a random gaussian noise, V NS is 
the backscatter coefficient due to the northern and 
southern points respectively, u, v, and w are spatial 
fringe frequencies due to the geometry of the stations 
and planets, and h N S are the local heights of the north- 
ern and southern hemispheres. 
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Figure 2: Fringe phase image for the DSS14-13 station 
pair. The cross-power image is colored such that the 
pixel brightness represents the power and the color the 
phase. Only the leading few pixels are easily viewable 
due to the predominance of the signal contained in the 
leading edge. 

Since there is an independent gaussian noise 
signal inherent to every measurement, no single solu- 
tion exists and the statistical solution is necessary. The 
iterative maximum likelihood function method 
(MLFM) has been used in the past and found to give 
the best results [5]. Equations 1-4 are represented such 
that one station is chosen as a phase-normalized refer- 
ence station (Equation 1), the phases are assumed to be 
zero. The phases at the remaining three stations are 
then established to be relative to that of the first. 

Furthermore, since the measured complex 
voltages at each station are zero-mean Gaussian ran- 
dom variables, as they are composed of the reflection 
from a large number of randomly oriented sub- 
resolution facets within each delay-doppler cell, the 
square of their magnitude is proportional to the back- 
scattered power, P NiS . The formulation of the maxi- 
mum likelihood estimator is as follows. 


Here, the variables are the same as for equa- 
tion (5) with the addition of a 2 which is the noise com- 
ponent of the signal. Since the terms on the right all 
contain R N and R s , the coupled equations must be in- 
terated to a stable solution. 

Unlike past experiments where, in addition to 
iterating the maximum likelihood solution for the re- 
flectivities, the altitude of the reflecting pixels also an 
unknown, we now possess an extremely accurate to- 
pographic map courtesy of MOLA. We use the 
MOLA topography as an input into the solution, thus 
eliminating the two height variables and simplifying 
the calculation. 

The actual application of the MOLA topogra- 
phy in the radar solution is complicated by both the 
size of the resolution cell (the observed reflected 
power is due to reflection from a large number of 
small facets) and the antenna response function [6]. 
Thus, the target response function (a combination of 
the altitude and orientation of facets within a resolu- 
tion cell) is convolved with the antenna ambiguity 
function to properly weight the altitude. 

Finally, with the completion of the MLFM 
processing, we have 2-d maps of the surface reflectiv- 
ity across the observation regions. We model the re- 
flectivity using Hagfors’ model, where R is the reflec- 
tivity, C ' 1 ' 2 is the rms slope, p is the fresnel reflectivity, 
and G is the incidence angle. 
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Where there are M stations, P(v) is the cross- 
power spectra between two stations and has a northern 
and southern component. K, A m , B m , C m , E, f, and D 
are from variable changes that, for clarity’s sake, are 
combinations of the power and cross-powers. Finally, 
Rn and R s are the northern and southern reflectivities, 
for which we want to solve. By differentiating with 
respect to P N and P s , setting the results equal to zero, 
and solving for R N and R s we arrive at the following 
coupled equations. 


_ A m f 2 + B lrl D + C m F 


( 6 ) 


R =~-(cos 4 9+ Csin 2 d) V2 (8 ) 
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Introduction: The spectral reflectance properties 
of the Moon are governed by the minerals and glasses 
composing the lunar regolith, their physical state, and 
the optical effects of soil maturity. Work by Hapke 
and others has provided all the tools necessary to pro- 
duce a model of near-IR spectra of the Moon within 
simplifying assumptions. We have produced such a 
model and are beginning to apply b to lunar science 
problems. Two classes of problems are amenable to 
immediate use: determining compositions of lunar 

surface regions using groundbased and Clementine 
data, and understanding the detection limits for miner- 
als and rock types using existing and planned data sets. 

Model: Our model is based on the equations of 
Hapke [1] who showed how the visible and near-IR 
spectra of mixtures of minerals could be computed 
from their optical constants at arbitrary grain sizes and 
relative abundances, and recently it uas shown how the 
method of [1] could be modified to include the effects 
of submicroscopic particles [2], 

Optical constants of all components are required 
for this model. An inversion of the Hapke equations 
was presented by [3] in order to compute optical con- 
stants from reflectance spectra. We used this inversion 
to compute the optical constants oi anorthite from a 
spectrum of the mineral in the USGS Denver spectrum 
library. For the optical constants of olivine we used the 
the optical constants computed by []]. For pyroxene 
we applied a new methodology wherein we fit the opti- 
cal constant spectra of a series of pyroxenes with Gaus- 
sians, then regressed the Gaussian par ameters on com- 
position, ignoring Gaussians we believed due to water 
or iron oxide contaminants. We simulated the spec- 
trum of shocked anorthite by interpolating across 
the 1250 ferrous iron band of the USGS plagioclase. 
For native iron, no measurements of these constants 
exist at the spectral sampling of lunar telescopic meas- 
urements, but from 600 to 2000 nm the measurements 
of [4] Johnson and Christy (1974) scatter randomly 
about straight lines for both n and k, so we used linear 
fits to the Johnson and Christy data to represent the 
optical constants of Fe. 

In the forward implementation of the model, the 
chemistry of minerals, grain sizes, modal abundances 
and abundance of submicroscopic iron are defined. 
From the mineral and glass chemistries optical con- 
stants are computed, which are then modified by sub- 
microscopic iron coatings. These modified optical 
constants are then converted to single scattering albedo 
and mixed according to their modal abundances. 


Results: Qualitatively, the model produces spectra 
which closely mimic the appearance of lunar spectra, 
sharing albedo, continuum slope, and spectral contrast, 
as well as the shape of the absorption features. Valida- 
tion of the forward model is proceeding use the spectra 
and analyses of [5]. The validation will determine, for 
example, when the model prescribes a certain grain 
size, or grain size distribution, how this optical grain 
size corresponds to that measured via sieving. 

An important unresolved measurement involves 
submicroscopic iron. Hapke showed that the optical 
effects of iron extend far beyond the range detected by 
ferromagnetic resonance. Work by Morris shows that 
the correlation between various small grain sizes of 
iron is relatively poor. Thus, measurements of Is/FeO 
can only weakly constrain the amount of submicro- 
scopic iron prescribed by our model to best fit lunar 
spectra. Future lunar soil validation measurements 
should measure a wider, or additional, iron size ranges. 

While the validation dprocess is proceeding, we 
have applied the model to lunar science problems that 
do not depend critically on absolute abundances. For 
example, we recently showed that areas on the Moon 
previously interpreted to be composed of anorthosite 
that had its 1.25 micron band erased by shock, can 
plausibly be mimicked by spectra of mature crystalline 
anorthosite [6]. 

A recently introduced application of this model is 
to determine detection limits for minerals as a function 
of soil maturity and data quality (signal to noise ratio). 
In this process we model a particular soil composition, 
ensuring that the model spectrum lies within the field 
of measured lunar spectral properties. We then vary 
modal and chemical properties to determine the magni- 
tude of the differential signal that can be detected by a 
remote sensor. Our preliminary results suggest that 
10% differences in mineralogy can only be detected for 
the most immature surfaces at the 1% precision of 
Clementine and groundbased data. Future lunar mis- 
sions should features sensors with much higher signal 
to noise ratios. 
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LPI Contribution No. 1129 47 


COHERENT BACKSCATTERING BY RANDOM PARTICULATE MEDIA IN THE SOLAR SYSTEM. K. Muinonen. 
Observatory, University of Helsinki, P.<>. Box 14, FIN-00014 U. Helsinki, Finland (Karri.Muinonen@Hel.sinki.Fi). 


Photometric observations for large numbers of asteroids 
indicate an opposition effect, a nonlinear increase of bright- 
ness at small solar phase angles, the angle between the Sun 
and the observer as seen from the asieroid. Vast sets of as- 
teroid polarimetric observations shov* negative polarization, 
a peculiar degree of linear polarization (I± - I\\)f{Ix + /;|) 
for unpolarized incident sunlight: at small phase angles, the 
brightness component Jp with the electric vector parallel to 
the scattering plane defined by the Sun, the asteroid, and the 
observer predominates over the perpendicular component /j_. 
For recent reviews, see Muinonen etai (2002ab). 



The coherent backscattering mechanism (CBM) for the 
opposition effect is described in Fig. la for second-order scat- 
tering. An incident electromagnetic plane wave (solid and 
dashed lines; wavelength A and wavenumber k = 2tt/A) in- 
teracts with two scatterers A and B, which are of the order of 
the wavelength to hundreds of wavelengths apart, and prop- 
agates to the observer to the left. The two scattered wave 
components due to the two opposite propagation directions 
between the scatterers interfere constructively in the conical 
directions defined by rotating the light source direction S about 
the axis L joining the two end scatterers. We illustrate a scat- 
tering direction on the cone precisely opposite to the light 
source direction. Thus, the exact backward direction (phase 
angle a = 0°) is on the constructive-interference cone for 
arbitrary locations of the two scatterers whereas, in other di- 
rections, interference varies from constructive to destructive. 
Three-dimensional averaging over scatterer locations results in 
a backscattering enhancement with decreasing angular width 
for increasing order of interactions, because the average dis- 
tance between the end scatterers is larger for higher orders of 
interactions. The scattering processes can be caused by any 


disorder or irregularity in the medium. 



b. 


Figure 1 : Coherent backscattering mechanism for (a) the op- 
position effect and (b) the negative polarization. See text. 

CBM for the negative degree of linear polarization is ex- 
plained for second-order scattering in Fig. lb. The incident ra- 
diation is unpolarized by definition, which requires the deriva- 
tion and proper averaging of the Stokes vectors corresponding 
to the scattered electromagnetic fields (E s ) for two linear po- 
larization states of an incident plane wave ( Ei ). In Fig. lb, 
incident polarizations parallel and perpendicular to the scatter- 
ing plane (here yz-plane) are treated in the two leftmost and 
two rightmost panels, respectively. 

Consequently, in Fig. lb, an incident electromagnetic plane 
wave interacts with two scatterers A and B at a distance d from 
one another aligned either on the x-axis or the y-axis. while 
the observer is in the yz-plane. For the present geometries, 
the constructive interference cones of Fig. la reduce to the 
yz and xz-planes, depending on the alignment of the scatter- 
ers. Since first-order scattering is typically positively polarized 
(e.g., Rayleigh scattering and Fresnel reflection), the scatter- 
ers sufficiently far away from each other (kd = 2izd/ A » 1) 
interact predominantly with the electric field vector perpendic- 
ular to the plane defined by the source and the scatterers (two 
upper panels), while interaction with the electric field vector 
parallel to that plane is suppressed (two lower panels). The 
observer in the yz-plane will detect negative polarization from 
the geometry in the upper left panel of Fig. 1, and positive 
polarization from the geometry in the upper nght panel. How- 
ever, the positive polarization suffers from the phase difference 
kd sin a, whereas the phase difference for the negative polar- 
ization is zero for all phase angles. Averaging over scatterer 
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locations will result in negative polarization near the backward 
direction. Scattering orders higher than the second experience 
similar preferential interaction geometries, and contribute to 
negative polarization. As above for the opposition effect, the 
contributions from increasing orders of scattering manifest 
themselves at decreasing phase angles. 

Figure 2 shows tentative modeling of the opposition effect 
and negative polarization of C, M, S, V, and E-class asteroids 
using Monte Carlo computations for coherent backscattering 
(Muinonen 2002; Muinonen et al. 2002ab). The novel nu- 
merical technique relies on the reciprocity principle for elec- 
tromagnetic scattering, allowing a renormalization of the co- 
herent backscattering contributions at each multiple-scattering 
event. 
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Coherent backscattering is here computed for a plane wave 
(wavelength A = 27r/ 10 /xm) normally incident on a semi- 
infinite medium of discrete Rayleigh scatterers with single- 
scattering albedos of 0.9 and medium mean-free-path length 
of £ — 1 , 2, 3, . . . , 50 /im. For Rayleigh scatterers, the single- 
scattering albedo chosen roughly results in a maximum op- 
position effect amplitude as well as a maximum depth of the 
negative polarization surge. 5 x 10 4 rays were used in each 
computation. 

In Fig. 2, an ad hoc exponential weighting scheme of 
exp (~£/L) has been applied to the numerical results in all 
cases with L = 1.0, 2.0, 2.0, 2.0, and 6.0 /im for the C, M, 
S, V, and E -cl ass curves, respectively. The weighting scheme 
derives from the idea that, rather than modeling the random 
scattering medium using a single mean free path only, we allow 


a distribution of mean free paths and thus physical character- 
istics on asteroid surfaces. Additionally, the CBM computa- 
tions have been multiplied by a linear function 1 — Ka, where 
K = 0.021, 0.017, 0.015, 0.013, and 0.009 /° forC, M, S, V, 
and E-class curves, respectively. For comparison with the po- 
larimetric observations, the negative polarization surges have 
been divided by the scaling “fudge factors” of 1.5, 2.3, 3.0, 
4.0, and 7.2. The numerous implications of the tentative fits 
are subject to future discussion. 
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Figure 2: Opposition effect (a) and negative linear polarization 
observations (b) for C, M, S, V, and E-class asteroids with 
heuristic modeling. For the references of the observations, see 
Muinonen et al. (2002a). 
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Introduction: The Compact Reconnaissance Imag- 
ing Spectrometer for Mars (CRISM) on the Mars Re- 
connaissance Orbiter (MRO) will conduct a compre- 
hensive series of investigations of the Martian surface 
and atmosphere. The investigations will be accom- 
plished using an instrument design that provides high 
spatial and spectral resolutions, extended wavelength 
range, and ability to gimbal through a range of orienta- 
tions. Baseline investigations include a near-global 
survey to find high science priority sites, full-resolution 
measurement of thousands of such sites, and tracking 
of seasonal variations in atmospheric and surface prop- 
erties. 

Science Overview: The Mars Exploration Payload 
Analysis Group [1] recommended specific hyperspec- 
tral imaging investigations to characterize Martian ge- 
ology, climate, and environments of present or past 
life. CRISM's three groups of investigations address all 
of MEPAG’s recommendations. The first two groups 
correspond to two primary objectives of MRO: to 
search for evidence of aqueous and/or hydrothermal 
activity ; and to map and characterize the composition , 
geology, and stratigraphy of surface features. These 
investigations are implemented by high-resolution hy- 
perspectral mapping of thousands of high priority tar- 
gets including candidate sedimentary deposits [2], vol- 
canic regions, crustal sections exposed in steep es- 
carpments, and sites which exhibit evidence in Mars 
Express/Omega data for concentrations of aqueously 
formed minerals. The third group of investigations ad- 
dresses the primary MRO objective to characterize 
seasonal variations in dust and ice aerosols and water 
content of surface materials , and the secondary objec- 
tive to provide information on the atmosphere com- 
plementary to other MRO instruments. These investi- 
gations are implemented using a systematic, global grid 
of measurements of the emission phase function (EPF) 
acquired repetitively throughout the Martian year. EPF 
measurements allow accurate determination of column 
abundances of water vapor, CO, dust and ice aerosols, 
and their seasonal variations [3]. At the same time, the 
grid's repetitive coverage will track seasonal variations 
in water content of surface materia!. Additional, tar- 


geted observations of the polar caps will investigate 
their inventory of water and C0 2 ices. 

When not taking targeted measurements, CRISM 
will conduct a -100 m/pixel, ~59- wavelength survey to 
search for evidence of aqueous activity that lacks mor- 
phologic expression and/or is below the resolution of 
previous spectral mapping. The survey addresses 
MRO’s secondary objective to identify new sites with 
high science potential for future investigation , and will 
be particularly important for identification of key Noa- 
chian deposits. For example, even at low Martian ero- 
sion rates [4], morphologic expressions of Noachian 
hot spring deposits would have been removed in ~10 9 
yrs, so that such deposits may now exist only as miner- 
alized spots in morphologically unremarkable eroded 
escarpments, crater ejecta, and talus. Much of the sur- 
vey will be completed before MRO’s highest downlink 
rates, so that newly discovered sites can be targeted 
with full-resolution coverage. 

Instrument Overview: An overview of the instru- 
ment is shown in Figure 1 . CRISM's two major subas- 
semblies are the Optical Sensor Unit (OSU) and the 
Data Processing Unit (DPU). The OSU contains visible 
(VIS) and infrared (IR) imaging spectrographs that 
share the same 2.05° field-of-view and cover the wave- 
length ranges 400-1050 run and 1050-4050 at 7 
nm/channel. The ER focal plane is cooled with doubly 
redundant cryogenic coolers. A side-facing radiator 
cools the spectrometer cavity to decrease instrument 
background. An onboard integrating sphere provides 
radiance calibration, and a shutter and dedicated dark 
pixels on each focal plane provide background meas- 
urements interleaved with Mars measurements. Optics 
and the focal planes are enclosed in a housing that is 
gimballed by a high-precision motor/encoder ±60° 
along-track from nadir. During measurement of a tar- 
get, along-track scanning takes out most groundtrack 
speed and allows long integration times and high sig- 
nal-to-noise ratio (SNR) data to be obtained. At the 
same time, scanning provides the capability to remeas- 
ure a spot repeatedly at multiple emission angles to 
characterize the EPF. High spatial resolution (—18 
m/pixel from 300 km altitude) allows characterization 
of the surface at the outcrop scale, and wavelength 
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coverage to >4000 nm provides sensitivity to carbon- 
ates even at low abundances. The DPU provides loss- 
less or lossy compression of the data stream in real 
time. Key design elements are adapted from the 
CONTOUR and MESSENGER optical instruments. 

Operations Overview: Figure 2 summarizes 
CRISM's three complementary data acquisition strate- 
gies, which are built around the three steps needed to 
characterize high scientific priority, small-scale depos- 
its: find the deposits, separate their signature from that 
of the atmosphere, and acquire high spectral and spatial 
resolution measurements with high SNR. Most of the 
time the OSU views nadir, building up coverage in 
multispectrai survey mode. In the second mode of ob- 
serving, for a group of ~36 orbits once every -30° of 
Ls, EPFs are measured on a 10° longitude x 15° lati- 
tude grid. The third mode is targeted observations, 
approximately 2200 of which will be obtained at key 
locations. At each site, gimballing is used to take out 
most along-track motion, so that the field-of-view is 
slowly scanned over a rectangular swath approximately 
11 km wide (when measured from 300 km altitude). 
Measurements of the swath are bracketed by incoming 
and outgoing measurements of the EPF of the center- 
point of the swath at emission angles up to ±60°. Lar- 
ger high-priority targets will be identified in Mars Ex- 
press/Omega data and other data sets. The multispec- 


tral survey provides redundancy to Omega data in lo- 
cating targets that have VIS-IR spectral signatures but 
lack obvious morphologic expressions, and it provides 
the spatial resolution needed to identify small targets 
not evident in Omega data. 

References: [11 MEPAG, Mars Exploration Pro- 
gram: Scientific Goals, Objectives , and Priorities , 
2000. [2] Cabrol, N. and E. Grin, Icarus , 742, 160- 
172, 1999. [3] Clancy, R.T., and S. Lee, Icarus , 93 , 
135-158, 1991. [4] Presley, M., and R. Arvidson, 
Icarus , 75, 499-517, 1988. 



Fig. 1 . Preliminary CAD renderings of CRISM's OSU, 
DPU, and control electronics. 


Measurements of the 
emission phase function 
at each target and as a 
global grid every $0 days 
track seasonal changes in 
aerosol properties, CO and 
water vapor abundances, 
and the hydration state 
of the surface. 


Gimballing takes out along-track motion of 
the field-of-view, allowing longer integration 
times and higher SNR. 


Targets are imaged at 24 m/pixel 
(from 400 km). High spectral resolution 
(7 nm/channel) and coverage of carbonate 
fundamentals at 3800-4000 nm allow low 
concentrations of aqueous minerals 
In small outcrops to be identified. 



When not measuring targeted 
sites, CRISM will survey Mars at 
~100 m/pixel In 50 key wavelengths 
to discover new high-priority sites 
below the resolution of previous data. 


Fig. 2. Schematic depiction of key attributes of CRISM's measurement plan for the Martian surface and atmosphere. 
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Introduction: Hapke s [1] radiative transfer 

model for particulate surfaces has been demonstrated 
to be accurate for describing the reflectance of surfaces 
where the particles are larger than the wavelength of 
light (geometric optics regime). However, for particles 
that are smaller than the wavelength of light, the be- 
havior of radiation in a particulate surface is different 
and not well described by geometric optics. Combin- 
ing Mie and Hapke theories has been shown as a po- 
tentially promising avenue in modeling the reflectance 
of fine particle surfaces [2,3] through these studies 
have also shown that important discrepancies between 
models and observations exist. Another important area 
of basic remote sensing research is m the changes that 
occur in reflectance due to cementation [4], Processes 
that occur in soils can lead to cementation of particu- 
lates such as the formation of caliche or hardpans, or in 
the case of evaporite basins. Observations of the re- 
flectance of particulate and cemenled materials show 
that the change due to cementation are systematic and 
distinct from compaction or other effects, but have not 
been modeled [4]. Here I describe experimental ob- 
servations and modeling of these two cases (fine parti- 
cles and cementation) and discuss possible paths to 
resolve these problems. 

Fine Particle Reflectance: Mustard and Hays [3] 
measured the reflectance of particle size separates of 5 
m intervals between 0 and 25 m of the minerals ol i- 
vine and quartz to examine the spectral properties of 
these materials when the particle size is greater than, 
approximately the same size as, and then much smaller 
than the wavelength. They used the wavelength range 
0.3-25 m which, when combined with the size sepa- 
rates used, allowed examination of this processes 
across a range of optical constants expressed by the 
minerals. The expectation was that for a given particle 
size separate at some point, and likely multiple points, 
along this wavelength range the transition from optical 
to Mie scattering should occur. 

The results of these measurements indeed demon- 
strated the expected effects. When the particle size 
approached, and then became smaller, than the wave- 
length of light, the reflectance exhibited a sharp de- 
crease. This is predicted to some extent by the size 
parameter which relates the particle size to the wave- 
length of light. However, the size parameter needs to 
be combined with concepts of extinction efficiency to 
fully explain the observations. Mustard and Hays 
therefore defined a critical diameter, Dc (equation 1, ) 


n=real part of the complex index of refraction), based 
on Mie and scattering theory [1]. When Dc is less than 
2, then the reflectance should experience a sharp de- 
crease. This relationship is only valid for optical re- 
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Figure 1. Comparison of measured and modeled reflectance 
spectra of fine particle separates of olivine 


gimes where the imaginary component of the index of 
refraction is small. Mustard and Hays [3] demonstrated 
that the wavelengths where quartz and olivine should 
exhibit a sharp decline in reflectance on the basis of 
the particle size and the index of refraction corre- 
sponded well with observations. 
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Modeling of these same observations had mixed re- 
sults. A Mie model was used to calculate single- 
scattering albedo and a Hapke [1] model to determine 
reflectance, and inputs to the model were optical con- 
stants and particle size distribution. For olivine (Fig- 
ure 1), the model agrees with observation in terms of 
trends: relative albedo of the separates for the various 
scattering regimes was correct, transitions and cross- 
overs of reflectance occurred at approximately the 
same wavelengths. Although the model and observa- 
tions differ in their details, this could be due to the op- 
tical constants used. For quartz, however, it was quite 
a different story. The model was consistent with the 
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observations only in the most crude sense, and the de- 
tails were significantly different. Here again, the sus- 
picion is that with appropriate optical constants, per- 
haps a better agreement could be attained. 

Cemented Soils: Cooper and Mustard [4] exam- 
ined the effects of cementation on the reflectance 
spectra of Mars analogs. The samples consisted of 
mixtures of sulfate and JSC-Mars 1, a martian soil 
analog derived from Hawaiian palagonite, with pro- 
portions that matched the amount of sulfur detected by 


Loose Powders 
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Figure 2a. Reflectance spectra of JSC Mars-1 and MgS0 4 
mixtures. Note the vertical axis for MgS0 4 is different. 

the Viking and Pathfinder landers. The samples were 
measured as powders, packed powders, and cemented 
(details of these preparations are presented in [5]). 

The spectra of the loose soils vs their cemented 
counterparts were largely unchanged in the visible to 
near infrared. However, in the thermal IR (7-20 m) 
significant changes did occur (Figure 2ab). The sulfate 
features around 9 m were barely identifiable in the 
powdered sample but became greatly enhanced in the 
cemented soils. This was understood from a physical 
perspective. The optical constants for sulfate become 
large at the 9 m features, creating a restrahlen band. 
When in the powdered form, the restrahlen band is 
weak due to multiple scattering. However, after ce- 
mentation, the radiation sees a more coherent surface 
and multiple scattering is reduced creating an enhanced 
restrahlen feature compared to the iso-chemical pow- 
dered counterpart. It is interesting to note, however, 
that the JSC-Mars 1 sample also exhibits strong 
changes in the spectral properties in the thermal IR, 


though the cause of these differences were not ex- 
plored by [5]. 

While these data showed drammatic changes that 
are readily understood, these effects have not yet been 
modeled. This could have much application for under- 
standing Mars as the new generation of high spatial 
and spectral data sets become available. It appears that 
the combination of visible-NIR with thermal IR data 
could provide compositional and textural information 
if properly modeled. 


Cemented Mixtures 
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Figure 2b. Same samples as in 3a but after cementation as 
described in text. Note change in vertical axis for mixtures 
relative to 2a and b. 

Challenges: Laboratory observations of very fine 
particle separates indicate that theory provides a basis 
for understanding the observations, but there are sig- 
nificant problems in applying quantitative models. 
Similarly, theory provides a sound foundation for un- 
derstanding the effects of cementation, but has not yet 
been applied or developed to treat this texture. Fine 
particles and cementation are two important textures 
expected on Mars and these effects will likely be pre- 
sent in new data sets to be collected. Thus it will be 
important to develop a better modeling base to analyze 
these problems. 

References: [l]°Hapke, B. 1993. Theory of Reflectance 
and Emittance Spectroscopy , Cambridge University Press, 
Cambridge. [2]‘Moersch, J. E and P. R. Chnstensen (1995) 
J. Geophys. Res., 100 , 7465-7477 [3]*Mustard, J. F. and J. E. 
Hays, (1997), Icarus , 125, 145-163. [4] Cooper, C. and J. F. 
Mustard (2002) Spectroscopy of loose and cemented sulfate- 
bearing soils: Implications for duncrust on Mars, Icarus, (in 
press). 
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The reflection of electromagnetic radiation from a 
planetary regolith involves a combination of geometric 
and physical optics processes which contribute to the 
signal returned to the remote observer. The geometric 
optics effects are the product of smgly and multiply 
scattered radiation from the surfaces of the regolith 
particles, combined with radiation which has under- 
gone various combinations of transmission through one 
or more regolith grains followed by one or more scat- 
terings from other particles. The physical optics ef- 
fects include diffraction of radiation around the edges 
of large irregular particles and cooperative coherent 
scattering between particles which are small when 
compared to the wavelength of the incident radiation. 
These effects produce measurable changes in the inten- 
sity and polarization of reflected light as a function of 
illumination and viewing geometry In particular, as 
phase angle becomes small, the reflectance of a par- 
ticulate material will increase non-lmearly and exhibit 
the ‘opposition effect’ [1]. In the planetary science 
context, the phase curve, and in particular the size of 
the opposition surge and the width of the phase curve 
near zero degrees, have been attributed to two proc- 
esses commonly called ‘shadow hiding 5 (SHOE) and 
‘coherent backscattering 5 (CBOE) 1 2,3,4] Understand- 
ing the contribution of the SHOE .md CBOE compo- 
nents to the integrated phase curve will permit models 
to be developed which can confidently determine the 
fundamental regolith textural properties such as parti- 
cle size and packing density from remote sensing data. 

Laboratory experiments have attempted to distin- 
guish the contribution of SHOE and CBOE by present- 
ing samples with circularly polarized monochromatic 
electromagnetic radiation and measuring the circular 
polarization ratio in the returned signal. If the returned 
signal is singly scattered then the ciu ularization polari- 
zation ratio should decrease as phase angle decreases. 
However, if the returned signal is multiply scattered 
then the circular polarization ratio should strongly in- 
crease as phase angle decreases. We have observed this 
increase in circular polarization rauo with decreasing 
phase angle in many highly reflective particulate media 
[5,6,7]. 

In the case of materials with low reflectance the ex- 
pectation is that most of the returned signal is singly 
scattered because of the high probability of a photon 
being absorbed at each scattering. Hence, for absorbing 
media the circular polarization ratio is not expected to 


sharply increase as phase angle decreases. In general 
we have found this to be true [5]. However, recently 
we have encountered an interesting counter example. 
Measurements on a suite of boron carbide samples 
(reflectance =~ 5 %) have found a significant increase 
in circular polarization ratio with decreasing phase 
angle, a result that is not consistent with our interpreta- 
tion of the process. This result suggests that albedo 
alone is not the principal regulator of the amount of 
multiple scattering in the medium. The explanation for 
this unusual behavior remains unknown; a conjecture 
might be that a unique particle shape may create a very 
unusual single scattering phase function in boron car- 
bide. 

It is important to understand this unusual behavior 
before we can proceed with development of models 
that can accurately retrieve textural properties from 
remote sensing data. 

This work performed at JPL under a contract from 
NASA’s Planetary Geology and Geophysics Program. 
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Introduction: Like the Moon, Mercury has no atmos- 
phere to protect it from the harsh space environment and 
therefore it is expected that it will incut the effects of space 
weathering [1]. These weathering processes are capable of 
both creating regolith and altering its optical properties 
[1,2,3]. However, there are many important differences be- 
tween the environments of Mercury and the Moon. These 
environmental differences will almost certainly affect the 
weathering processes as well as the products of those proc- 
esses. It should be possible to observe the effects of these 
differences in Vis/NIR spectra of the type expected to be 
returned by MESSENGER. More importantly, understand- 
ing these weathering processes and their consequences is 
essential for evaluating the spectral data returned from 
MESSENGER and other missions in order to determine the 
mineralogy and the iron content of the Mercurian surface. 
Theoretical and experimental work has been undertaken in 
order to better understand these consequences [4,5]. 

Mercurian Environment: Due to its higher impactor 
flux and greater density. Mercury will produce 13.5 times the 
melt and 19.5 times the vapor than is produced on the Moon 
per unit area [6]. Additionally, Mercury has a magnetic field 
that protects its surface from charged particles which se- 
verely reduces the solar wind flux at the surface [7]. There- 
fore, melting and vaporization due to micrometeorites are 
expected to dominate space weathering on Mercury with 
little or no solar wind sputtering effects [1]. Furthermore, 
agglutinitic glass-like deposits and vapor deposited coatings 
should be created faster and more efficiently on Mercury. 
The Mercurian environment is also notable for its extreme 
temperature range. Due to its slow rotation and proximity to 
the sun, equatorial regions of Mercury achieve temperatures 
above 700K during the day, while nighttime temperatures 
fall below 100K. 

The nanometer-scale metallic Fe particles (npFe 0 ) that 
are ubiquitous in the rims and agglutinates of lunar soil [8] 
(figure 1) should also be present on Mercury. Even for the 
endmember case where the surface of Mercury has no native 
FeO, the iron brought in by meteorites should be sufficient to 
make the formation of npFe 0 through \apor fractionation an 
important process on the planet. Amounts as small as 0.05 
wt % npFe 0 are enough to affect the optical properties [2]. 



Figure 1. TEM bright field image of npFe°-rich rim on a 
lunar soil grain. 


Size of npFe 0 Particles: The size distribution of metallic 
Fe particles in a soil strongly controls the optical effects on 
the Vis/NIR spectrum. In lunar rims, these npFe 0 particles 
range from ~l-12nm in dia. with an average of ~3nm [9]. 
The smallest particles (3-5nm) will tend to redden the soil 


while larger particles (>10nm) will simply cause darkening 
[9,10]. 

Ostwald Ripening: In figure 2 are shown the results of a 
npFe 0 study in which Fe° particles averaging about 8 nm in 
diameter were created and then heated for 10 hours at a range 
of temperatures [12]. The graph plots the final size of the 
particles vs. the heating temperature. This study is not di- 
rectly applicable to Mercury because here iron particles are 
in grain to grain contact rather than occurring as isolated 
particles suspended in a glass matrix as in space weathering 
products. However, the experiment is useful because it dem- 
onstrates that this size of nanophase Fe° particles may only 
be stable to about 200°C. The majority of Mercury’s surface 
reaches daytime temperatures significantly above 200° C. In 
fact, during the course of a Mercurian day, the soil at the 
hottest parts of Mercury will stay above 400°C for about 2 
weeks. This increased temperature may be enough to allow 
the npFe 0 particles to grow significantly via Ostwald ripen- 
ing, a process well known in material sciences [1 1]. Ostwald 
ripening is a process by which npFe 0 particles in a glass ma- 
trix will tend to coarsen due to differences in free energy 
between curved surfaces. 



Figure 2. Size of npFe 0 particles after annealing for 10 
hours. Modified from Glieter, 1989 [12]. The shaded areas 
show the upper extent of the temperature regimes for the 
Moon and Mercury, respectively. 

A vapor deposition experiment of Hapke et al. [13] dem- 
onstrated that heating npFe°-rich vapor coatings to a tempera- 
ture of 650°C for just one hour is sufficient to remove a fer- 
romagnetic resonance. This is hypothesized to occur because 
particles of npFe 0 have grown to be larger than the range that 
is measured by FMR techniques [4 - 33 nm in dia.[14]]. 
Thus, it appears that even for npFe 15 particles suspended in a 
glass matrix, increased temperatures can result in significant 
grain growth. 

Determining the rate of Ostwald ripening on Mercury is 
difficult due to a lack of directly relevant experimental data. 
The least constrained, and most important, variables are the 
diffusion coefficient of Fe° into the glass matnx, which is 
strongly temperature dependent, and the surface energy (i.e. 
boundary between npFe 0 and matrix), also somewhat tem- 
perature dependent. By considering a wide range of values 
for these, we have attempted to bound the possible range of 
grain growth through time that might occur on Mercury (fig- 
ure 3). 
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Figure 3. Possible range of effects of Ostwald ripening with 
time near the equator of Mercury. Lower bound assumes D 
(diffusion coef.) =10' 19 m 2 /s, cr (surface energy) =0.1 mj/m 2 . 
Upper bound assumes D=10~ 17 m 2 /s, a =100 mJ/m 2 . 


Obviously, these results cover a wide range of possibili- 
ties, however, even our most conservative estimates (lower 
bound) indicate that Ostwald ripening should have a signifi- 
cant effect on equitorial Mercurian soils, doubling the size of 
the npFe 0 in a matter of centuries. Of course, with increasing 
latitude, less solar heating occurs and the thermal regime 
becomes much more lunar-like where Ostwald ripening will 
have little to no effect. Certainly polar regions that do not 
reach temperatures above 200° C (beyond 75°N or S for the 
hot poles, and 65° for the warm poles) are not expected to be 
affected by Ostwald ripening as npFe 0 particles appear to be 
stable at those temperatures [12]. 

Discussion: The size of npFe 0 particles will be reflected 
in remotely acquired data. Ostwald ripening should result in 
larger Fe° particles, on average, near the equator where the 
highest temperatures are reached. As noted earlier, small 
npFe 0 particles cause reddening of the reflectance spectrum 
and larger particles result in darkening. If Ostwald ripening 
dominates over npFe 0 production near the equator, we expect 
the spectral continuum to become darker at lower latitudes as 
one approaches the equator and the increased heat allows for 
larger npFe 0 particles. The continuum should be reddest at 
high latitudes, where Ostwald ripening has little or no effect 
and npFe 0 particles remain small. 

Our current spectral data set for Mercury, is limited. 
Most of our spectral data is telescopic [ 1 5], largely providing 
an integrated disk view, masking possible latitudinal varia- 
tions, as well as regional differences. Also hidden are matur- 
ity differences that might be expected at young craters. Two 
bands of spectral data taken during the Manner 10 flyby 
confirms that spectral differences exist on regional and local 
scales [16], however this dataset does not have the spectral or 
spatial resolution necessary to see the effects discussed 
above. Recently, the surface was mapped over the wave- 
length range 550-940 nm at roughly 200 km resolution by the 
Swedish Vacuum Solar Telescope [17]. Originally, latitu- 
dinal variations were removed in the calibration of this 
dataset, however, some of the data has been reprocessed to 
allow spectral variations with latitude to be observed [18] 
(figure 4). Despite significant scatter, it does appear that, as 
predicted, there is a positive correlation between spectral 
redness and increasing latitude. 



Absolute Latitude 

Figure 4. Spectral slope from 650 to 850nm vs. absolute 
latitude for 16 features on Mercury. Reprocessed Mercury 
data provided by J. Warell. 

Future Work: To address some of the questions raised 
here more quantitatively, we have begun two sets of experi- 
ments. Both use porous silica gel powders impregnated with 
nanophase metallic iron, following the technique of Allen et 
al [19). The first of these experiments will investigate the 
optical effects of particle size of npFe 0 . Because the pore 
size, and by analogy the npFe 0 panicle size, can be con- 
trolled, these powders are well suited for studying the optical 
properties of differing sizes of npFe 0 particles. We expect to 
derive from this experiment the particle size at which npFe 0 
transitions from causing spectral reddening to darkening. 
The second experiment involves heating the npFe 0 -bearing 
powders for extended periods of time (days to weeks) and 
measuring the growth rate of the npFe 0 particles. This will 
allow us to better constrain the range of expected effects of 
Ostwald ripening. In combination, these experiments will 
provide a foundation for understanding the optical effects of 
npFe 0 in Mercury’s environment. This knowledge, in turn, is 
essential to evaluate mineralogy and iron content from the 
Vis/NIR spectra expected to be obtained by MESSENGER. 
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Introduction: Several techniques have been used 
to estimate the opposition spike parameters of the 
lunar surface by Clementine data [1 3]. The simplest 
way is averaging many images that contain the space- 
craft shadow point (zero-phase-angle point) neglect- 
ing albedo variations, which anyhow remain after the 
averaging [1]. More accurate way is to use a phase- 
ratio image, one component (one initial image) of 
which contains the zero-phase-angle point and the 
other one is obtained at larger phase angles [2]. This 
technique is more precise as the phase-ratio images 
almost do not contain albedo variations. A drawback 
of both these approaches is strong dependence of 
final results on camera characteristics drift that was 
considerable during the mission. The third technique 
suggests to use images acquired almost simultane- 
ously and containing the zero-phase-angle point [3]. 
It allows to avoid problems with the drift of camera 
parameters. With that technique logarithmic deriva- 
tives of the phase function for each point of imaged 
surface can be calculated. For this purpose D/C color- 
ratio (from UWIS camera images) is used [3]. A 
shortcoming of the technique is availability of slight 
variations of the D/C color-index (D = 0.95 pm and C 
= 0.90 pm) over the ratio image. We estimate here an 
opposition spike parameter for several small mare 
and adjacent highland regions of the lunar surface 
suggesting a new technique (that develops the third 
one), which is also free of the UWIS camera char- 
acteristics drift. We have studied a near-equatorial 
region placed in the latitude range 8°N - 8°S at the 
0.5°E longitude (the lunar near-side) 

Technique and results: The Clementine space- 
craft was moving along the polar orbit and imaging 
the lunar surface near the nadir point. Thus images 
taken successively at the same spectral band overlap 
almost by 30%. The phase-ratio images for the over- 
lapping areas do not burden with slight variations of 
the D/C color-index. We use the overlapping images 
containing the zero-phase-angle point or its neighbor. 
This point is placed at different sites of the images (or 
even beyond the images), not close to each other, as it 
was in case of the D/C color-ratio images. Therefore, 
we may not determine the derivatives of the phase 
function directly. Nevertheless, these phase-ratio im- 
ages are very informative and can be used to estimate 
lunar opposition spike characteristics. 

Initially we followed the algorithm of Clementine 
data calibration described in [4]. We used the radio- 


metric calibration only (including the flat field cor- 
rection). Thus the counts in the images used are pro- 
portional to the reflectivity values. Totally eight pairs 
of UWIS camera images in each filter covered the 
lunar surface region indicated above were used. The 
ratio of two initial images in one filter represents a set 
of quotients F(a 1 )// r (a 2 ) for a set Aa - cq - a 2 , where 
F( a) is the phase function, a being the phase angle. 
Having taken the equation for the phase function in 
the form F( a) - exp(P(a)), where P(a) = k x a + k 2 a 2 
+ . . . , the coefficients k, can be found. For this pur- 
pose we use the experimental quotients and the least- 
squares method. 

Consider the suggested method using as an exam- 
ple a pair of images comprising the north-eastern 
portion of Sinus Medii (North is below). Initial im- 
ages we used is presented in Fig. 1 a and b (the cra- 
ters Bruce and Blagg are seen in Fig. 1 a). The white 
lines show borders of the common portion of the im- 
ages. The ratio of the initial images is shown in Fig. 
2. The levels of the equal values of |Aa| for each im- 
age element are shown in Fig. 3. As it turns out the 
simplest presentation for the phase function, F(a) = 
exp(A:a), is quite enough, since the values |Aa| are 
rather small, they do not exceed a few degree. We 


minimize the sum 


l f(a 2 ) 


-*(c 


in order to 


estimate k. We have estimated the phase function 
slope k for eight pairs of images covering both mare 
surface (Sinus Medii) and highlands southward of the 
mare region. We use images taken in four spectral 
bands. A, B, C, and D. All these provides data to 
study correlation between the phase function slope k 
(at first several degrees of phase angle) and reflectiv- 
ity of the lunar surface. This correlation is shown in 
Fig. 4. As one can see there is no clear dependence of 
k on reflectivity. 

Discussion: The brightness opposition spike of 
regolith-like surfaces is formed by two mechanisms: 
the shadow-hiding mechanism and the mechanism of 
coherent backscatter enhancement. The slope of 
phase angle function formed by the first mechanism 
decreases with growth of albedo, as the shadowing 
effect is suppressed by multiple scattering, whose role 
increases with albedo increase. Unlike this, the coher- 
ent backscatter mechanism predicts increasing the 
slope with growth of albedo. Our results show that 
there is no dependence at all. Three explanations for 
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this can be suggested: (a) the radiometric calibration 
is not perfect enough to obtain reliable estimates of 
photometric slopes; (b) the lunar surface in the region 
under study is too dark (actually albedo varies in the 
range 7 - 1 8%) to provide a detectable demonstration 
of both mechanisms; and (c) as the mechanisms affect 
the photometric slope in the opposite manner, they 
can mutually be compensated. It is difficult to make a 
choice between this three cases. In principle one can 
have surprises with calibration of Clementine UVVIS 
data, but we consider that the probability to obtain a 
“zero’ 1 result due to instrumental errors is very low\ 
The second reason is also unlikely, as almost ten 
years ago B. Hapke [5] showed that some sample of 
the lunar regolith definitely reveal the coherent back- 
scatter effect. The third reason seems to be more at- 
tractive, though it demands studies, in particular, with 
laboratory or computer modeling. 

Anyway, in future it is challenging to use the 
above-described technique for data, which will be 
obtained during ESA Smart- 1 mission to construct 
phase ratio images using overlapping Clementine and 
Smart- 1 frames. 



Fig. 1 a,b. The region under study. Crosses indi- 
cate the zero-phase-angle points. Overlapping regions 
have white outlines. 




Fig. 3. The isolines of the equal |Aa| = a t - a 2 values 
(scale is given in degrees) for the phase ratio image. 
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Fig. 4. Correlation between the phase function 
slope k and albedo of the lunar surface (data for 
UVVIS filters A, B, C, D). 
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Radar is a powerful source of information about the physical and dynamical properties of solar system bodies. 
Radar-detected targets include the Moon, Mercury, Mars, Venus, Phobos, lo, Europa, Ganymede, Callisto, Titan, 
Iapetus, Saturn’s rings, eight comets, and 179 asteroids (75 main-belt and 104 near-Earth). This talk offers a 
perspective on the disc-integrated radar properties of solar system bodies and then turns to what radar remote sensing 
can tell us about asteroids using spatially-resolved measurements. 

In dual-polarization observations, which have been routine for two decades [1], echo power is received in the 
same circular polarization as transmitted (the SC sense) as well as the opposite (OC) sense. SC/OC is a measure of 
the near-surface structural complexity at scales near the wavelength (usually 3.5 cm for Goldstone and 13 cm for 
Arecibo). If SC/OC is near zero, then single reflections from smooth surface facets dominate the scattering process, 
but constraints on the scattering process and physical characteristics become increasingly ambiguous and 
model-dependent as SC/OC increases to a few tens of percent or larger, because such values indicate increasing 
contributions from single reflections from rough surfaces and/or multiple scattering. 

Among all the radar- detected planetary bodies in the solar system, Europa, Ganymede, and Callisto have the 
most unusual radar properties [2,31- Their reflectivities are huge and their circular polarization ratios exceed unity 
(Figure I); that is, in contrast to the situation with other targets, the scattering largely preserves the handedness, or 
helicity, of the transmission. As Hapke [4,5] first realized, the key to understanding the satellites’ radar behavior is 
the coherent backscatter effect. The icy satellites’ echoes are due not to external surface reflections but to volume 
scattering from within the satellites’ regoliths, which the radar sees as an extremely transparent, disordered random 
medium. The high radar transparency of ice compared with that of silicates permits deeper radar sounding, longer 
photon path lengths, and higher order scattering from regolith heterogeneities. Coherent backscattering yields strong 
echoes with SC/OC larger than unit} because the incident wave’s direction is randomized before its helicity is 
randomized and before its power is absorbed. 

The SC/OC estimates for several near-Earth asteroids (NBAs) approach Callisto’s value [6], but their radar 
albedos are not large, so those asteroids’ echoes do not suggest significant coherent backscattering, much less the 
presence of ice. Large values of SC 'OC probably can be explained by many, very different structural configurations. 
On the other hand, the fact that SC/( )C estimates for NBAs range from near zero to near unity establishes that the 
surfaces of these objects are extremely variegated and in some cases much more complex than any seen by spacecraft 
that have landed on Eros (whose SC/OC = 0.3, near the NEA average [7]), the Moon, Venus, or Mars. 

For most asteroids, radar is the only groundbased remote sensing technique that can achieve useful spatial 
resolution. Delay-Doppler imaging with ~20-m resolution is now routine whenever the echoes are strong enough, 
and images with adequate orientational coverage can be used to make three-dimensional models. As reviewed by [8] 
radar has revealed both stony and metallic objects, principal-axis and complex rotators, very smooth and extremely 
rough surfaces, objects that must be monolithic and objects that probably are not, spheroids and highly elongated 
shapes, contact-binary shapes, and binary systems (which apparently constitute - 1/6 of NEAs as large as 0.2 km [9]). 

Recently it has been shown that knowledge of physical properties is tightly coupled to long-term prediction of 
close approaches of NEAs to Earth, because of “Yarkovsky” accelerations from thermal re -radiation of absorbed 
solar energy[10]. Someday our ability to ascertain whether some asteroids could threaten Earth centuries in the 
future may rest on how well we understand their global distributions of photometric and thermal properties. 
References. Ostro S. J. (2002) Planetary radar astronomy, in The Encyclopedia of Physical Science and Technology , 
3rd Edition (R. A. Meyers, ed.), Academic Press, 12, 295-328. [2] Ostro S. J. et al. (1992) JGR 97, 18227-18244. 
[3] Black G. J. et al. (2001) Icarus 151, 160-166. [4] Hapke B. (1990) Icarus 88, 407-417. [5] Hapke B. and Blewett 
D. (1991) Nature 352, 46-47. [6] Benner L. A. M. et al. (1997) Icarus 130, 296-312. [7] C. Magri C. et al. (2001) 
MAPS 36, 1697-1709. [8] Ostro S. J. et al. (2002) Asteroid radar astronomy, in Asteroids III (W. Bottke, A. Cellino, 
P. Paolicchi, and R. P. Binzel, eds ), Univ. of Arizona. [9] Margot J. L. et al. (2002) Science 296, 1445-1448. 
[10] Giorgini J. D .et al. (2002) Science 296, 132-136. 
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Introduction: Particles in a planetary regolith can be 
treated as independent scatterers when their size is 
large compared with the wavelength of light. The na- 
ture of this scattering, however, is poorly understood 
as the particle size approaches and becomes smaller 
than the wavelength. In order to understand this inter- 
action better, the reflectances of we 1 1 -sorted particulate 
samples of known composition were measured. This 
reflectance data was then analyzed to determine how 
scattering properties change with particle size. 

Method: The samples used in this study were 13 alu- 
minum oxide abrasive powders, 6 iron oxide pigment 
powders, and 4 precipitated calcium carbonate pow- 
ders. The aluminum oxides range in size from 0. 1 to 
30.09 pm, the iron oxides range in size from 0.15 to 
2.5 pm, and the calcium carbonates range in size from 
0.06 to 3.0 pm. SEM analysis has been done to verify 
not only the particle size of the samples, but the size 
distribution and the particle morphology. In addition, 
x-ray diffraction has been used to identify the iron ox- 
ide mineralogy. The instruments used to obtain the 
reflectance data are the long and short arm goniometric 
photopolarimeters (or goniometers] at the Jet Propul- 
sion Laboratory. Both instruments measure the inten- 
sity of red laser light reflected from a sample in 8 
senses of polarization (4 linear, 4 circular). The long 
arm goniometer is capable of measuring phase angles 
from 0.05 to 5° with a minimum phase angle ('g') in- 
crement of 0.01°. Nelson et al. (2000) [1] has previ- 
ously presented the data for the aluminum oxide sam- 
ples, as well as describing the instrument. The short 
arm goniometer is capable of measuring phase angles 
from 4 to 160° with a minimum phase angle increment 
of 0.5°. For the data presented here, the incidence an- 
gle (Y) is 60°, making the maximum possible phase 
angle 140°. Data taken with this instrument for a num- 
ber of samples has previously been published [2]. 

For each instrument, a run begins by measuring a 
99% reflectance standard for calibration purposes. The 
standard is measured at i - 0°, g = 5° in all 8 polariza- 
tions. The sample is then measured at the incidence 
angle used for data collection and g = 5° in all polari- 
zations. The results are used to calculate normal albedo 
and to calibrate the sample data. The resulting phase 
curves from both instruments cover both very large 
and very small phase angles. These curves are then 
subjected to a series of theoretical fits to attempt to 
determine scattering parameters such as single particle 
scattering function, single scattering albedo, the inten- 
sities of the coherent backscatter and shadow-hiding 
opposition effects, the half-width half-max (HWHM) 
of the coherent backscatter opposition surge, the trans- 
port mean free path, and the scattering cross section of 
the particles. An improved suite of bidirectional re- 


flectance equations from Hapke, 2002 [3] are used to 
obtain the fits, and an example is shown in Figure 1. 
The single particle scattering function 0p(g)’) and the 
single scattering albedo (‘w’) are first estimated by 
assuming p(g) can be approximated by a 2nd order 
Legendre polynomial expansion, and then fitting the 
phase curve with a Hapke equation fit. If a satisfactory 
fit cannot be obtained, additional approximations for 
p(g) are explored, including higher order Legendre 
polynomial expansions and double Henyey-Greenstein 
functions (expanded in Legendre polynomials). The fit 
is allowed to vary from the data where the coherent 
backscatter opposition surge is apparent in the circular 
polarization ratio data. These fits are then refined by 
adding in terms to account for both the shadow-hiding 
and coherent backscatter opposition surges: the coher- 
ent backscatter term contains a variable that measures 
the HWHM of this opposition surge. Previous work 
has noted that this HWHM varies with particle size 
[4,5]. The HWHM is used to then calculate the trans- 
port mean free path, and this used to determine the 
scattering cross section of the particles. The single 
particle scattering function is used to obtain the aver- 
age scattering angle <0>, which can then be related to 
the extinction and absorption cross-sections. These 
two parameters are normally used as inputs for Mie 
scattering calculations, and can be used to compare 
these lab results with parameters commonly used in 
planetary remote sensing. 

Results: Contrary to theoretical expectations, the alu- 
minum oxide powders all exhibit coherent backscatter 
peaks with HWHM's of less than 2° for both very large 
and very small particles. Coherent backscatter opposi- 
tion surges are apparent in the short ami data only in 
the 2.1, 3.2, and 4.0 pm sized aluminum oxides. On the 
other hand, the iron oxide powders all exhibit coherent 
backscatter peaks with HWHM’s of about 20 degrees, 
showing no dependence on particle size. Likewise, the 
calcium carbonate samples also exhibit coherent back- 
scatter peaks widths with no apparent dependence on 
particle size. The single particle phase functions of all 
the powders exhibit well developed backscatter lobes, 
regardless of what approximation is used to represent 
pig) (see Figures 2,3 and 4). Over the range of phase 
angles studied, none of the powders had the strong 
forward scattering lobes predicted by theoretical mod- 
els based on scattering by spheroids. As the particle 
size decreases, there is no indication of the change 
from Mie-type forward scattering to Rayleigh scatter- 
ing predicted by models that assume the particles 
scatter independently. 
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Phase Curve and Hapke Fit: Calcium Carbonate, 3 0 um 



Figure 1: Example of Hapke equation fit to calibrated 
phase curve data. Single panicle scattering function is 
approximated by a 3 order Legendre polynomial ex- 
pansion. 


Single Scattering Phase Functions Aluminum Oxide 



Figure 3: Single particle scattering functions for alu- 
minum oxide samples, using 5 th order double Henyey- 
Greenstein function expanded in Legendre polynomi- 
als. Curves are offset, and the y-axis scale is the same 
for each. 


Single Scattering Phase Functions Iron Oxide 



Figure 2: Single particle scattering functions for iron 
oxide samples, using 2 nd order Legendre polynomial 
expansion. Curves are offset, and the y-axis scale is 
the same for each. 


Single Scattering Phase Functions: Calcium Carbonate 



Figure 4: Single particle scattering functions for cal- 
cium carbonate samples, using 3 rd order Legendre 
polynomial expansion. Curves are offset, and the y- 
axis scale is the same for each. 
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The Early Years (1960’s and 1970’s). It has often 
been said that exploration of the Moon was imple- 
mented backwards. The Apollo national initiative of 
the 1960’s resulted in landing men on the Moon at six 
sites and returning samples to earth-based laboratories 
for study. These unique and valuable samples brought 
a quantum leap in understanding the age and evolution 
of the solar system. Simultaneously, the former Soviet 
Union developed unmanned sample return methods 
and successfully returned lunar samples from three 
additional locations. Although these remarkable ac- 
complishments were envisioned as the dawning of a 
golden era of space exploration, the emphasis in the 
US quickly turned to developing the Shuttle and using 
unmanned missions to explore more distant horizons. 
Management essentially considered lunar exploration 
closed. In spite of continuous efforts from an ever- 
dwindling community of lunar scientists, no phase of 
detailed study of the Moon with advanced sensors was 
implemented. A well -instrumented lunar polar orbiter 
was proposed countless times, but never flown. Lunar 
exploration went into hibernation for decades. 

The Reawakening (1990’s). Lunar exploration 
was re-initiated almost by accidem. [In reality, it was 
sneaked in by never-give-up lunar scientists.] The 
Galileo spacecraft flew by the Moon twice to get extra 
energy for its trip to Jupiter. Both times the teams were 
convinced it was valuable to test and calibrate systems 
using the Moon. This was the first time relatively 
modem CCD sensors were used in lunar orbit and 
Galileo provided the first new orbital spectroscopic 
data for the Moon [1]. A few years later DoD decided 
to test a group of miniaturized cameras designed to 
measure fast moving targets. For political reasons they 
needed to go outside near-Earth space and approached 
NASA to ask if there was scientific interest in the data 
to be produced for the Moon and an asteroid. The lunar 
science diehards were elated, and soon a Clementine 
science team was formed to plan and implement sci- 
ence data [2]. After countless calibration issues (some 
still ongoing) the first global topography [3] and global 
multispectral imagery data sets were produced for the 
Moon [e.g. 4], The only planned mission to the Moon 
in this second wave of exploration was Lunar Pros- 
pector, the cheapest Discovery mission ever selected 
by NASA. This little mission provided the first global 
assessment of elemental abundance (especially Th, K, 
Fe, Ti, H) [5]. 

The Surprises. The combined new pulse of lunar 
data knocked the socks off complacent planetary sci- 


entists who thought they understood the Moon. This is 
not to say important research had not continued via 
lunar sample studies and telescopic measurements of 
the lunar nearside during the dark decades. The origin 
of the Moon and the Earth/Moon system through a 
violent planetismal encounter was well on the way to 
being understood [6] But our perception and under- 
standing of the Moon changed forever with this little 
bit of new data from relatively modem sensors. Three 
particularly important surprises are highlighted below. 

SPA, the Biggest Basin in the Solar System . Al- 
though several of the diehard lunar enthusiasts had 
long argued for a big basin on the lunar farside [e.g. 7], 
the South Pole-Aitken basin (SPA) is missing from 
most commercial lunar maps produced using Apollo 
era data. This remarkable oversight reflects the limited 
character of the data available. The first digital images 
from Galileo changed everything, however. Not only 
was the basin rim readily recognized, but the basin 
interior had a distinctively low albedo and exhibited 
spectral properties indicating that it was rich in mafic 
minerals [1]. Subsequent data from Clementine illus- 
trated the dramatic topography (~12 km deep!) [3] and 
unusual mineralogy exposed in the interior [8]. It be- 
came apparent that SPA may be the largest docu- 
mented basin in the solar system [3, 9], and the deep- 
seated material of the interior represents a terraine on 
the Moon that is quite distinctive and unlike the 
lithologies documented with samples returned from the 
nearside [10]. 

Thorium Hot Spot. The global elemental data from 
Lunar Prospector confirmed that the SPA interior is 
iron-rich (midway between feldspathic highlands and 
mafic mare basalt). One of the most unexpected re- 
sults, however, was the very well-defined pattern of 
radiogenic materials, almost all of which are associated 
with the nearside Imbrium basin and its deposits. Since 
such elements (KREEP material) are expected to be 
concentrated in the residual of primordial lunar differ- 
entiation, it is not surprising that a large basin tapped 
such a concentration. The surprise is that the pattern is 
only associated with the Imbrium basin, which is not 
only smaller than SPA but also occurred near the end 
of the period of basin forming events. This highly lo- 
calized concentration of radiogenic (incompatible) 
elements indicates there are enormous gaps in our un- 
derstanding of the first few hundred million years of 
planetary evolution and differentiation. 

Polar H (water?). The extremely low volatile con- 
tent of lunar samples was one of the first characteris- 
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tics noted for returned lunar samples (along with their 
ancient age). Thus, the hint of polar ice by Clementine 
[11] and the unquestioned discovery of hydrogen con* 
centration at the lunar poles by Prospector [12] has led 
to great optimism that water resources may indeed ex- 
ist in the permanently shadowed regions of the Moon. 
This counter-intuitive concept, although reasonable in 
theory [13], is made more credible by the discovery of 
extensive volatile (perhaps ice) deposits in craters at 
Mercury’s poles [14]. 

Meanwhile, back at the ranch.... One of the more 
perplexing puzzles presented by the returned lunar 
samples was that lunar soils don’t really look like lunar 
rocks at visible to midIR wavelengths. This is a key 
concept for remote sensing studies, since data obtained 
by remote sensing provides the only information on 
local, regional, and global properties of the Moon and 
soils invariably dominate remotely acquired data. Un- 
derstanding soil processes is thus critical and it has 
taken close to 30 years to resolve. Reinvestigating lu- 
nar soils with modem instruments in earth-based labo- 
ratories identified the importance of nano-scale proc- 
esses that occur on soil grains when exposed to the 
space environment [15, 16]. Although the seeds were 
planted from the beginning [17], it is the cumulative 
insight from detailed laboratory and theory investiga- 
tions that allowed identification of the dominant role of 
nanophase reduced iron in altering optical properties 
[18] and provided the modem paradigm of regolith 
processes active in the space environment. 

The Critical Role of Samples. What can we leam 
from our experience of using samples in exploration of 
the Moon? Some aspects seem deja vu. Mars dust and 
soil doesn’t look anything like crustal rocks (SNC 
meteorites); nor do typical meteorites look like the 
surfaces of most asteroids, that might represent likely 
parent bodies. Nevertheless, it appears to be assumed 
that remote sensing and in-si tu measurements of other 
bodies will provide sufficient information on the dis- 
tribution and type of materials to provide key geologic 
context needed to unravel planetary evolution. The 
lunar experience provides a few cautionary insights. 

The first is that without samples of lunar soil in our 
laboratories we would not have been able to under- 
stand (read interpret) the available remote sensing data 
sufficiently to discern the primary character of the 
Moon. Try to imagine our understanding without the 
samples. We would have detected the presence of py- 
roxene and would have noted that it appears to be more 
abundant at fresh craters for some reason. The distri- 
bution of iron from Lunar Prospector and the Apollo 
orbital experiments might have led to the supposition 
that all the maria are basaltic, but there would be no 
conclusive evidence. We would puzzle about the 


weakness of pyroxene features in soil, perhaps hy- 
pothesizing it is due to a high degree of mysterious 
contamination. The highlands are older than the maria 
(based on superposition and crater density), but we 
would have no reference for actual age. We would 
know the highlands are low in iron, but might think the 
higher albedo is due to smaller particles resulting from 
more impacts. The weaker absorptions in the highlands 
might also be due to the greater age, and hence con- 
tamination. Such musings would all be nonsense, a 
gross misconception based on ignorance of funda- 
mental information about lunar soils and their relation 
to primary rocks. 

A second related example is that without samples 
of key lunar lithologies (high-Ti and low-Ti basalt, 
feldspathic breccias, KREEP, Mg-suite), we cannot 
discern the inter-relations between different rock types, 
so essential to formulating the basic geologic history 
(magma ocean with mafic mantle cumulates and 
crustal anorthositic cumulates). In addition to petro- 
logic relations, the precision needed for diagnostic mi- 
nor and trace elements and isotopic abundances re- 
quires Earth-based laboratories. The importance of the 
Thorium Hot Spot has special meaning to lunar scien- 
tists because it most likely represents concentration of 
the magma ocean residuum in a very localized region 
of the Moon. The unique mafic-rich intenor of SPA 
presents a challenge to what we think we know about 
the structure and composition of the lunar lower crust 
and mantle. If we knew in the 1960’s what we know 
now, we would certainly have argued long and hard 
with the engineers to obtain a representative sample of 
SPA. As a footnote, it is also recognized that, at least 
for the Moon, soils provide a remarkably representa- 
tive suite of materials of local and foreign origin [19]. 
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Abstract: We describe the fust intercontinental 
planetary radar test performed in Italy observing the 
near Earth asteroid (NEA) 33342 (1998 WT24) in De- 
cember 2001 by means of the bist itic configurations 
Goldstone (California, USA)-Medidna (Italy) and 
Evpatoria (Ukraine)-Medicina. 

The experiment goal was to characterise the system for 
realtime radar follow-up observations of NEAs and 
artificial orbiting debris, in the frame work of a feasibil- 
ity study which aims at using the Sardinia Radio Tele- 
scope, at present under construction, also as a planetary 
radar facility. 

We report the preliminary results of the radar observa- 
tions carried out by the IRA-CNR (Istituto di Radio- 
astronomia - Consiglio Nazionale cielle Ricerche) and 
the OATo (Osservatorio Astronomico di Torino) 
groups, aimed at exploring the scientific potentials of a 
new space radar program, using the existing facilities 
in Italy. 

The planetary radar technique is uniquely capable of 
investigating geometry and surface properties of vari- 
ous solar system objects [1], demon strating advantages 
over the optical methods in its high spatial resolution 
and ability to obtain three-dimens ional images. 

A single radar detection allows to obtain extremely 
accurate orbital elements, improving the instantaneous 
positional uncertainties by orders of magnitude with 
respect to an optically determined orbit. Radar is a 
powerful means to spatially resolve NEAs by measur- 
ing the distribution of the echo power in time delay 
(range) and Doppler frequency (line-of-sight velocity) 
with extreme precision in each coordinate, as it pro- 
vides detailed information about the target physical 
properties like size, shape, rotation, near-surface bulk 
density and roughness and internal density distribu- 
tion 1]. 

The Medicina 32m antenna had been successfully used 
for the first time as the receiving part of a bistatic con- 
figuration during a test experiment (September 2001) 
held to check the capabilities of the entire data acquisi- 
tion system. This test was possible thanks to the col- 
laboration undertaken with the Evpatoria radar station, 
and consisted in the observation of the ETALON-1 low 
orbit satellite. 



Realtime echo received in Medicina (Italy) 

References: [1] S.J. Ostro (1993) “Planetary Radar 
Astronomy” Rev. Modern Phys vol.65, no .4 
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Introduction: The peak of the solar energy distribu- 
tion occurs at visual wavelengths and falls off rapidly 
in the infrared. This fact, improvements in infrared 
detector technology, and the low surface temperatures 
for most icy objects in the outer solar system have re- 
sulted in the bulk of telescopic and spacecraft observa- 
tions being performed at visual and near-infrared 
wavelengths. Such observations, begun in the early 
1970’s and continuing to present, have provided com- 
positional information regarding the surfaces of the 
satellites of Saturn and Uranus, Neptune’s moon Tri- 
ton, Pluto, Pluto’s moon Charon, Centaur objects, and 
Kuiper belt objects. Because the incident sunlight 
penetrates the surface and interacts with the materials 
present there, the measured reflected sunlight contains 
information regarding the surface materials, and the 
ratio of the reflected to incident sunlight provides a 
mechanism of identifying the materials that are pre- 
sent. 

Fink and Sill [1] and Roush [2] present the reflectance 
spectra of various ices, An initial identification of the 
species present on solar system objects can be provided by 
comparison with these laboratory spectra. More detailed 
information, such as relative abundances grain sizes of 
surface constituents, requires more detailed models. 



Figure 1 . The cartoon on the right fjf illustrates spatial 
mixing of spectra where the reflect am v is dominated by 
photons encountering individual materials. Along some 
pathways a component of intimate mix ing remains . : The 
upper- left cartoon illustrates intimate spectral mixing 
where the reflectance is dominated by photons that 
encounter several materials. The lower-left cartoon 
illustrates Intimate spectral mixing with a surface layer. 

Various mathematical models are used to character- 
ize the chemical and physical nature of the surfaces 
under investigation. They can also be used to extend 
the results of laboratory studies to different mixtures 
and particle sizes that are not directly investigated. 


Spatial mixtures (see Fig. 1) describe photons re- 
flected from the surface that are dominated by interac- 
tions with individual surface species. The total sur- 
face reflectance is given by F r = X/XR/ + XjxRj + ••• 
+ XkxR m where X is the areal extent, and R the reflec- 
tance of each component, respectively, and it is re- 
quired that LA/ = 1. Examples of these types of sur- 
faces are individual rocks or outcrops. 

In intimate, or "granular", mixtures (see Fig. 1) the 
surface reflectance is a non-linear function of the re- 
flectances of each component because a photon en- 
counters several grains, due to multiple scattering, 
before exiting the surface. Examples are individual 
minerals that are contained in rocks and planetary 
soils and regoliths. 

A series of papers summarized in Hap/ce [4] devel- 
oped an approximation to the intimate mixture case 
allowi the reflectance to be calculated from the optical 
properties of the surface components. This important 
contribution allowed the quantitative estimation of the 
relative abundances and grain sizes of surface compo- 
nents and, given the appropriate observations, some 
information regarding the surface physical structure. 
Table 1 summarizes some of these efforts applied to 
objects in the outer solar system. 


Table I. Outer Solar System Surfaces 

Modeled Using Hapke Theory 

Object Property Reference 


Europa 

T 

5 

Ganymede 

T 

5 

Callisto 

T, C 

5.16-8] 

Enceladus 

T 

5 

Tethys 

T 

5 

Dione 

T, C 

5,9 

Rhea 

T, C 

5,9 

lapetus 

T, C 

5, [10-11] 

Rings 

T 

5 

Miranda 

C 

12 

Ariel 

T, C 

5, 13 

Umbriel 

T 

5 

Titania 

T, C 

5, [13-14] 

Oberon 

T,C 

5,(13-14] 

Triton 

T, C 

[15-16], [17-19] 

Pluto 

T, C 

20, [21-25] 

Charon 

Trans- 

T, C 

26, [23,26-30] 

Neptune 

objects 

C 

[31-33] 


T= Temperature, C= composition 
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Figure 2, after Cruikshank et a 1 [32], is an illustra- 
tion of the application of Hapke modeling to the inter- 
pretation of the spectrum of the Centaur Pholus. In 
the mode! the surface is represented by a spatial mix- 
ture with one area being carbon black and the other an 
intimate mixture of components. 



Figure 2. Comparison of a Hapke model of the sur- 
face wiih the measured spectrum of Pholus. At the top 
is a calculated spectrum of each pure material 

Molecular mixtures occur at atomic scales where 
one atom or molecule is substituted into the crystalline 
structure of a host material, for example, CH, in N^. 
On the basis of subtle variations in the positions of 
methane bands, Owen et al [21] and Cruikshank et al 
[23] suggested that on Pluto the methane ice exists in 
two reservoirs: one relatively pure methane and the 
other mixed in nitrogen. Doute et al [24] modeled 
Pluto’s surface as a complex mixture including verti- 
cal layering (see Fig. 2). 

The long list of Table 1 may erronously lead one to 
conclude that compositional interpretation of icy sur- 
faces in the outer solar system is readily available via 
application of Hapke's [4] theory. Most ices are non- 
absorbing at visual wavelengths and the visual colors 
of the various icy objects located in the outer solar 
system clearly demonstrate that non-ice materials are 
present. Some optical constants are available for lim- 
ited materials exhibiting a range of colors at visual 
wavelengths [e.g. 34-35]. However, the majority of 


these have been determined with the material at ambi- 
ent terrestrial conditions. How these characteristics 
change at conditions more appropriate to the outer 
solar system is largly unknown, although some infor- 
mation is becoming available [36]. Additionaly, opti- 
cal constants for a variety of other candidate materials, 
e.g. silicates, are sorely lacking at visual and near-ir 
wavelengths. It is clear that there is still ample labo- 
ratory work to undertake. 

From an observational perspective, the ability to 
achieve more unique interpretative solutions would 
benefit from spectral observations over a wide range of 
viewing geometries. Such observations are essentially 
impossible from ground-based observatories and await 
spacecraft measurements. 

I gratefully acknowledge research support from 
NASA's Planetary Geology & Geophysics Program 
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Introduction: In 1975, while examining why lunar 
regolith becomes darker and redder with time, Bruce 
Hapke with colleagues advocated that the lunar soil 
grains should have been coated with vapor-deposited 
rim containing nanometer-sized iron particles [11. 
High-velocity dust impacts and/or solar wind irradia- 
tions might produce ferrous silicate vapor. This process 
is called “space weathering”. However, this hypothesis 
had been ignored until nanophase iron particles were 
confirmed in lunar regolith grains b> TEM [2,3]. 

S-type asteroids, majority in asteroids, are believed 
to be parent bodies of ordinary chondrites, which are a 
large majority in meteorites. However, S-type steroids 
exhibit more overall depletion and reddening of spectra, 
and more weakening of absorption bands relative to 
ordinary chondrites. This spectral mismatch is also 
considered being caused by space weathering process. 
Recent observation by Galileo showed that fresh crater 
interior and ejecta on S-type Ida have the spectra closer 
to that of ordinary chondrites [4]. Some small asteroids 
have intermediate spectra between S-type and Q-type 
(ordinary-chondrite like) spectra [5]. NEAR’s X-ray 
spectrometer showed that S-type 433 Eros has elemental 
composition close to ordinary chondrites [6]. 

Pulse Laser Irradiation Simulating Dust Im- 
pacts: Although van de Graaf type electrostatic dust 
accelerators can accelerate submicron dust particles to 
be more than 10 km/s, they can hardly irradiate enough 
number of dust impacts to simulate space weathering. 
Moreover such accelerators use conductive dust mate- 
rials such as Fe, Ag, and C, which would affect the 
optical property of targets. Therefore, to simulate the 
high-velocity dust impact heating, we irradiate powder 
and pellet samples of silicates by a pulse laser beam 
(1064 nm) with pulse duration 6-8 ns, which is compa- 
rable with a real dust impact [7,8]. The total irradiation 
energy in unit area is 240 mJ mm 1 for 30 mJ pulse 
energy. Energy deposition rate is about 10 1(> WcnV 2 , 
where evaporation and ion formation is strongly ex- 
pected [9]. 

Laser- irradiated samples show significant redden- 
ing: reduction of spectra is much larger in a visible 
region than in a near-infrared region [7,8]. We use 
olivine from San Carlos with 8.97 vvt% FeO, enstatite 
from Bamble, Norway, with 9.88 wt% FeO and hyper- 
sthene from Mantyharju, Finland, with 16.7 wt% FeO. 


Reflectances of olivines are more easily changed than 
those of pyroxenes [8] This is compatible with the 
compiled asteroid data: using areas of 1 jum and 2 
absorption bands, olivine-rich asteroids have more 
reddened spectra than pyroxene-rich asteroids [10,11]. 

In Fig. 1, spectra are compared with observed 
spectra of olivine-rich asteroids. 446 Aetemitas spec- 
trum is close to the olivine spectrum after 30mJ irra- 
diation. The scaled spectra of irradiated olivine look 
like compatible with rather flat asteroid spectrum at 
wavelengths longward of 1.5 /xm [12]. A more reddened 
slope of 863 Benkoela spectrum is reproduced by that of 
five times 30 mJ irradiations at wavelengths shorter than 
1.6 jxm although repeated irradiation produced con- 
tinuous reddened slope at longer wavelengths. 

We also obtained following results: (1) After twice 
irradiation of 30 mJ pulse laser scanning, 510 nm band 
of enstatite was weakened [13]. This corresponds to the 
change of iron oxidation state. This effect is also dis- 
cussed on lunar sample data. (2) We irradiate the 
mixtures of olivine and pyroxene with various mixing 
ratio by the pulse laser. Although the irradiated spectra 
of Ol-Px mixtures depend on the mixing ratio, the al- 
teration is not proportional to the ratio. 
Space-weathered spectra of Ol-Px mixtures can be pre- 
dicted using the end-member spectra of them [14]. 
(3) We have confirmed spectral changes (darkening and 
reddening) on an ordinary chondrite sample and a ser- 
pentine sample. 

Nanophase Iron Particles: To clarify microscopic 
cause of reflectance change, we observed the irradiated 
samples by high-resolution scanning electron micros- 
copy (SEM) and transmission electron microscopy 
(TEM). Thin sections (about 80-100 nm thickness) of 
irradiated grains were observed by 200keV 
High-Resolution TEM equipped with EDS. 

In the SEM image of irradiated olivine, there are 
observed a number of small holes which were formed 
through surface evaporation (Fig. 2). In the TEM image 
of a thin section of an olivine grain, a large number of 
nanophase particles (several to 20 nm) are observed in 
the amorphous rim (Fig. 3(a)). These are crystalline iron 
particles [15,16]. These particles are responsible for 
optical properties changes such as darkening and red- 
dening as proposed [1] and calculated [17]. These 
particles may be related to nanophase iron in lunar soils 
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[2,3]. In the TEM image from a laser-irradiated ensta- 
tite sample, nanophase iron particles are observed in 
amorphous aggregates (Fig. 3(b)). Our results show that 
hydrogen is not necessary to form nanophase iron. 

In order to quantify the amount of nanophase iron, 
ESR (electron spin resonance) is measured for irradiated 
pellet samples by a microwave cavity method where 
spins in the surface depth of 100 micron can be detected 
nondestructive^ [18]. Strong ESR signals, which derive 
from nanophase iron particles, are observed both in 
irradiated olivine and pyroxene samples (Fig. 4). 
Moreover, ESR intensities increase with the space 
weathering degree simulated as laser irradiation time. 

References: [1] Hapke B. et al., (1975) Moon 13, 
339-353, [2] Keller L. P. and McKay D. S. (1993) 
Science 261, 1305-1307. [3] Keller L. P. and McKay D. 
S. (1997) GCA. 61, 2331-2341, [4] Chapman C. R. 
(1996) Meteoritics & Planet . ScL, 31, 699-725. 

[5] Binzel R, et al. (1996) Science 273, 946-848. 

[6] Trombka J. I. et al. (2000) Science 289, 2101-2105. 

[7] Yamada M. et al. (1998) Antarctic Meteorites 
XXIII, 173-176. [8] Yamada M. et al. (1999) Earth 
Planets Space 51, 1255-1265. [9] Kissel J. and Krueger 
F. R. (1987) Appl. Phys . A 42, 69-85. [10] Hiroi T. and 
Sasaki S. (1999) LPS XXX, #1444. [11] Hiroi T. and 
Sasaki S. (2001) Meteoritics &. Planet. Sci. 36, 
1587-1596. [12] Clark B. E. and Hiroi T. (1994) Bull . 
Am. Astron. Soc. 26, 1172. [13] Hiroi T. et al. (2001) 
Earth Planets Space 53, 1071-1075. [14] Kurahashi E. 
and Sasaki S. (2002) LPS XXXIII, #1479. [15] Sasaki, 
S. et al. (2001) Nature 410, 555-557. [16] Nakamura K. 
et al. (2001) LPS XXXII, #1547 [17] Hapke, B. (2001) 
J. Geophys. Res. 106, 10039-10073. [18] Ikeya M. etal. 
(1994) Rev. Sci. Instrum . 65, 3670-3672. 
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Figure 1. Reflectance spectra of olivine pellet samples 
before and after pulse laser irradiation. Olivine spectra 
scaled at 550nm. The scaled spectra are compared with 
asteroids 446 Aetemitas and 863 Benkoela. 



Figure 2. Backscattered SEM image of a pulse-laser 
irradiated olivine pellet surface (30mJ x 10). Width of 
the figure is 40micron. 



Figure 3 (a) A bright-field TEM image of a thin section 
from a laser irradiated olivine sample (30mJ x 20). 
Nanophase iron particles whose size is several to 20 nm 
are within the olivine rim. (b) A bright-field TEM image 
of a thin section from a pulse-laser irradiated enstatite 
pellet sample (30mJ x 20). Nanophase iron particles are 
embedded in amorphous aggregates. 
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Figure 4. ESR spectra of olivine samples, which were 
irradiated by 0, 1,5, 10 times at 30mJ. 
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Introduction: The exploration of the surfaces of 
solar system bodies (planets, satellite, asteroids, ...) now 
systematically involves observations by imagers and 
imaging spectrometers over the UV, visible and near- 
infrared ranges. All along the missions the light reflected 
by the surface is recorded under various 
illumination/observation geometries, often constrained 
by the spacecraft orbit. Large changes of reflectance and 
spectral contrast with phase angle are generally observed 
for a given area [1], Also, contrary to terrestrial snow or 
ice fields, a tendency for backscattering behavior is 
widely found for the surfaces of the icy satellites [2]. 

In particular, lo’s Jupiter satellite displays very 
strong contrast reversals between surface units at low 
and high phase angles [3]. This behavior seems to be 
linked to very different BRDF for the various materials 
constituting its volcanic surface dominated by sulfur 
compounds and S0 2 frost. Several “a priori” interpreta- 
tions have been proposed for this peculiar behavior but 
none has been theoretically modeled or tested experi- 
mentally. 

The correct interpretation of the spectra produced by 
the imaging spectrometers under varied geometries first 
requires to be able to “correct” for these bidirectional 
effects. If sufficiently sampled, the BRDF of a planetary 
surface also potentially contain interesting textural 
information on its constituting materials. However, first 
of all laboratory investigations are clearly needed in 
order to understand the physical and textural parameters 
that control the BRDF. For that purpose a spectrophoto- 
goniometer has been designed at LPG to measure the 
spectral bidirectional reflectance distribution function 
(SBRDF) of granular and compact geophysical surfaces. 
The instrument covers the solar spectrum, from 0.3 to 5 
pm, and enables a wide range of bidirectional 
configurations [4]. 

BRDF Measurements: We will first present the 
instrument and its performances. Then we will lay out 
the BRDF measurements performed over the visible 
spectral range on sulfur samples with various textures 
obtained by different physical processes that may occur 
on planetary surfaces: sublimation, thermal sintering, gas 
condensation, melting/solidification (5, 6\. 


Several sulfur powders exhibit a noticeable backscat- 
tering behavior, similar to some of the bright areas on Io. 
Such behavior, not predicted by the current single 
scattering models of large particles [7], seems to be 
linked with the sub-micron texture of the particles rather 
than with their average size and shape. During thermal 
sintering of these powders a more isotropic behavior 
develops. A very strong specular reflection is observed 
for melted/solidified sulfur but also, to a lower extend, 
for crystalline powders obtained by slow gas condensa- 
tion at low pressure. All these directional behaviors are 
modulated by the strongly wavelength dependent absorp- 
tion of sulfur. 

Conclusion: These experiments and their compari- 
son show that the BRDF signature, coupled with the 
spectral information, should help to discern between 
different textural states of a given surface material, and 
thus can give some clues on its formation and evolution 
processes. However adequate bidirectional sampling of 
the observations is requisite to correctly retrieve this 
information. 

References: [1] Doute, S., Schmitt B., Lopez - 
Gautier R., Carlson R., Soderblom L., Shirley J. and the 
Galileo NIMS Team (2001) Icarus, 149, 107-132. [2] 
Verbiscer A., and Helfenstein P. ( 1 998) In Solar System 
Ices, (B. Schmitt, et al. eds.) Kluwer Acad. Publ., 
Dordrecht, Astrophys. Space Sci. Lib., To/. 22? \ 157- 
198. [3] Geissler, P.E.. , McEwen A., Phillips C., 
Simonelli D., Lopes-Gautier R. and Doute S. (2001) J. 
Geophys. Res £ \ in press. [4] Brissaud O., Schmitt B., 
Bonnefoy N., Doute S., Fily M., Grundy W., and 
Rabou P. (2002) App/. Optics, in preparation. [5] 
Bonnefoy N., Brissaud O., Schmitt B., Doute S.,FilyM., 
Grundy W., and Rabou P. (2000) Remote Sensing. Rev., 
19, 59-74. [6] Bonnefoy N. (2001) Thesis , UJF-LPG 
Grenoble. [7] Grundy, W.M., Doute S. and Schmitt B. 
(2000) J. Geophys. Res £ f 105, 29,291-29,314. 
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Figure l: Bidirectional Reflectance Distribution- 

Function (BRDF) in the principal illuminatin plane 
(source incidence: 60°) at A = 492 nm for 4 sulfur 
samples: a) (top left) powders with “nominal” grain size 
of~40 and 160 pm, b) (top right) thermaly sintered 40 
pm powder ( 1 2 days at 90°C), c) (bottom left) gas phase 
condensed crystals under low pressure (2 10 2 mbar), d) 
(bottom right) melted (125 °C) and reffozen sample. 
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Introduction: We report on a multifaceted inves- 
tigation of remote sensing data for the Terra Meridiani 
region. Our study area extends ±10 degrees about the 
prime meridian and from -5 to +10 degrees in latitude, 
encompassing the largest of the hematite bearing re- 
gions identified on the surface of Mars [1]. A detailed 
investigation of this area was motivated by the likely 
selection of the hematite-bearing unit as the destination 
for one of the twin MER rovers. This study makes use 
of a variety of remote sensing data sets with the goal of 
characterizing the region to the fullest possible extent. 
The initial results from this body of work have con- 
tributed to the MER landing site selection process, and 
will aid in mission planning and operations if one of 
the rovers is sent to this region. 

Remote Sensing Units Mapping: Through extern 
sive examination and interpretation of MOC narrow 
angle and wide angle images [2], assisted by analysis 
of TES-derived bolometric albedo [3], TES-derived 
thermal inertia [4], MOLA-derived pulse width data (a 
surface roughness measure) [5], and the TES-based 
hematite index map [1], our study area was partitioned 
into five primary surface units (Figure 1). 

The MOC image analysis also established that 
these units are set in a clear stratigraphic sequence [6]. 
The dissected cratered terrain exposed across the 
southern portion of the study area is mapped as the 
typical Noachian cratered terrain that pervades the 
southern hemisphere. The etched unit overlies the dis- 
sected terrain and is interpreted to be a severely de- 
graded volcanic complex based on morphologic com- 
parisons to known volcanic regions of Mars such as 
Athabasca. Notably, the etched unit has both an inter- 
mediate albedo and a very high thermal inertia, a com- 
bination that runs counter to the global norm. The 
hematite unit and plains unit overlie the etched unit 
and exhibit identical morphologic character. However, 
the plains unit lacks the elevated hematite index values 
that distinguish the hematite unit. The mantled cratered 
terrain extends to the north and is characterized by a 
dust mantling several meters thick. Our current model 
of the geologic history of this region terminates with 
the stripping back of the dust mantle still present to the 
north from the remainder of the study area, exhuming 
the plains, hematite, and etched units. 

TES Hyperspectral Data: To facilitate the spec- 
tral analysis of the study region, a data cube with 256 
pixels/degree and 138 spectral channels was con- 
structed. Each TES spectrum was localized in the data 


cube according to the coordinates of the observation 
center, with each spectrum being relegated to a single 
spatial pixel. While this is not an accurate reproduction 
of the TES footprint on the surface, this configuration 
does allow for both spatially orientated analyses and 
the interpretation of individual spectra. Further, if the 
intent is to pass the spectral data to a deconvolution or 
analysis routine (as is the case here) this method of 
data localization guarantees that each spectrum is 
processed only once, instead of multiple times as 
would be the case if the spectrum were replicated in all 
of the spatial pixels covered by the actual footprint. 

Spectral deconvolution/unmixing. The TES data 
cube was processed using a variant of the linear de- 
convolution algorithm advocated by the TES team [7]. 
This technique is capable of simultaneously decon- 
volving the atmospheric and surface signals and un- 
mixing the derived surface spectrum in terms of a set 
of mineralogic endmembers. The dual deconvolu- 
tion/unmixing procedure was performed with an end- 
member set consisting of six atmospheric endmember 
spectra provided by Bandfield, one blackbody spec- 
trum to account for variations in spectral contrast, and 
-40 pure mineral spectra that collectively cover a wide 
range of possible surface mineralogies. 

The deconvolution algorithm generates a series of 
resultant data cubes, all of which have the same spatial 
dimensions as the input. Resultant data cubes include: 
1) the best-fit parameters for each spectral endmember 
(atmospheric and surface); 2) the parameter value as- 
signed to the blackbody endmember and the RMS er- 
ror of the spectral fit; 3) the best fit model spectrum 
over the wavenumber interval for which the linearity 
of the convolution holds; 4) the derived atmospheric 
spectrum and 5) the derived surface spectrum over the 
same wavenumber interval; 6) the normalized abun- 
dance of the atmospheric endmembers; and 7) the 
normalized abundance of the surface endmembers. 

Data interpolation and map generation. Prior to 
the spatial interpolation of the deconvolution results, 
the RMS error band is used to eliminate low quality 
surface observations from consideration. While the 
data used to populate the original data cube was fil- 
tered for quality (e.g. day side, nadir observations, 
without flagged errors), the deconvolution error for 
some orbits is anomalously large. This can be attrib- 
uted to high atmospheric opacity at the time of obser- 
vation as the linearity of the spectral convolution 
breaks down for high opacity observations [7]. High 
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RMS errors are also common for dusty regions, an 
occurrence that can be explained by the nonlinear par- 
ticle size dependency of emission spectra [8]. How- 
ever, reducing the data density for dusty surface re- 
gions is an acceptable trade-off for eliminating high- 
opacity orbits from consideration. The actual interpola- 
tion scheme that generates well-populated endmember 
maps also takes into account the RMS error of the de- 
convolution, and gives observations with low RMS 
errors greater weight. As a result, the statistical analy- 
sis of the interpolated maps will also be effectively 
error-weighted, with low-error observations attaining 
greater significance. Interpolated endmember maps can 
be appropriately combined to produce component 
maps of, for example, hematite abundance (Figure 2), 
or olivine abundance (Figure 3), for the entire region. 

Results: The deconvolution/unmixing results and 
interpolated component maps provide valuable insight 
into the mineralogical variation among the surface 
units. Specific results include: 1) The mapping of the 
surface mineralogy served to bolster the validity of the 
unit boundaries, as all of the units are evident in the 
interpolated component maps. 2) Spectrally (and min- 
eralogically) the etched unit differs significantly from 
the dissected unit only in clinopyroxene and glass 
abundance. This is a surprising result, given the differ- 
ences these units exhibit in other characterizing quanti- 
ties (e.g. albedo and thermal inertia). 3) The most 
striking correlation among the component maps is the 
association of iron-rich (~Fo 3 i) olivine with hematite 
(Figures 2-3). These components exhibit a positive 
correlation within all of the individual surface units, as 
well as across the study area as a whole. Table 1 gives 
the correlation coefficients for these components over 
various spatial regions. 


Normalized Hematite Component 



Figure 2. Normalized hematite component for the Terra 
Meridiani study region. 


References: [1] Christensen P. R. et al. (2000) 
JGR, 105 , 9623-9642. [2] Caplinger M. A. (2002) 
LPSC XXX111 , 1405. [3] Christensen P. R. Personal 
communication. [4] Jakosky B. M. et al. (2000) JGR , 
105, 9643-9652. [5] Garvin J. B. et al. (1999) GRL, 
26, 381-384. [6] Arvidson R. E. et al. (2002) in 
preparation . [7] Smith M. D. and Bandfield J. L. 
(2000) JGR , 105, 9589-9607. [8] Hamilton V. E. 
(1999) LPSC XXX, 2001. 


Table 1. Correlation coefficient of olivine and 
hematite components for selected regions. 


Region 

Correlation Coefficient 

Full scene 

0.624 

Full scene excepting MCT 

0.481 

Ph only 

0.317 
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Figure 1. Units map for the Terra Meridiani study re- 
gion. MCT: mantled cratered terrain; P: plains unit; Ph: 
hematite-bearing plains unit; E: etched unit; DCT: dis- 
sected cratered terrain. 
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Figure 3. Normalized olivine component for the Terra 
Meridiani study region. 
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Introduction: The Bloomsburg University Goniometer 
Laboratory (B.U.G. Lab) consists of three systems for 
studying the photometric properties of samples The primary 
system is an automated goniometer capable of measuring the 
entire bi-directional reflectance distribution function (BRDF) 
of samples. Secondary systems include a reflectance 
spectrometer and digital video camera with macro zoom lens 
for characterizing and documenting other physical properties 
of measured samples. Works completed or in progress 
include the characterization of the BRDF of calibration 
surfaces for the 2003 Mars Exploration Rovers (MER03), 
Martian analog soils including JSC-Mars-l, and tests of 
photometric models. 

B.U.G. Lab: The B.U. goniometer is an automated 3- 
axis system, providing for motion in incidence (-60° to +60°, 
emission (-80° to +80°), and azimuth angle (0°-180°), thus 
allowing an entire half-hemisphere of scattering geometries 
to be measured (Fig. 1). Assuming the sample has a random 
(or at least azimuthal ly uniform) surface texture, the other 
half-hemisphere is the same by symmetry across the 
scattering plane. Maximum phase angle is 140°, minimum 
phase angle is 3°. Computer controlled stepper motors move 
all three axes with a precision of ±0.00 V . The light source is 
a 100W quartz-tungsten halogen bulb and is powered by a 
radiometric-grade power supply. The output light is 
collimated, filtered with 50nm FWHM interference filters, 
chopped, and fed to the incidence ami along a 1.3m fiber 
optic bundle. Collimating optics at the top of the incidence 
arm provide a 2cm spot on the sample from a distance of 
56cm. The detector, sensitive to light from 400nm to 
lOOOnm, consists of a silicon chip 1cm 2 , amplifiers, and 
lock-in electronics that detect only the light chopped at the 
frequency and phase determined by the chopper on the 
source. It is located at the top of the emission arm a distance 
of 77cm from the sample. The detector’s field-of-view is 
restricted by a light shield giving an un-vignetted FOV of 
2.8° and vignetted FOV of 25°. The sample dish is 57mm in 
diameter and 10mm deep, sits at the intersection of the three 
axes, and is stationary. Secondary reflections from the 
equipment surrounding the sample are minimized by 
flocking everything within the FOV of the detector. Control 
of the stepper motors and detector readout is completely 
automated and achieved using custom software written with 
LabView™. Axes motion sequences art preprogrammed and 
saved. A typical run at a single wavelength measures the 
scattering behavior from 0° to 60° incidence (/) in 15° steps, 
emission (e) from 0° to 80° in 10° slips, and azimuth (a) 
from 0° to 180° in variable steps to maintain a consistent 
solid angle between measurements of -0.05 steradian. In all. 


each run of this type returns 680 data points and takes 82 
minutes. Radiance factors are determined by dividing all 
sample radiances by that of a spectrahalon disk at ^=0°, z=5°, 
and multiplying the result by cos(/)/cos(^). The RMS 
precision of the 680 radiance factors are typically <1%, while 
absolute accuracy is estimated to be <10%. 

Spectrometer and Digital yideo Camera. A single 
grating spectrometer documents the reflectance spectrum 
(relative to spectrahalon) of samples from 400- lOOOnm with 
a spectral resolution of ~6nm A Panasonic color video 
camera fitted with a macro-zoom lens is used to 
photographically document the samples. Typically two 
images are taken - one showing the entire sample dish (FOV 
~6cm, resolution of -100pm per pixel), the second zoomed 
in (FOV of l cm, resolution of- 1 5pm per pixel). 

Examples: Figures 2a-c show an example of the data 
obtained by B.U.G. Lab on a sample of JSC Mars 1 analog 
soil provided by E. Guinness. Fig 2a is a close-up of the soil. 
Fig 2b is its spectrum, and Fig 2c is a 3D cylindrical plot of 
the scattering behavior observed for /=45° at 860nm. The z- 
axis is in units of radiance factor. The radial axis is in units 
of emission angle, 0° in the center to 80° around the rim. 
Azimuth (0°) starts at the far right center. The scattering 
plane runs through the middle of the plot, forward scattering 
to the left, backscattering to the right. Note the overall 
backscattering behavior (a lambertian plot would be 
uniformly flat across the top), and the opposition surge. 

Work Completed and in Progress: A number of 

projects have been completed or are underway. 

Complete BRDF characterization of the calibration 
materials to be flown on the MER03 rovers. Samples of 
material and the flight spare set of PAN CAM filters were 
provided by J. Bell (Cornell). Since the PANCAM is fixed 
relative to the calibration target on the rover, the emission 
angle was fixed at 53.5° and incidence and azimuth were 
varied. 

BRDF characterization of Martian analog soils. E. 

Guinness (Wash. U.) provided samples of JSC-Mars-1, and 
R. Morris (JSC) has provided four additional analogs. All of 
these samples were measured using the flight spare 
PANCAM filters and the data will be made available through 
the PDS after a brief analysis. 

Controlled experiments to test photometric models . Well 
characterized particulate samples are being measured in a 
variety of particle sizes, packing densities, and surface 
roughness to test the ability of current photometric models to 
extract surface physical properties. 
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Figure 1. Goniometer Figure 2b. Spectrum of JSC Mars 1 
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Introduction: Photometric observations of plane- 
tary regoliths demand a sophisticated interpreting basis 
that includes theoretical modeling, computer and labo- 
ratory simulations of light scattering. The high-level 
theoretical models were produced by B. Hapke, e.g., 
[1] and they are widely used in planetology. At present 
time quick progress is observed in interpreting pho- 
tometry. We overview here this progress and problems 
in understanding of physical mechanisms forming 
photometric properties of the planetary regoliths. 

Physical mechanisms: The mechanisms important 
for interpreting photometry of the planetary regoliths 
are: shadowing (a) in particulate media and (b) on ran- 
dom topography, as well as (c) incoherent multiple 
scattering and (d) coherent backscatter enhancement. 

(a) Shadowing in particulate media : This is the 
main effect contributing to photometric properties of 
the planetary regoliths at all phase angles even in cases 
when the regolith surfaces are not very absorbing. The 
shadow-hiding effect is studied vciy well [1], never- 
theless, a set of problems related to the effect remain 
unresolved. Recently, for instance, shadow-hiding ef- 
fect was calculated for a particulate structure arranged 
as a fractal that is rather typical fot the planetary sur- 
faces [2]. This development seems to be progressive 
for interpreting photometry as it provides good fits for 
observations and operates with a new parameter, the 
number of hierarchical surface levels that is treated in 
the fractal manner as a fractional value. Note also that 
analytical calculations of the shadowing effect meet 
problems with densely packed particulate media. This 
is accounting for correlation in spatial distribution of 
the medium particles. The most effective way to over- 
come the problem is computer modeling, see, e.g., [3]. 

(b) Shadowing effect on random topography : At 
scales more than particle size, the planetary regoliths 
are complicated with random topography that also in- 
fluences phase curves at large angles. A first approxi- 
mation model accounting for the influence of random 
topography is presented in [1]. Improvements of this 
model are difficult, as the rigorous description of shad- 
owing for single-valued random surfaces demands 
continual integrations that is hardly possible analyti- 
cally [4]. Again the best way to calculate this effect is 
now a ray-tracing modeling [4]. Hie domain, corre- 
sponding to the scale of surface topography is compa- 
rable with size of surface particles, is studied very 
poor. Almost all theoretical models deal with large- 
scale topographies described with gaussian statistics, 


whereas so-called rocky topography being very useful 
for practice has not been studied at all. 

A promising mathematically rigorous approach to 
calculate the influence of surface topography on pho- 
tometric properties of the planetary surfaces can be 
applied for rather dark surface arranged hierarchically, 
when owing to large difference in scales a preceded 
hierarchical level of topography can be considered as a 
reference surface for the next topographical level [5]. 
This approach gives the following photometric disk 
function: D(a,bj)~ cos[7t(/ -a/2)/(7c-a)](cos^)r /(,Hl) , a, /, 
and b being the phase angle, photometric longitude, 
and latitude, respectively. This formula describes very 
well the surface of the Moon and Mercury. 

(c) Incoherent multiple scatter : If a regolith-like 
surface is weakly absorbing, multiple scattering can 
dominate photometric properties in the whole range of 
phase angles. The multiple scattering manifests itself 
inside and between particles as well as between ele- 
ments of surface topography. Multiple scattering inside 
a particle can be considered when size of the particle is 
significantly larger that the wavelength. This forms the 
single-particle albedo and scattering indicatrix, which 
are often treated separately that is somewhat eclectic. 
Moreover, due to finite size of particles the shape of 
the indicatrix depends on the distance between the par- 
ticles, if the distance is small enough. This is important 
for calculations of inter-particle multiple scattering. 

To find phase curves for all orders of scattering, the 
radiative transfer equation is used. However this equa- 
tion does not take into account correlation in ray 
propagation in media, whereas this correlation is im- 
portant as it originates the shadowing effect. Our com- 
puter modeling shows that this factor is rather impor- 
tant even for higher scattering orders [6]. Another gen- 
eral problem of available theoretical models is treating 
particulate media as continual that is bounded with a 
plane, whereas real particulate media are discrete and 
have a transition upper layer, which is formed due to 
vertical fluctuations of particle number near the aver- 
age level and can strongly influence the phase curves 
especially in the first order of scattering. This problem 
is effectively resolved by computer modeling [3,6]. 

(d) Coherent backscatter enhancement is an im- 
portant factor effecting at small phase angles on optical 
properties of regoliths even with low albedo. This en- 
hancement is due to constructive interference of double 
passing light rays along the same trajectories in a par- 
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ticulate medium, e.g., [1,2,7]. The width of the coher- 
ent backscatter opposition spike depends on the char- 
acteristic distance between scatterers. For spikes with 
the width of a few degrees, the distance is of the order 
of 10 wavelengths. If the distance decreases the width 
increases. That is in principle the distance can be so 
small that the coherent backscatter effect will cover 
phase angles more than 20°. Nobody knows is it possi- 
ble or not, as at small distance the ray approach should 
not work. Another question concerns the coherent 
backscatter effect for aggregated particles that are typi- 
cal for the regoliths. What is the minimal size of the 
aggregates, for which the effect is still observable? 
This and many other questions have no answers yet. 

There are several attempts to develop a combined 
model of the shadow-hiding and coherent backscatter 
effects, however, serious progress in this direction has 
not been yet achieved. The problem is that the scatter- 
ers of a medium should be considered simultaneously 
as large enough to cast shadows and small enough to 
form the coherent backscatter enhancement. A pro- 
spective way to develop this problem is to use two 
scales model when a medium is considered to consist 
of large enough (-100 pm) aggregates of small (-1 
pm) graines. This allows one to consider separately the 
shadow-hiding effect (aggregate scale) and the coher- 
ent backscatter (grain scale) [2,7]. 

Laboratory simulations at small phase angles: 
Modeling photometric studies of many different sam- 
ples, structural analogs of planetary regolith, with con- 
trolled parameters were recently carried out in Kharkov 
[7]. One Kharkov’s photometer/polarimeter allows 
measurements at the phase angle range (0.2-17.5°). The 
measurements provide a basis to verify theoretical 
models. For instance it was found that for dielectric 
powders, both the opposition effect and negative po- 
larization become more prominent when the samples 
are compressed. It was also found that the opposition 
effect and negative polarization parameters depend 
strongly on micro- and macro -structured optical inho- 
mogeneities of scattering surfaces. These findings are 
not satisfactory explained by available theoretical 
models of the coherent backscatter enhancement. We 
plan to continue measurements paying attention to 
complicated samples that are made up by different 
materials and surface structures in order to model large 
scale inhomogeneity of planetary surface. 

Geological aspects: A set of images acquired for 
the same scene of a planetary surface at different phase 
angles can be used to determine parameters of model 
photometric function, e.g., phase function slopes for 
each point within the scene. This allows one to use 
images of the surface in the photometric parameter for 


geologic interpretation. For example, phase function 
slopes for a few small lunar regions with Clementine 
UVVis images were mapped in [8]. The maps show 
several anomalous sites, e.g., small fresh impact craters 
have anomalous halos in these maps. The phase func- 
tion within the halos is less steep than for surrounding 
mare surface. We interpret these halos to be due to 
geologically recent impact-induced alteration of the 
equilibrium millimeter-scale regolith structure. This 
structure is established due to micrometeorite bom- 
bardment at geologically short time scale. An anomaly 
of the same signature was found at the Apollo- 15 
landing site. We interpret it as a result of the regolith 
structure alteration with the lander jets. A unique pho- 
tometric anomaly not correlated with albedo was found 
within the Reiner-y formation. This favors the impact 
hypothesis for the nature of the Reiner-y formation. 
Thus mapping the photometric characteristics is a new 
powerful tool in studies of the surface of atmosphere- 
less celestial bodies, e.g., asteroids and Martian sur- 
face. This tool allows studies of recent seismic events 
and population of meteoroids in the inner Solar system. 
The ESA Smart- 1 mission to the Moon will be a good 
opportunity to develop this direction, as during this 
mission many images of some regions of the lunar sur- 
face will be obtained at different phase angles. 

Conclusion: The areas of interpreting photometry 
that seem to be prospective are: 

(a) Accounting for the border effects in studies of the 
shadowing in particulate media; 

(b) Consideration of fractal-like structures to model 
photometric properties of the planetary surfaces; 

(c) Laboratory and computer simulations of photomet- 
ric properties of micro- and macro-inhomogeneous 
particulate surfaces; 

(d) Geological study of the images that present pa- 
rameters of photometric function. 
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Introduction: During our 15 year quest to aug- 
ment vis-near infrared spectroscopic (0.4 - 2.5 pm) 
discoveries with mid-infrared (2.5 - 13.5 pm) spectro- 
scopic discoveries at Mercury, we have achieved much 
and learned even more. This talk (and paper) will de- 
scribe some of those discoveries and lessons learned. 

Methods used to interpret mid-infrared spectra in- 
clude 1) identification of key spectral features diagnos- 
tic of composition. This has been achieved in the labo- 
ratory using terrestrial and lunar rocks, minerals, pow- 
ders, and glasses. 2) comparison of laboratory and 
telescopic mid-infrared spectra of lunar soils from 
similar locations on the Moon. 3) use of the same 
spectrograph to obtain spectra of rocks, minerals, and 
powders to compare to spectra of the Moon’s and Mer- 
cury’s surface and thus calibrate spectrograph perform- 
ance and resulting spectral character. 4) comparison of 
spectra obtained from spacecraft above all of the 
Earth’s atmosphere of objects in the solar system (Jupi- 
ter, Saturn, asteroids) to those obtained from mid- 
infrared instruments from ground-based observatories. 

Near-infrared and mid-infrared diagnostics con- 
trasted: The power of visible and near-infrared spec- 
troscopy (0.4 - 2.5 pm) has l»een demonstrated 
throughout the solar system for decades. For example, 
the size and shape of absorption bands in reflected light 
from the surface caused by crystal field transitions, 
metal-metal intervalence charge transfer transitions, 
and oxygen-metal charge transfer transitions of the 
materials in the regolith, are very diagnostic. The best 
known of these are probably the lunar olivine and py- 
roxene bands which permitted mapping of near-side 
abundance of these minerals and identification of the 
lunar mantle in Copernicus [1]. 

From 2.5 - 7 pm, volume scattering of light from 
regoliths becomes important. There are absorption 
(and emittance) features associated with photons scat- 
tering in individual grains. Many silicates, sulfates, 
and carbonates have diagnostic features in this spectral 
region. Most of this spectral region is available to 
spectrographs at high altitude at ground-based observa- 
tories and from stratospheric observatories like the 
Kuiper Airborne Observatory (KAO— re tired) the 
Stratospheric Observatory for Infrared Astronomy 
(SOFIA-to be commissioned in 2005). 

Major rock-forming minerals have their fundamen- 
tal molecular vibration bands in the region from 7.5 - 
1 1 pm (the Reststrahlen bands). The transparency fea- 
ture between 11 and 13 pm is associated with the 


change from surface to volume scattering. There are 
also major features associated with the bending, twist- 
ing modes of silicates, and other solar system materials 
occurring in the region from 13-40 pm. 

An emissivity maximum (EM) associated with the 
principal Christiansen frequency, usually between 7 
and 9 pm is a good diagnostic of specific mineral iden- 
tity and also bulk regolith or rock type in mixed miner- 
alogic and textural assemblages common in regoliths. 

Two good reviews for more details regarding the 
above highly-condensed and incomplete discussion 
above may be found in [2] and [3]. 

YIS/Near-infrared Discoveries at Mercury: 
Visible and near-infrared (0.4 - 1.5 pm) spectroscopy 
has been successful in documenting the extremely low 
probability of presence of the Fe 2 * charge-transfer ab- 
sorption band in reflected light from Mercury’s sur- 
face. This means that FeO in the regolith is very low in 
abundance, if present at alL [4], [5], [6], [7]. In addi- 
tion, the regolith, while mature, is more transparent 
than the regolith of the Moon, thus indicating low 
abundance of Fe blebs and other opaques that are char- 
acteristic of the lunar mature soils [8], [9]. 

Reflectance spectroscopy is useful on Mercury 
wherever there is reflected light. Of course spatial 
resolution is limited by the spectrograph aperture size, 
the plate scale of the telescope and detector, and 
Earth’s atmospheric turbulance. 

Mid-infrared Discoveries at Mercury: This 

spectral range is dominated by thermal emission. The 
location on the planet, for all data, is biased toward the 
hottest regions in the footprint of the spectrograph ap- 
erture. The other caveats of spatial resolution men- 
tioned above for vis-near-ir observations also hold. In 
addition, spatial resolution is limited by the diffraction 
limit of the telescope. 

Volume Scattering Region : Emissivity peaks at 5.7 
and 6 pm, in a spectrum of Mercury from 100 - 160° 
longitude [10], resemble prominent emissivity maxima 
at 5.7 and 6 pm exhibited by low-iron olivine powders 
in laboratory spectra. The exact wavelength of the pair 
of peaks is seen, in laboratory studies, to shift by small 
amounts according to FeO content in the olivine. 

A strong 5.5 pm emission feature in a spectrum 
from 45-85° longitude closely resembles that of labo- 
ratory c lino-pyroxene powders [11]. The best fit is to 
diopside and the low-FeO abundance indicated by 
near-ir reflectance spectroscopy supports a low-iron 
bearing clino-pyroxene. Spectroscopic observations of 
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Mercury have been made in this spectral region from 
the KAO [10], and between 3 and 7 pm from 13,786 ft 
altitude [11] at Mauna Kea (the NASA Infrared Tele- 
scope Facility, IRTF). 

Emission Maximum in region of Principal 
Christiansen Frequency : Emissivity maxima (EM) at 
or close to 7.9 - 8.0 pm and indicative of intermediate 
silica content (~50 - 57% Si0 2 ) occur in spectra from 
12 - 32° and 22 - 44° [12], 40 - 45°[13], 45 - 85° 
[11], 10 - 75° [10], and 110 - 120° [12] longitude. 
These locations fall in the inter crater plains east of the 
crater Homer. Spectra from 68 - 108° [11] and 100 - 
160° [10] longitude have multiple EM indicating a 
more complicated bulk composition and or mixed min- 
eralogy of more basic composition (45 - 49% Si0 2 ). 

Reststrahlen features: Reststrahlen features in 

spectra from 110 — 120° longitude have been modeled 
with simple linear spectral mixing of laboratory spectra 
from 0 - 74 pm powders. A mixture of labradorite (a 
mineral of the plagioclase feldspar solid solution with 
Na-rich plagioclase, albite— Ab 4 o molar abundance) and 
low-iron ortho-pyroxene powders matches some but 
not all features in the Mercury spectrum. [14]. The 
Reststrahlen features in a spectrum from 68 - 108° 
longitude have not yet been fitted with any model but 
show three emissivity peaks between 7.8 and 9.3 pm, 
some of which may be caused by Reststrahlen features. 
A composite spectrum from CVF spectral imaging at 
205 - 240° longitude [15] indicates a probability of 
picrite based upon comparison to a laboratory spectrum 
of one powdered picrite sample. In this case the sam- 
ple was of a rock type with dominant low-iron olivine 
and other mineralogy. 

Transparency minima: The Mercury spectrum 

from 110 - 120° longitude has a clear and strong 
transparency minimum at 12.3 pm that is at the same 
location as the transparency minimum in a laboratory 
spectrum of labradorite powders [12]. This is consis- 
tent with the location of the EM in the same spectrum 
as described above. Spectra from longitudes centered 
on 80°, 256° and 266° have probable transparency 
minima at 12 pm [16]. The bulk composition associ- 
ated with a transparency feature at this wavelength is 
intermediate to basic (45 - 57% Si0 2 ). Spectra from a 
region centered on 15° has a minimum at 12.5 pm [16] 
indicative of about 44% Si0 2 or an ultra-basic compo- 
sition. A spectrum from a region centered on 229° 
longitude has a doublet transparency minimum with 
one at 12.2 pm and another at 12.5 pm. This is indica- 
tive of more complex mineralogy in that region. 


Summary: About 40% of Mercury's surface has 
been measured spectroscopically. Roughly speaking 
the coverage is of the equatorial and low latitude re- 
gions at most, but not all, longitudes. 

Mercury's surface composition is heterogeneous. 
Regions near Homer and the Murasaki Crater Complex 
appear feldspathic, trending toward Ab^ - Ab^, more 
Na-rich than the lunar anorthosites. Bulk compositions 
are of intermediate silica content. Some mixed compo- 
sitions of more basic silica content are present in the 
regions from 68 - 160° longitude but are not at all lo- 
cations measured. Thus, it appears that low-iron ba- 
salts (mixtures of feldspar, pyroxene, and minor oli- 
vine) and Na-rich feldspathic units dominate the inter- 
crater plains. 

Regions west of Caloris have mixed mineralogy 
and a more complex bulk composition with some basic 
and ultra-basic regolith types. There is evidence for 
picrite-like soils at 205 - 240° longitude. According to 
transparency minima, an ultra-basic composition falls 
east of Homer crater near 15°. Two measurements 
indicate that ultra-basic regolith types are located far 
west of Caloris from about 205° longitude perhaps to 
as far as 15° longitude. 

The picture above is painted in the broadest of 
strokes with the areal extent of the spatial footprint no 
smaller than 200 km by 200 km for the very best spa- 
tially resolved observations [12] and as much as 1000 
km by 1000 km for the least spatially resolved [16]. 
More observations are called for. 
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Introduction: Significant increase of the spatial 
resolution of optical experiments on and near the 
planetary surfaces, e.g., in-situ measurements for Mars 
Pathfinder mission and soon-to-come orbital high 
resolution multispectral observations for Mars Express 
mission, demands efforts in laboratory and theoretical 
simulations of imaging spectrophotometric studies to 
interpret physical and miner a logical characteristics of 
the planetary surfaces. For this reason a new spectral 
imaging facility has been designed and settled since 
April 2000 at the Observatory Midi Pyrenees [l]. This 
is intended to help the definition of new spacebome 
instruments, to simulate planetary observations, and to 
improve the understanding of the effects of observa- 
tional conditions and of physical properties on the bi- 
directional reflectance of rocky surfaces and soils. At- 
tempts to interpret laboratory measurements carried out 
with the facility using the standard photometric model 
by Hapke [2]. Sometimes it may result in non-unique 
solution with ambiguous physical interpretation [3,6]. 
Hapke’ s model implicitly considers that the scale of 
surface topography is much larger than the size of re- 
golith particles handled as a particulate medium. Be- 
sides, this topography is described with the Gaussian 
statistics. However, in laboratory and in-situ planetary 
measurements one deals with surface topographies that 
have the characteristic scale comparable with the size 
of regolith particles, and often such meso-scale topog- 
raphies are formed by rock fragments [1,3, 6, 7] (see in 
this volume an abstract by Cord et al.). To study this 
effect, we carried out a computer modeling of multiple 
light scattering in particulate media complicated by 
rocky topographies and by a single-valued random 
surface with the characteristic scale of order of particle 
size. The geometric optics approximation is used. 

Computer model for particulate media and dif- 
ferent surface topography. We produce a new com- 
puter model for calculations of multiple scattering in 
arbitrary systems of opaque diffusely scattering parti- 
cles. Besides, we carry out simulations of the shadow- 
hiding effect for surfaces with different topographies, 
using this in combination with the model for particulate 
media. Our computer experiments simulate what hap- 
pens in nature or laboratory studies: a system of parti- 
cles generated in a computer memory is “illuminated” 
by rays, which are traced in the system from the light 
source to the observer. Unlike our previous models 
[4,5] the model used here is able to simulate effectively 
multiple scattering in semi-infinite particulate media 
with arbitrary surface topographies having the charac- 


teristic scale comparable with the size of medium par- 
ticles. For the sake of simplicity the particles are sup- 
posed to be spherical with the same size and Lamber- 
tian surface. We use albedo particle surface equal 0.8. 
This allows us to restrict calculations by the first 5 
scattering orders. As in previous studies [4,5], our 
computer modeling consists of two separate stages: (1) 
generation of a random system of particles with certain 
volume density p and (2) Monte Carlo ray tracing. We 
use a system of particles randomly distributed in a rec- 
tangular parallelepiped. The upper border of this par- 
allelepiped (the model frame) can be undulated (or 
complicated) by a topography. The lateral and bottom 
planes are considered to be cyclically closed. That is, if 
a ray leaves the parallelepiped volume, e.g., through 
the lateral plane, it comes into the same volume from 
the opposite side. This is an effective way to simulate a 
semi-infinite medium with a finite volume. 

Results and discussion: Although our model al- 
lows calculations at arbitrary illumination/observation 
geometry we focused here on two of them. We use the 
“mirror” geometry, i = £ = a/2, 5 = it, where i and £ are 
the angles of incidence and emergence, respectively, 6 
is the azimuthal angle, and a is the phase angle. Be- 
sides we use an azimuthal geometry at / = e = 85°, 6 = 
a. We compare three border topographies; (a) flat-in- 
average border of a particulate medium; (b) border a 
particulate medium undulated by a random topography 
with the Gaussian distribution of heights and Gaussian 
correlation function; (c) border of a particulate medium 
complicated by a thin layer of “rocks” randomly dis- 
tributed (with random embedding) over the border of 
the fine-grained particulate medium. 

Before to study media with complicated borders we 
model the shadow-hiding effect for a Gaussian (G) 
topography with RMS slope equals 1 and for rocky 
topography formed with spherical stones randomly 
embedded in a substrate with the flat boundary (Fig. 1). 
Thus the latter is formed by the boundary and the sur- 
face of stones laid upper the boundary. The volume 
density of stone layer is taken 0.3. This density pro- 
vides almost the same number of details in the model 
frame as in case of G-topography. 

Fig. 1 shows phase curves of the probability that a 
ray is not interrupted by G- and rocky (R-) topogra- 
phies at both the mirror (-M) and azimuthal (-A) 
geometries. As one can see at the M-geometry the R- 
topography shows much steeper curve at small phase 
angles than in case of G-relief. This is quite reasonable 
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as R-topographies have angles of surface slopes greater 
than 90°. While at the A-geometry the shadow-hiding 
probability for these so different types of topography 
reveals very similar phase angle behaviors. 

Calculations of phase-angle curves for semi-infinite 
media (with p = 0.1), whose borders are flat-in- 
average (F) and complicated by G- and R- (at p st0 ncs = 
0.3) topographies, are presented in Figs. 2 and 3 for the 
A- and M -geometries, respectively. The main parame- 
ter of G-topography is the RMS slope that is equal to 
the ratio of RMS heights deviation and horizontal cor- 
relation radius L. Another important parameter is the 
mean size of topography details, which is about L. We 
used L = 8r, where r is the medium particle radius. As 
a characteristics of R-relief we use the radius of 
“stones” R = 1 Or. These values are chosen to obtain the 
shadowing properties of R-topography maximally close 
to G-one. We consider also that the thickness of the 
stone layer equals 2R, the stone surface is characterized 
by Lambertian indicatrix, and multiple scattering taken 
place for R-topography and in the basic medium. 
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Figs. 2 and 3 clearly show the opposition effect. The 
brightness surge seen in Fig. 2 is related to correlations 
between incident and emergent rays in both the par- 
ticulate medium and the surface topography at small 


phase angles. We note that for M-geometry (Fig, 3) 
curves F and G are almost coincident. They go sepa- 
rately in the A-geometry (Fig. 2), though at targe phase 
angles they tend to be close to each other. This close- 
ness of the F and G curves is due to that details of the 
G-topography consist of too small numbers of medium 
particles: the details are transparent enough at large 
phase angles in both the A- and M-geometries. At 
small phase angles the G-topography can affect results 
only at A-geometry (as it should be). For comparison 
Figs. 2 and 3 show curves H obtained by multiplication 
of curves F by corresponding shadow-hiding prob- 
abilities shown in Fig. 2. Such a factorized form is 
used, when details of topography are not transparent, 
i.e. when L » r [2]. As one can see at large phase an- 
gles this factorization gives results that are very distinct 
from the cases when multiple scattering between topog- 
raphy elements is taken into account. Curves R in Fig. 
2 and 3 show that the surface particulate topography 
noticeably influences on the photometric properties. 

Conclusion; We showed that particulate media 
with borders complicated by G- and R-topographies 
reveal different photometric behavior in comparison 
with the case of the factorized approximation usually 
used in theoretical models. This points at the need for 
an improved theoretical modeling for describing the 
real case of natural regolithic surfaces, for which the 
meso-scale organization of the textural relief influences 
the surface photometric behavior. 
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Introduction: Over the last several decades, re- 
mote spectroscopic studies have very successfully been 
used to explore the composition of various solid sur- 
faces in the solar system. Building on these initial suc- 
cesses with mineral identification, recent efforts have 
focused on more detailed spectral inferences including 
assessing the composition and abundance of mafic 
minerals. Such efforts are complicated by the fact that 
most silicate-rich surfaces include lithologic combina- 
tions of several mafic minerals which have overlapping 
absorption bands. 

Two of the more quantitative appioaches to spectral 
interpretation are non-linear mixture modeling based 
on Hapke’s radiative transfer theory 1 1] and absorption 
band modeling [2]. Both have been used successfully, 
yet as described below, these two methods are very 
complimentary offer even greater promise if they can 
be successfully combined. 

Absorption Band Modeling: The Modified Gaus- 
sian Model (MGM, [2]) of absorption bands was devel- 
oped to specifically resolve overlapping absorption 
bands in spectra. Under the MGM. spectra are mod- 
eled as a sum of absorption bands (each represented by 
a Modified Gaussian distribution) superimposed onto a 
baseline continuum. The shape of ihese distributions 
was derived theoretically, and empirically shown to 
accurately model isolated absorption bands. Each ab- 
sorption band is characterized by a center, width, and 
strength, and the continuum is described by a slope and 
offset (all as function of energy). Using a simple Beer- 
Lambert approximation for absorbance (e.g., "apparent 
absorbance" [3]), modeling takes place as a function of 
log reflectance and energy. The M( iM has been used 
to characterize continuous changes m individual ab- 
sortpions as a function of composition in several com- 
mon mafic minerals including olivine and pyroxenes 
and successfully used to establish constraints on the 
composition of asteroidal, martian, and lunar surfaces. 

Non-Linear (Hapke) Mixture Modeling: Ab- 

sorption band modeling has many strengths, yet one of 
its drawbacks is that is does not account for albedo. In 
contrast, results from mixture models specifically rely 
on the albedo of various minerals. In addition to al- 
bedo (and single-scattering albedo in particular), the 
Hapke model combines mineral spectra based on radia- 
tive transfer equations that also include surface rough- 
ness and average single particle phase functions. This 
mixture framework allows the addition of mineral 


phases (metal and other opaques) that do not have spe- 
cific absorption features that can be modeled with the 
MGM. In addition, the effects of space weathering can 
be explicitly included within the framework of Hapke 
mixture models [Hapke, 2001]. 

Mixture models also have disadvantages, including 
the fact that they are very dependent on the choice of 
endmember spectral components [4]. Typically, labo- 
ratory analogs (terrestrial or meteoritic samples) are 
used as endmembers in modeling remote spectra. Re- 
cently these "library” materials have also be aug- 
mented with spectra calculated from optical constants 
[5, 6]. However, it is difficult to find endmembers that 
have the same major element, minor element, and 
physical state (grain size, temperature, texture) as those 
observed on a remote surface. 

Integration of the Two Methods: Many of these 
issues in mixture modeling can be overcome by using 
information derived from MGM absorption band mod- 
eling, which derives information directly from the ac- 
tual spectrum as measured. Absorption band and 
Hapke mixture models are therefore quite complimen- 
tary. In particular, results from MGM on mineral com- 
positions can be used to constrain endmember compo- 
sitions and modal abundances for use in Hapke-based 
mixture models. These Hapke models in-tum will pro- 
vide better estimates of abundance than the MGM 
alone. 

In addition, many of the first order approximations 
within the MGM, including a polynomial continuum 
(representing first surface reflectance) and its simple 
absorbance model, can clearly be improved with Hapke 
theory. 

We will explore various approaches to integrating 
these two model and assess their utility by applying 
them to laboratory spectra of controlled mineral mix- 
tures, well characterized natural samples (e.g., meteor- 
ites), and space-altered surfaces (the moon and silicate- 
rich asteroids). By combining the MGM with Hapke 
theory we hope to improve our ability to remotely con- 
strain mineralogy from reflectance spectra. 
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Introduction; It has become accepted lore 
that the myriad of grains of nanophase Fe° 
(abbrev. npFe 0 ) in lunar agglutinates are the 
result of “auto-reduction of impact-melted 
lunar soil in the presence of solar-wind 
hydrogen” [1]. However, recent studies 
have demonstrated other sourc es of npFe 0 in 
lunar soils [2-3], as present in thin patinas 
(~0.1 pm) on the surfaces of most soil 
particles [4]. The major portion of this 
npFe 0 formed by deposition of vapor 
produced by abundant inicrometeorite 
impacts, as documented by the presence of 
multiple and overlapping patinas. A smaller 
portion may have formed by radiation 
sputtering [5]. The vapor deposition of this 
patena npFe 0 is not actually at the surface, 
sensu stricto, but just below [3]. It would 
appear that the presence of npFe 0 in the 
vapor-deposited patinas (rims) on virtually 
all grains of a mature soil [6| provides an 
additional and abundant source for the 
greatly increased 7 s /FeO values observed [7]. 
[At this point, I would like to see just how 
red with fury Bruce’s face must be, as the 
above statements appear to steal his “I told 
you so’s.”.] You see, it was Bruce and his 
Pittsburgh colleagues that had astutely 
predicted all this way back in the Apollo 
days, but its significance went largely over- 
looked for decades. For an elegant and 
thorough review, albeit with a bit of bias, 
please read Hapke [8]. 

Agglutinates by Auto-reduction of FeO: 

The reaction that supposedly occurs 
involves in the auto-reduction paradigm 
involves the presence of solar-wind protons 


during impact-melting of lunar soil — 
FeO ( ]) + 2H + = Fe° = H 2 0. Obviously 
much of the water was lost to space, but 
with the fine-grain size of the resultant np- 
Fe 0 intimately intergrown and dispersed in 
the glass, one would suspect that there 
should be distinct remnants of the w r ater 
trapped within the glass. Yet, detailed FTIR 
studies by Taylor et al. [9] failed to detect 
any traces (>5 ppm) of water in numerous 
agglutinitic glasses. This fact alone should 
make one suspicious of the hydrogen auto- 
reduction theory. 

Surface-Correlated npFe 0 : The Lunar Soil 
Characterization Consortium [e.g., 10] was 
established to coordinate detailed study of 
the physical and chemical properties of the 
finest fractions of numerous lunar soils 
(<45-20, 20-10, and <10 pm). They have 
produced hard data and have come to the 
following conclusions: 1) the abundances of 
agglutinitic glass and the I s /FeO values 
increase with decreasing grain size; 2) if the 
increase in Is/FeO that is attributable to the 
increase in agglutinitic glass is accounted for 
in each change in grain size, the “residual” 
is the possible surface-correlated I s /FeO 
contribution; and 3) on average, there is 
about a 100% increase in I s /FeO value 
between the finest two size fractions, with 
only about 10% increase in agglutinitic 
glass. These general large increases in 
I/FeO versus nominal increases in 
agglutinitic glass contents are also 
demonstrated in the <10 pm fraction of the 
mare soils, as shown below [7]. 
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TABLE 1. Characteristics of Mare Soils. 


Low-Ti Soils 

High-Ti Soils 

Ig/F eO(250|im) 

14 94 

14 

81 

(Maturity) 
IsTF eO(<ionm) 

35 160 

30 

160 

Agglut. Glass 55% 70 

45 

62 

Origin of Agglutinitic Glass: It has been 


demonstrated conclusively by the research 
of the LSCC [7] that it is the fusion of the 
finest fraction (F 3 model; 11] of lunar soil 
that forms the majority of the agglutinitic 
glass. But, it is exactly this <10 p.m fraction 
that also contains the vast majority of the 
surface-correlated npFe 0 . It is considered 
entirely possible that most of the npFe 0 in 
the “agglutinitic glass” had its ultimate 
origin as vapor-deposited npFe 0 in patinas 
on soil particles, which were subsequently 
melted and became part of the glass 
component. Effectively, as the <10 urn 
fraction contains the most vapor-deposited 
npFe 0 , this is melted and incorporated into 
impact glass. The average size of the npFe 0 
in this recycled glass (agglutinate) increases. 
This correlates with the data of Keller et al. 
[12] and James et al [13], where the average 
grain size npFe 0 in “agglutinitic glass” is 
larger than in vapor-deposited npFe 0 . 

New Minerals: Another significance 
of the process of vapor-deposition of npFe 0 
was recently discovered in a regolith breccia 
clast in lunar meteorite Dhofar 280 [14], 
This clast contains numerous, dispersed 
micron size (0.5- 10pm) FeNi metal grains, 
similar to an agglutinate-glass rich soil. In 
addition, grains of three new Fe-Si phases 
are present that vary in size from 2-30pm. 
We have reported, for the first time, the 
occurrence of iron-silicides in lunar rocks. 


These Fe-Si compounds have approximate 
compositions of; FeSi, Fe 2 Siand FeSi 2 . 

A likely possible scenario for the 
formation of these phases involves the 
melting and vaporization of lunar soil by 
micro-meteorite impact. It is thought that 
the ultimate origin of most of this npFe 0 is 
from the vaporization of FeO and Si0 2 from 
the melting of the soil. The FeO dissociates 
to Fe and O and the Fe is condensed onto the 
soil particles. However, no consideration of 
the Si0 2 in the vapor-phase has ever been 
made. We suggest that the Si0 2 in the vapor 
also dissociates to SiO and Si. It is 
thermodynamically possible that the Fe in 
the vapor combined with various 
proportions of Si, also in the vapor and was 
condensed as small Fe-Si grains, in a crude 
way similar to the condensation of Calcium- 
Aluminum-rich Inclusions (CAIs). If these 
Fe-Si phases are the result of impact- 
produced vaporization of lunar soil 
components, it is predicted that many of the 
small, micron-sized Fe metal grains that 
occur in the lunar soil may well contain Si 
and Fe-Si phases as well. 
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Introduction: Spectral mixture analysis has been 
shown to be a powerful, multifaceted tool for the 
analysis of spectral image data [e.g., I, 2]. In this ap- 
proach small areas thought to represent distinct litholo- 
gies are chosen within the image. The spectral proper- 
ties of these representative "endmembers" are then 
combined in a least-squares mixing algorithm to pro- 
vide the best fit for each pixel within the image. Each 
pixel is thus represented by the proportion of each 
endmember required to minimize the difference be- 
tween the model and the actual spectra. Applications 
to both lunar and terrestrial remote sensing include 
mapping soils and exposed bedrock and modeling sur- 
face processes [e.g., 3, 4]. 

Linear vs. Non-Linear Mixing: Depending on the 
application, spectral image data can be “unmixed” as- 
suming linear mixing [5], or non-linew mixing [6, 7]. 
Under the former, the reflectance sped rum of a mix- 
ture is assumed to be a linear combinaiion of the end- 
member reflectance spectra. Non-linear mixing, on the 
other hand, is more appropriate in multiple scattering 
regimes, and each endmember ’s spectral contribution 
to the mixture spectrum is weighted by the endmem- 
ber’s optical extinction coeffient [7]. It has been dem- 
onstrated that a non-linear mixing model more accu- 
rately predicts abundances for the Moon's surface - 
particularly of low albedo materials - and is more ap- 
propriate for most lunar science questions [8]. 

Examples: Recent examples of lunar applications 
of mixture analysis have focused on: ui) merging high 
spectral resolution data with high spatial resolution 
imagery [9] (Figure 1), (b) on estimating and removing 
the contribution of shade/shadow from a scene to allow 
more accurate compositional estimates [10] (Figure 2), 
and (c) on estimating and modeling soil transport 
across mare-highland boundaries [8, 1 1]. In the data 
merging example, a linear mixing model has been ap- 
plied with promising results; future work, however, 
will include nonlinear mixing as part of the data model- 
ing. With shade removal, both linear and non-linear 
models have been applied. The non-! i near model ap- 
pears to better separate low-albedo materials from 
shade, although further validation is needed to assess 
the results more fully. Finally, in modeling mare- 
highland mixing [II], a non-linear model was used to 
estimate the amount of mare and highland soils distrib- 
uted across major contacts, and to develop a physical 
model for transport. 



Figure 1: From [9j Top : Mafic 75pm/.95pm ratio (a 
widely used multispectral surrogate for lpm band depth) of 
Copernicus crater ; using Clementine UWIS camera data. 
Bottom : 1pm band depth image derived from calculated N1R 
data (created by combining multispectral data and tele- 
scopic N1R spectra for the same locations on the Moon via a 
mixture model). Albedo artifacts are eliminated in the band 
depth image. 
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Figure 2: Tog: Original Clementine UWIS camera im- 
age at 0.75 pm of Stevinus Crater. Note shadows within 
craters. Bottom : The same image after shade removal using 
a mixing model. Black spots in small fresh craters are pixels 
that are saturated in the original calibrated data. 


In transforming remote sensing data into quantita- 
tive physical parameters, linear spectral mixture analy- 
sis has proven to be both diverse in its applicability and 
robust. The addition of Hapke’s photometric model to 
allow more realistic separation of intimate mixtures has 
served to make it an even more powerful tool, to which 
we currently owe a portion of our understanding of 
lunar surface composition and evolution. 
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The innermost of Saturn’s classical satellites, Mi- 
mas and Enceladus, have historically proven chal- 
lenging observing targets for Earth -based telescopes 
due to their proximity to the bright planet and rings. 
Since the ring plane crossings in 1995 and 1996, in 
preparation for the upcoming Cassini mission in 2004, 
the Hubble Space Telescope’s Witie-Field Planetary 
Camera (WFPC2) has been monitoring the Saturn 
system at each opposition as well as near quadra- 
ture each year (R. French, PI). This observing pro- 
gram has provided numerous images of Saturn’s clas- 
sical satellites in each of the five WFPC2 wideband 
UVBRI filters (F336W, F439W, F555W, F675W, and 
F814W) and occasionally in the F255W, F785LP, and 
F1042M filters at solar phase angles between 0.26° 
and 6.4°. 

We have produced mutually consistent solar phase 
and rotation curves for Mimas, Enceladus, Tethys, 
Dione, and Rhea and report the first observation of 
the steep, narrow opposition effects of Mimas and 
Enceladus at visible wavelengths, necessitating an in- 
crease in currently accepted values of their geomet- 
ric albedos. Among these satellites, only Rhea was 
known to have a strong opposition effect in the vis- 
ible [1], and 1995 ring plane crossing observations 
revealed Enceladus’ large opposition surge in the 
near-infrared [2]. Combining these new observations 
with existing Voyager and telescopic data enables a 
more precise determination of the physical character- 
istics of regolith particles comprising the surfaces of 
these satellites. Application of the Hapke photomet- 
ric equation [3] to phase curves which include near 
opposition data yields information on the directional 
scattering behavior of surface particles, regolith com- 
paction, and the degree to which coherent backscatter 
plays a role in the scattering of light from the par- 
ticulate surface of an icy satellite [4]. Comparisons 
between the derived photometric parameters for each 
satellite suggest the degree to which location - either 
interior (Mimas and Enceladus) or exterior (Tethys, 
Dione, and Rhea) to the densest point in the E-ring 


- affects the physical properties of regolith particles. 
Mimas and Enceladus have brighter trailing hemi- 
spheres [5-8], while Tethys, Dione, and Rhea display 
the opposite hemispherical albedo dichotomy [5, 8, 

91 - 

Multiple observations during a single HST orbit 
provide a series of measurements of the brightness of 
a satellite at constant phase angle, longitude, and lat- 
itude in each of the five filters, UVBRI, occasionally 
spanning a range from 255 - 1042nm. These spectral 
data indicate the satellite’s color and demonstrate 
wavelength dependence on albedo. 
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Introduction: Thermal infrared data from the Ad- 
vanced Spacebome and Thermal Emission and Reflec- 
tion Radiometer (ASTER) is used to identify the distri- 
bution of three lithologies within the ejecta blanket of 
Meteor Crater, Arizona. This is possible because ther- 
mal infrared (TIR) data have been shown to add 
linearly and thus be easily interpretable [1]. As all 
small, simple impact craters are initially created with 
the same morphology, this has many applications to 
new data from the Thermal Emission Imaging System 
(THEMIS) currently in orbit at Mars and will use the 
premise of Ramsey [2] concerning the use of terrestrial 
multispectral images as proxies for thermal infrared 
data sets of Mars. 

Background: Non-degraded, small impact craters 
outnumber large craters on unmodified planetary sur- 
faces and are of younger age. Younger impact craters 
will have less erosion, a well-preserved ejecta blanket, 
and less dust at the surface, making them excellent sites 
for thermal infrared (TIR) observations of the surface 
composition. Of further interest is being able to iden- 
tify inner-crater stratigraphy and the distribution of 
those units in the ejecta blanket in order to infer the 
geology and local climatic conditions of the region 
before and after crater formation. 

ASTER TIR images are used as a proxy for upcom- 
ing THEMIS data by extracting emissivity spectra of 
the lithologic units representing the Meteor Crater rim 
and ejecta blanket. ASTER TIR has a comparable 
spectral resolution and spatial resolution (90 meters) to 
that of THEMIS TIR (100 meters). In order to test 
methodologies to be used with THEMIS data, image 
end member deconvolutions are performed using 
ASTER TIR spectra of the Meteor ('rater ejecta blan- 
ket. Additionally, hand samples will be collected dur- 
ing field work (summer 2002) and the emission spectra 
of each lithology will be used as sample end members. 
The sample end member analysis will be compared to 
the image end member analysis to validate the image 
end member analysis. 

It should be noted that any comparisons between 
Meteor Crater and a martian impact crater will be 
strictly a comparison concerning morphology and the 
geologic process of impact cratering The sedimentary 


lithologies found at Meteor Crater are not likely to be 
found on the surface of Mars and therefore not in any 
martian crater or ejecta blanket. However, as all small, 
simple impact craters are initially created with the same 
morphology, martian impact craters should show a 
similar pattern of the distribution of ejecta, which will 
be deposited in inverse stratigraphic order with respect 
to the pre-impact stratigraphy. An objective of the 
project is to attempt to quantify this relationship for 
select small, simple impact craters. In addition to 
stratigraphy, other crater properties of Meteor Crater 
will be investigated to determine if these properties can 
be identified by TIR remote sensing. 

It has been shown that an image end member de- 
convolution using the thermal infrared wavelengths can 
be used to identify the lithologies distributed in the 
Meteor Crater ejecta blanket [3]. However, the air- 
borne Thermal Infrared Multispectral Scanner (TIMS) 
and its 4 meter spatial resolution was used in the earlier 
study. Still in question and what this project seeks to 
answer is if THEMIS data sets, with 100 meter spatial 
resolution, can be used to identify the lithologic mem- 
bers of small, simple impact craters on Mars. ASTER, 
rather than TIMS, will be used in this study due to the 
comparable spatial resolution to THEMIS over the TIR 
wavelengths. The ASTER image end member analyses 
will be compared to the earlier TIMS image end mem- 
ber deconvolution for validation. 

Study Site: Meteor Crater, Arizona is located in 
semi-arid north-central Arizona and is estimated to be 
49,000 to 50,000 years old, making it one of the most 
recent and well-preserved impact sites on Earth [4]. 
The simple, bowl-shaped crater is 1 80 meters deep and 
1.2 km in diameter, with an eroded rim standing 30-60 
meters high [5]. The local geology consists of three 
flat-lying sedimentary members with known lithologies 
and contrasting spectral signatures. The oldest unit 
sampled by impact is the Permian Coconino Sandstone. 
Above the Coconino are the Permian Kaibab Lime- 
stone and a thin veneer of the Triassic Moenkopi For- 
mation. The crater rim and ejecta blanket consists of 
the same lithologies as described above, but inverted 
due to impact. Within the ejecta blanket, the units are 
preserved up to 1.2 km (or 2 crater radii) away. Addi- 
tionally, the ejecta blanket has a noticeable wind-streak 
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to the northeast due to climatic conditions since crater 
formation. 

Meteor Crater has been described as a perfect example 
of a well-preserved, simple impact crater by many ge- 
ologists [3, 4, 5, 6, 7]. Meteor Crater is convenient to use 
for remote sensing applications and associated ground- 
truthing with on-site field research. The amount of 
previous work, mineralogical contrast between the 
units [8], low amount of erosion since crater formation 
[3,7] due to a semi-arid environment and young age, 
and the mixing of the ejecta make Meteor Crater the 
logical choice to be an analog for similar craters on the 
martian surface. 

Instrumentation: ASTER is a multispectral 

imager on the Terra spacecraft and has been opera- 
tional since October 2000. ASTER has fourteen bands 
from the visible to the thermal wavelengths, including 
five in the TTR at 90 meter spatial resolution and three 
in the visual and near infrared (VNIR) wavelengths at 
15 meter resolution [9]. ASTER’S TIR bandwidths 
range from 8.13 p.mto 1 1.65 jun. 

THEMIS is a multispectral instrument on the Mars 
Odyssy orbiter and operates over nine bands from 6.5 
pm to 14.5 pm. THEMIS also has five bands in the 
VNIR wavelengths at 20 meter resolution [10]. The 
release of a portion of THEMIS data representing a 
small (3%-4%) fraction of the martian surface is tenta- 
tively set for October 2002. 

Conclusion and Application to Mars: The results 
of the study will lead to a greater understanding of 
THEMIS TIR data of young, small, simple impact cra- 
ters. 

THEMIS TIR data sets of the martian surface will 
yield the emissivity of each pixel at the highest spatial 
resolution ever collected of the martian surface [10]. 
Because a sample end member analysis cannot be per- 
formed, the limits of an image end member analysis 
must be known. An image end member analysis could 
reveal the distribution of various lithologies within the 
ejecta blanket and mineral spectral libraries can be 
used to extract the mineralogical content of each lithol- 
ogy. Impact craters on Mars will be chosen to be ex- 
amined using THEMIS TIR data based on 1) age, as 
young inpact craters will have less erosion and dust, 2) 
size and shape, as this study concerns what can be 
learned from the study of small (<3-4 km diameter), 
simple craters, and 3) the TIR spectra, as the impacted 
area requires spectrally different stratigraphy. 
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Introduction: The linear mixing of thermal infra- 
red (TIR) emission spectra in multi-inineralic mixtures 
has been proven, and its limits and applicability have 
been quantitatively investigated [l,2j. Limiting factors 
in the accuracy of any linear retrieval (spectral decon- 
volution) algorithm include the spectral precision of 
the instrumentation as well as the fact that the number 
of end-members must be < the number of spectral 
bands in the TER [1]. Because of this end-member con- 
straint, there is no way to examine a multispectral TIR 
image using a large, spectral end-member library [3]. 
A possible solution is the implementation of an auto- 
mated, blind end-member algorithm in analyze all pos- 
sible subsets (of arbitrary size k) of minerals present in 
the mixture from within a mineral library (of size n). 
For example, in the only such study to employ this 
technique, the Kelso Dunes, California were examined 
using 375 unique combinations (M) of the most likely 
(n=) 15 minerals present in the dunes [4]. The decon- 
volution model results were analyzed for their “good- 
ness of fit” to laboratory spectra of collected sand sam- 
ples. This investigation proved highly successful, fur- 
ther supporting the capabilities of linear retrieval in the 
TIR [4], However, the process was computationally 
intensive because it was entirely manual. The devel- 
opment of an automated algorithm to accomplish this 
task would decrease the amount of time required for 
such an investigation by many orders of magnitude. 

From a remote sensing perspective, an automated 
blind end-member spectral deconvolution algorithm 
could be useful to determine mineral abundance for 
any mixture (or pixel-by-pixel in a 1 IR image) using a 
large emissivity library of minerals as input. This ap- 
proach is particularly timely as there multispectral TIR 
mapping instruments now orbiting Earth and Mars. 

Algorithm: Combinatorial mathematics states that 
the number of unique k - subsets of n elements is given 
by the equation: 

n\ 

n C k = 

(rt — k)\k\ 

It is observed that as n increases for any given k. the 
result grows significantly. For example, 15 C 4 = 1,365 
and 20 C 4 = 4,845. Algorithms have been written by 
combinatorial mathematicians to generate lists of these 
combinations with different ordering methods. In the 


techniques developed here, a lexicographic (integer) 
ordering algorithm was implemented [5]. 

All routines were written in the MATLAB® envi- 
ronment (Version 6.1.0.450). This software was chosen 
due to its ability to efficiently manipulate large matri- 
ces. The algorithm accepts the following as input: (1) 
an “unknown” mixture spectrum, at high spectral reso- 
lution (518 data points); (2) a library composed of n 
individual mineral spectra, also at high resolution and 
possibly present in the unknown; (3) a user-selected 
maximum output subset size k\ and (4) a user- selected 
threshold abundance below which minerals should not 
appear in the modeled output (e.g. 0-5%). Concerning 
(4), the mathematical techniques utilized in the linear 
deconvolution (and explained in detail in [1]) involve a 
linear regression analysis that can result in negative 
values in the output matrix of end-member fractions. 
Because negative abundances are physically impossi- 
ble, model results containing these are removed from 
further consideration. Linear deconvolution can also 
produce spurious results below the realistic accuracy of 
the model (~5%) [1]. 

A combinations list is then generated for n C 2, and 
each unique combination is used as indices for the 
modeled reconstruction of the unknown emissivity. For 
example, where /i=20, the 190 resultant combinations 
(1,2), (1,3), 1 1,4}, ... {18,20} and {19,20} are each 
used to construct a narrowed end-member library, and 
linear deconvolution is performed on this narrowed 
library, resulting in a vector of abundance fractions for 
each combination. If any fraction is less than the low 
threshold parameter, the combination is removed. Then 
a candidate emissivity is constructed using each frac- 
tion vector, and the root-mean-squared (RMS) error (or 
“goodness of fit”) of the candidate is calculated [1,4]. 
The combination resulting in the lowest RMS error is 
considered the best fitting combination and stored. 
Next, the subset size is increased by one, and a new 
combination list is generated (n C 3 long). If a better- 
fitting combination of three elements is found within 
this list, that combination becomes the new best fit. 
This process is repeated until (k) is reached. Theoreti- 
cally, the combination that results in the lowest RMS 
error will be the one that actually comprises the mix- 
ture. 
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Preliminary Results: Thermal emission spectra of 
several mixtures of known mineralic composition were 
used to initially test the capabilities of this algorithm. 
The mixtures were composed of hornblende, micro- 
cline, oligoclase and quartz, and have been analyzed in 
prior investigations of spectral deconvolution [1]. 

Is RMS error sufficient? To conserve processor 
time and rigorously test the algorithm, a blind end- 
member library of size n=20 was initially used, in con- 
junction with the TES spectral library of 170 minerals 
[3]. An end-member library subset was constructed of 
16 arbitrarily selected mineral spectra and the four 
known end-members. The k value was set at 4 for each 
execution, and the low threshold parameter was set at 
0%. Therefore, the number of combinations examined 
was (20 C 2) + (20 C 3) + (20 C 4) = 6,175. The results 
for three mixtures are tabulated in Table 1. 


Mixture 

2 

3 

4 


Real 


Real 


Real 


Mineral 

% 

Mod. % 

% 

Mod. % 

% 

Mod. % 

Hornblende 

10 

“o!oT 

20 

22.4 

70 

74.4 

Microcline 

10 

12.9 

20 

20.5 

5 

9.1 

Oligoclase 

20 

20.3 

20 

27.8 

5 

0.0 

Quartz 

60 

53.0 

40 

30.2 

20 

15.6 

Pyroxmangite 

0 

14.4 

0 

0.0 

0 

0.0 

Biotite 

0 

0.0 

0 

0.0 

0 

1.6 


Table l . Algorithm output for three mixtures 
using high-resolution thermal emission spectra. 


The modeled mixtures above exhibited the lowest RMS 
error, and were selected as the best-fitting combina- 
tions. The RMS errors for modeled mixtures 2, 3 and 4 
were 0.0040, 0.0039 and 0.0025, respectively. The 
algorithm at times finds other minerals that compose a 
mathematically better fit to the unknown than the ac- 
tual combination of minerals. This suggests that RMS 
error alone may not be capable of judging absolute 
“goodness of fit.’' Also, it is observed that some miner- 
als with abundance <15% are not always detected, or 
other minerals are detected in their place. This was also 
noted in previous investigations and appears to depend 
not only on abundance, but on spectral features as well 
[1,2]. These and other characteristics of the algorithm 
are currently being studied in greater detail. 

Algorithm runtime considerations. Performing 
thousands of spectral deconvolutions per unknown 
spectrum clearly requires a fast and capable computer. 
The trials performed in this study were done on a Win- 
dows® 2000 machine with a 1400 MHz Intel® Pen- 
tium® 4 CPU and 128 MB RAM. The total runtime of 
the algorithm was measured for end-member libraries 
of different sizes, ranging from 4 to 50, and k= 4. In the 


20 end-member library trials (6,175 total combina- 
tions) the algorithm ran for -8 seconds. In a 50 end- 
member library trial (251,125 total combinations) the 
algorithm ran for -4,032 seconds (67 minutes). Further 
investigation into the effects of increasing the number 
of end-members in the spectral library is underway. 
Ideally, the algorithm should be able to choose from 
the entire library of recorded thermal emission spectra 
and quickly produce a resultant fit. 

Effects of degraded spectral resolution. The effect 
of degrading the spectral resolution of the mineral and 
mixture spectra is also of interest, given that the Mars- 
orbiting Thermal Emission Imaging System (THEMIS) 
instrument only has 10 spectral bands in the TER (6.5 
and 14.5 pm). The capabilities of a linear retrieval al- 
gorithm at this resolution were tested by degrading the 
spectral resolution of the end-member and mixture 
spectra to the approximated THEMIS band centers. 
The resulting fractional abundances are tabulated in 
Table 2 (compare to the actual compositions of mix- 
tures 2-4 in Table 1). The RMS errors for modeled 
mixtures 2, 3 and 4 were 0.0041, 0.0033 and 0.0016, 
respectively. Further analysis of the accuracy at low 
spectral resolution is also underway. If successful, this 
algorithm will be implemented for imaging instruments 
such as THEMIS in order to identify mineral percent- 
ages and track sediment transport over time on the sur- 
face of Mars and Earth [3]. 


Mixture 

2 

3 

A 

Mineral 

Mod. % 

Mod. % 

Mod. % 

Hornblende 

5.4 

21.5 

70.0 

Microcline 

10.3 

19.3 

18 4 

Oligoclase 

32.3 

29.2 

0.0 

Quartz 

52.9 

30.7 

119 


Table 2. Spectral deconvolution results 
using 10-band thermal emission spectra. 
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Introduction: The planetary regoliths and dust of 
comets exhibit negative polarization branch (NPB) of 
light scattered at small phase angles. This is demon- 
strated in Figs. 1 and 2, where data for the 1 P/Hal ley 
comet (dark points in Fig. 1) [1] and asteroid 1 Ceres 
(opened circles in Fig. 2) [2] are presented. One can 
see good resemblance of NPB for such different ob- 



jects; in both cases the NPB have a parabolic shape, 
and close values of the inversion angle (near 20°) and 
depth (near 1%). Thus it is reasonable to suppose that 
this resemblance can evidence a common origin of the 
dust and regolith NPB. Nevertheless, here we demon- 



strate that the mechanisms of NPB for comets and 
planetary surfaces can be quite different. In case of 
comets this can be due to scattering by individual par- 
ticles, and for planetary regoliths this could be a com- 
bination of the “comet” mechanism and the coherent 
backscatter effect. 

Light scattering by cometary dust: Individual 
particles with size compared to the wavelength can in 
principle produce NPB. This is well known from calcu- 


lations of scattering properties of spheres and irregular 
particles by, respectively, the Mie formulas [3] and 
DDA approach [4]. Spheres as a model of cometary 
particles are a too rough approximation, since perfect 
balls, even after particle size averaging, can demon- 
strate resonant features (both in intensity and in polari- 
zation) that are not observed for comets. Results of 
DDA modeling for irregular particles are rather contra- 
dictory. In particular, results obtained in [4] and with 
the well-known code [5] show bad coincidence. That is 
why we start DDA modeling using our own codes [6]. 

Here we simulate scattering properties of dust par- 
ticles regarding them to scatter independently of each 
other. Besides we consider that they are randomly 
shaped and have the following mass distribution: m' u , 
where u = 1.5 [7] at a constant density of particles. To 
calculate scattering properties we used an improved 
DDA method presented in [6]. To simulate cometary 
dust particles, a set of 137062 cells is used to form an 
approximately spherical volume. Then two randomly 
chosen cells of the volume are marked as seeds of a 
particle material and two cells are considered as seeds 
of the empty space. Each cell different from the seed 
cells is marked as the nearest seed cell. Examples of 
cross-sections of particles generated are shown on 
Fig. 3. As can be seen we generate particles with very 
irregular shapes that are undoubtedly closer to real 



Fig. 3 

particles of comets than perfect spheres. 

Phase curves of linear polarization degree for a 
cloud consisting of particles with random shapes are 
given in Fig. 1. In this figure curve 1 presents data for 
particles with m - 1.65 + 0.051 at averaging over sizes 
from jc = 0.1 to 5.5 (jc = 2nr/X, where X is the wave- 
length and r is the radius of the sphere with the volume 
that is equal to the particle volume). The chosen refrac- 
tive index can simulate silicate materials with additions 
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of carbon. Curve 2 in Fig. 1 corresponds to dirty icy 
particles (m = 1.31 + 0.0 It); this curve is obtained at 
averaging over x = 0.1 - 10.5. As one can see the 
curves 1 and 2 are strongly distinct from the lP/Halley 
comet data. However, the mixture that is 33% of sili- 
cate and 67% of icy particles produces curve 3 in 
Fig. 1, which is a good model of the comet lP/Halley 
NPB. The mixture reproduces the parabolic shape of 
the NPB, gives almost the same inversion angle and 
rather close value of the NPB depth. The difference of 
values of the latter parameter is perhaps due to that we 
did not take into account in our modeling the contribu- 
tion of the gas component of the 1 P/Halley coma. 

Thus, NPB of comets can be caused with single 
particle scattering by rather small particles ( x< 10), 
size and compositional distribution of which give the 
observed parabolic shape of NPB. This is in consis- 
tence with recent laboratory investigations [8], where it 
was shown that clouds of silicate particles of random 
shapes with size of order of the wavelength do give 
symmetric NPB with the depth near 2% and inversion 
angle of about 20°. At the same time these experiments 
did not show brightness opposition effect that often 
accompanies NPB. 

Our calculations reveal that despite small number 
of large particles with 5 < x < 10, their contribution to 
the NPB is rather high. In particular, NPB of symmet- 
ric shape can be obtained for two component mixtures 
of particles with different sizes, even if each compo- 
nent has significantly asymmetric NPB (see Fig. 1). 

Light scattering by regoliths: The planetary re- 
goliths, for instance the lunar soils, exhibit conspicuous 
NPB that are studied in many works. There have been 
considered many mechanisms of the regolith NPB but 
the coherent backscatter effect is perhaps the most pro- 
spective among them [9], Even for dark surfaces this 
mechanism is valid [10]. We model here the NPB of 
asteroid 1 Ceres, which is a rather dark object. 

We treat the planetary regoliths as a particulate me- 
dia consisting of particles large compared to the wave- 
length. To calculate scattering properties we used an 
approach presented in [11]. This approach allows cal- 
culations of multiple scattering in a particulate medium 
with accounting for the coherent backscatter enhance- 
ment and shadow-hiding effect. 

We consider that the particles of modeling medium 
have the same size and are randomly shaped, like it is 
shown in Fig. 3. They are characterized by m = 1.5 + 
0. If x = 2, and single particle albedo co = 0.66. Note 
that such a particle being isolated does not produce 
NPB at all. Fig. 2 shows results of our simulation for 
the particulate medium with volume density p = 0.3. 
Curves 1 and 2 in this figure present sums of two and 
six orders of scattering, respectively. As can be seen, 


the NPB of the medium in the approximation of sixth 
order is in a good agreement with corresponding data 
for asteroid Ceres 1. This example clearly shows that 
NPB is effectively contributed not only by the second 
scattering order but high orders too. 

We emphasize again that the NPB in this specific 
case is formed only due to the coherent backscatter 
effect and definitely not due to polarization produced 
by individual particles. On the other hand, for planetary 
regoliths the two considered mechanisms can be impor- 
tant. They are able even to give NPB of complicated 
forms, in particular, with double minimum. Such NPB 
were recently discovered for bright asteroids and bright 
planetary satellites [12] and found in laboratory 
measurements [13]. 

Conclusion: Thus, NPB of dust cometary comas 
and planetary regoliths can generally be given rise by 
two different mechanisms. One of them is related to 
single particle scattering. This can be the main cause of 
comets NPB. We find symmetric (parabolic-like) NPB 
to be formed due to superposition of very asymmetric 
NPB corresponding to irregular particles with different 
sizes. The other mechanism is a collective electromag- 
netic effect called the coherent backscatter enhance- 
ment. This mechanism can form NPB independently of 
the fact whether the particles have own NPB or not. 
Unlike the case of scattering by individual particles, the 
coherent backscattering produce also the prominent 
brightness opposition effect. We note that cometary 
comas have actually the brightness opposition surges 
much less than that for planetary regolith. 
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